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1.0  HISTORICAL  INTRODUCTION 


Efforts  to  produce  high  power,  high  quality  electron  beams  for  use  in  Navy  applications 
have  focused  in  recent  years  on  techniques  to  improve  the  current  carrying  capacity  of  the 
betatron  accelerator,  first  operated  by  D.  Kerst1,2  in  1940.  It  was  recognized  early  by 
Kerst3  that  the  current  limit  in  a  betatron  is  due  to  the  defocusing  effects  of  space  charge 
at  injection  time;  space  charge  effects  become  rapidly  less  important  as  beam  energy  is 
increased.  One  obvious  way,  therefore,  to  increase  the  current  is  to  increase  the  energy  of 
the  injected  beam.  Of  course,  there  are  limits  to  this  strategy;  the  injector  rapidly  becomes 
large  and  expensive  as  the  energy  is  increased.  Still,  this  technique  has  been  successfully 
used  to  produce  high  current  (~  100  A)  beams  in  a  betatron.4,5  By  contrast,  the  largest 
Kerst  betatron6  achieved  currents  of  about  100  mA  using  a  100  keV  injector. 


Besides  increasing  the  energy  of  injection,  other  methods  have  been  tried  to  improve 
the  current  limit,  including  increasing  the  external  focusing  strength  in  various  ways  and 
by  decreasing  the  defocusing  effect  of  space  charge  by  the  introduction  of  a  background 
plasma.  Devices  using  this  latter  approach,  called  plasma  betatrons,7  never  worked  very 
well,  even  with  the  addition  of  an  applied  toroidal  field.8 


Early  attempts  to  increase  the  focusing  strength  of  applied  magnetic  fields  had  mixed 
results.  Among  the  earliest  were  the  fixed  field,  alternating  gradient  (FFAG)  betatron,9 
based  upon  the  (then)  newly  discovered  strong  focusing  principle10  and  the  betatron  with 
added  longitudinal  field.11 

Two  variations  of  the  FFAG  betatron  were  designed,  a  radial  sector  type  and  a  spiral 
sector  type;3  both  use  fields  that  are  fixed  in  time,  allowing  particle  orbits  to  move  radially 
as  energy  is  increased.  The  current  in  these  devices  is  thought  to  be  limited,  as  in  other 
strong  focusing  machines,  by  the  tune  shift  due  to  space  charge.  The  magnet  design  in 
these  devices  is  significantly  more  complex  than  for  a  conventional,  weak  focusing  design.3 

The  addition  of  a  toroidal  field  to  a  betatron  appears  to  have  been  first  tried  in 
England  shortly  after  World  War  II,  at  the  wartime  radar  establishment  at  Malvern.11  It 
was  found  that  such  an  additional  field  decreased  the  current  outpi  t.  of  *he  device,  the 
reason  being  that  the  electrons  tended  to  intercept  the  injector  after  o.  or  a  few  turns 
around  the  machine.  This  early  experiment  was  soon  abandoned  but  the  basic  idea,  that  a 
longitudinal  field  could  be  used  to  control  space  charge  defocusing,  was  nonetheless  sound; 
the  fundamental  problem  was  to  inject  and  trap  a  high  current  beam  in  the  presence  of  a 
longitudinal  field. 

In  1972  Rostoker12  described  theoretically  how  to  form  a  high  density  beam  of  elec¬ 
trons  in  a  torus,  using  the  inductive  charging  technique  pioneered  on  the  HIPAC  device.13 


1 


A  few  years  later  he  proposed14  accelerating  the  resulting  high  density  beam  using  a  flux 
changing  coil,  as  in  a  conventional  betatron.  A  similar  ‘modified  betatron’  for  high  current 
proton  beams  (not  using  inductive  charging)  was  proposed  by  Sprangle  and  Kapetanakos.15 
On  the  basis  of  this  theoretical  work  and  some  additional  studies16,17  two  experimental 
programs  on  the  modified  betatron  were  begun,  one  at  the  University  of  California  at 
Irvine  and  another  at  the  Naval  Research  Laboratory. 

The  program  at  Irvine  was  an  outgrowth  of  some  earlier  work  on  electron  cloud  trap¬ 
ping  at  Maxwell  Laboratories.18  It  was  found  that  a  large  fraction  of  the  cloud  injected 
by  inductive  charging  could  not  be  accelerated  by  the  betatron  field;  the  suggestion  was 
that  electrons  were  trapped  in  field  bumps  produced  by  the  discrete  coils  used  to  make 
the  toroidal  field.  By  using  a  variant  of  inductive  charging  the  Irvine  group  was  able  to 
produce  currents  up  to  200  A  and  energies  up  to  1  MeV  (but  not  both,  simultaneously)  in 
a  modified  betatron.19 

At  NRL,  extensive  analyses  were  made  of  injection,20,21,22  beam  equilibiia23  and 
stability,24,25,26  and  extraction27  during  the  design  and  construction  of  a  larger  scale  modi¬ 
fied  betatron,  designed  to  produce  a  50  MeV,  10  kA  electron  beam.  The  technical  obstacles 
to  injection  and  trapping  and  to  chamber  design  have  proven  to  be  difficult  to  overcome, 
however.  Recently,  the  successful  injection  of  a  2-3  kA  beam  has  been  achieved  using  an 
internal  diode  injector  arrangement.22 

One  of  the  problems  with  both  the  conventional  and  modified  betatrons  is  their  sensi¬ 
tivity  of  the  radial  location  of  the  equilibrium  orbit  to  the  particle  energy,  a  feature  which 
complicates  the  injection  and  trapping  of  large  currents.  It  was  first  recognized  by  Mon- 
delli  and  Roberson,28  in  the  course  of  their  study  of  the  Racetrack  Induction  Accelerator, 
that  the  use  of  an  t— 2  stellarator  field  will  significantly  reduce  this  sensitivity.  It  was 
subsequently  suggested29  that  the  addition  of  an  t= 2  stellarator  field  to  a  modified  beta¬ 
tron  could  improve  its  current  carrying  capacity.  In  a  test  of  this  prediction,  stellarator 
windings  were  added  to  the  betatron  at  Irvine.  A  200  A  beam  was  accelerated  to  2  MeV 
in  the  device;30  smaller  currents  were  subsequently  accelerated  to  10  MeV. 

The  present  report  presents  theoretical  result^  produced  with  ONR  support  at  SAIC, 
>y  on  the  stability  and  scaling  of  compact  electron  accelerators  using  continuous  strong  focus¬ 
ing  in  the  form  of  stellarator  windings.  The  work  reported  here  has  been  mainly  analytical 
in  nature  but  has  also  involved  the  development  and  use  of  a  variety  of  computer  programs 
needed  to  assist  various  aspects  of  the  study.  (See  Appendix  C.)  During  the  contract  we 
have' attempted  to  address  in  some  considerable  detail  every  major  physics  issue  which 
would  affect  the  design  of  accelerators  using  stellarator  fields.  Generally  speaking, ’these 
issues  divide  into  two  breadj  categories:  beam  equilibria  and  matching  and  beam  stability; 
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each  of  these  may  be  further  sub-divided  into  single  particle  and  collective  phenomena. 

Beam  equilibria  are  states  in  which  focusing  and  defocusing  forces  are  in  balance. 
Collections  of  non-interacting  single  particles,  launched  with  a  distribution  of  initial  con¬ 
ditions,  are  defocused  by  the  finite  beam  emittance  and  focused  by  external  fields.  For 
strong  focusing  systems  there  typically  exist  regions  of  parameter  space  in  which  particle 
orbits  are  stable,  separated  by  so-called  stopbands,  in  which  the  motion  is  unstable.  For 
non-interacting  particles,  every  collection  of  initial  conditions  in  a  stable  system  represents 
an  equilibrium.  It  is  possible  to  construct  special  equilibria,  so-called  matched  solutions, 
whose  macroscopic  properties  (e.g.  beam  shape)  repeat  periodically  with  the  same  pe¬ 
riod  as  the  focusing  system;  these  have  been  constructed  generally  for  1-D  systems  by 
Courant  and  Snyder.31  In  the  case  of  realistic  transport  systems  in  which  the  focusing 
system  changes  qualitatively  from  section  to  section,  numerical  analysis  is  necessary  in 
order  to  analyze  the  change  in  beam  properties.  A  single  particle  code,  SPIRAL,  has  been 
constructed  for  this  purpose.  This  code  and  its  use  in  performing  beam  matching  calcu¬ 
lations  from  a  diode,  through  a  focusing  (compression)  section  and  into  £  =  2  stellarator 
transport  section  are  described  below. 

Once  space  charge  forces  are  considered,  the  situation  becomes  more  complex;  only 
in  the  simple  case  of  a  uniform  elliptical  beam  can  equilibria  be  constructed  analytically. 
In  the  case  that  the  two  transverse  degrees  of  freedom  are  decoupled  this  was  shown  by 
Kapchinskii  and  Vladimirskii;32  the  resulting  K-V  equilibria  are  represented  by  the  usual 
K-V  envelope  equations.  The  stability  of  these  envelopes  in  a  bumpy-torus  (£=0)  stellatron, 
has  been  studied  under  this  program.33  ^  ^^rr~ 

When  the  transverse  degrees  of  freedom  of  particle  motion  are  coupled,  as  in  the  £=2 
stellatron,  the  analysis  of  beam  equilibria  and  stability  becomes  more  complex. 
Gluckstern34  appears  to  have  been  the  first  to  construct  a  K-V  type  of  equilibrium  for 
a  continuously  rotated  quadrupole  field;  Chernin35  carried  out  the  analysis  including  a 
longitudinal  field.  There  remained,  however,  the  question  of  how  small  departures  from 
these  equilibria  behaved.  To  answer  this  question,  an  envelope  type  of  equation  govern¬ 
ing  these  systems  was  required.  Such  an  equation  was  developed  and  the  stability  of  the 
K-V  £—2  stellatron  distribution  was  studied  under  this  program.36  A  code,  XYENV,  was 
written  to  integrate  the  generalized  envelope  equation  in  an  arbitrary  focusing  system, 
including  the  effect  of  space  charge  for  a  uniform  elliptical  beam.  It  has  been  used  to 
study  the  matching  problem  from  longitudinal  to  stellarator  fields  in  the  presence  of  space 
charge  forces  and  will  continue  to  be  used  as  a  design  tool  for  transport  line  design. 

Beam  stability  issues,  like  the  issue  of  envelope  stability,  naturally  divide  themselves 
into  single  particle  and  collective  effects.  The  most  basic  stability  question  for  the  stella- 


tron,  that  of  the  nature  of  the  stopband  for  single  particle  motion,  was  answered  in  the 
original  paper  of  Roberson,  Mondelli  and  Chernin.29  Other  single  particle  instabilities  are 
inevitably  present  in  real  machines,  however,  in  the  form  of  orbital  resonances  due  to  the 
presence  of  field  errors.  Both  linear  and  non-linear  resonances  have  been  investigated  in 
some  detail;37,38  our  conclusion  is  that  they  are  potentially  rather  serious  for  the  stellatron 
and  may  act  to  limit  the  achievable  current  and  energy  of  the  device.  Some  evidence  of 
the  ability  to  cross  high  order  resonances  has  been  reported  by  the  Irvine  group.39  One 
interesting  theoretical  finding  is  the  tendency  of  a  strong  longitudinal  field  to  reduce  the 
tune  shift  due  to  space  charge. 

Collective  beam  modes  have  also  been  investigated  during  this  effort.  In  particular, 
an  extensive  study  was  reported40  on  the  negative  mass  instability  in  the  stellatron,  and  its 
coupling  to  transverse  modes.  An  eigenvalue  problem  for  the  growth  rate  was  formulated 
and  solved  numerically;  in  addition,  the  eigenvalue  problem  was  solved  approximately,  in 
the  smooth  approximation,  and  the  resulting  dispersion  relation  reproduced  the  numerical 
results  to  very  good  accuracy.  It  was  shown  that  the  correct  dispersion  relation  was  just 
that  for  the  modified  betatron,26  under  the  replacement  of  a  single  term,  to  include  the 
smoothed  effect  of  the  stellator  field. 

Another  collective  mode  of  a  beam  in  a  stellarator  field  was  identified  by  Hughes 
and  Godfrey.41  This  is  a  short  wavelength,  transverse  electromagnetic  mode  involving  a 
resonance  between  the  betatron  oscillations  of  the  beam  and  a  vacuum  waveguide  mode. 
Linear  growth  rates  can  be  significant  under  some  circumstances  and  computer  simulations 
have  shown41  that  this  instability  can  be  disruptive.  There  is  some  evidence  that  self 
beam  energy  spread  due  to  space  charge  will  help  stabilize  this  mode.  Avoidance  of  this 
instability  requires  careful  design  of  the  beam  transport  line. 

During  the  course  of  this  work,  it  became  clear  that  recirculating  induction  linacs 
hold  promise  for  solving  some  of  the  problems  facing  a  betatron,  even  one  using  strong 
focusing:  for  the  stellatron,  as  for  the  modified  betatron,  injection  and  extraction  remain 
key  problems.  Also,  the  relatively  slow  acceleration  cycle  of  a  betatron  allow  time  for  the 
growth  of  beam  instabilities.  Finally,  the  weight  of  a  high  power  betatron  and  associated 
power  supplies  and  shielding  looked  formidable. 

Based  on  the  early  theoretical  work  of  Mondelli  and  Roberson28  on  a  racetrack  re¬ 
circulator  using  stellarator  fields,  on  the  experimental  efforts  on  recirculation  of  Wilson,42 
and  on  the  work  of  Birx43  on  branched  magnetics,  a  spiral  line  inductor  accelerator  (SLIA) 
was  proposed44  in  1985.  The  SLIA  uses  a  continuous  beam  line  which  passes  through  an 
accelerator  module  multiple  times.  This  design  vastly  reduces  the  difficulty  of  injection  and 
extraction  and,  since  the  effective  accelerating  gradient  is  much  larger  than  in  a  betatron, 
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instabilities  have  much  less  time  to  grow;  one  finds,  for  example,  that  the  long  wavelength 
negative  mass  instability  in  the  bends  is  not  nearly  as  important  as  in  a  betatron.  Also  in 
the  SLIA,  the  weight  of  the  power  supplies  is  reduced  by  resetting  and  re-using  the  induc¬ 
tion  module  cores  between  beam  passes.  Finally,  the  use  of  static  fields  in  the  transport 
line  allow  very  low  field  errors  to  be  achieved. 

While  solving  many  problems  of  the  closed  orbit,  betatron- type  accelerators  the  SLIA 
introduces  at  least  one  important  new  physics/design  issue:  the  beam  break-up  instability. 
This  instability  is  due  to  the  off-axis  passage  of  a  beam  through  an  accelerating  module 
which  excites  an  electromagnetic  cavity  mode  which,  in  turn,  further  deflects  the  beam. 
It  has  proven  to  be  the  single  most  important  current  limiting  instability  in  the  ATA 
accelerator.45  Techniques  for  its  control,  including  careful  cavity  design  and  strong  focusing 
are  under  investigation.  SAIC  has  developed  a  numerical  code,  BBUSH,  which  tracks  a 
beam  through  an  arbitrary  configuration  of  gaps  and  transport  lines.  It  was  used  in  our 
analysis  of  the  effect  of  sextupole  windings  on  stabilization  of  the  bbu.46 

The  Spiral  Line  Induction  Accelerator  was  favorably  reviewed  by  the  JASON  commit¬ 
tee  in  June  1986  and  was  very  highly  rated  for  Navy  and  other  applications  in  the  extensive, 
DARPA/Services  sponsored  review,  the  Net  Technical  Assessment  of  the  Tactical  Appli¬ 
cations  of  Relativistic  Charged  Particle  Beam  Weapons,  conducted  during  1986-87. 

The  remainder  of  this  report  provides  technical  details  on  the  work  performed  at 
SAIC  under  ONR  support.  Three  main  sections,  on  beam  equilibria  and  matching,  on 
beam  stability  and  on  the  Spiral  Line  Induction  Accelerator  follow.  Each  section  includes 
a  discussion  of  work  performed  and  summarizes  our  results.  Some  details  are  given  in 
Appendices:  Appendix  A  includes  reprints  of  published  papers,  Appendix  B  is  a  list  of 
conference  presentations  and  Appendix  C  is  a  list  of  brief  description  of  computer  codes 
written  in  the  course  of  this  work.  A  final  section  gives  our  conclusions  and  recommenda¬ 
tions. 
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2.0  TECHNICAL  ACCOMPLISHMENTS 


This  section  describes  the  technical  work  done  with  ONR  support  on  the  physics 
and  design  issues  affecting  high  current  electron  accelerators  using  stellarator  focusing. 
Two  broad  areas  are  discussed:  beam  equilibria  and  matching  (Section  2.1),  and  beam 
stability  (Section  2.2).  A  third  section  (2.3)  describes  how  our  studies  of  these  issues  led 
to  the  proposal  of  the  Spiral  Line  Induction  Accelerator,  now  presently  under  assessment 
for  Naval  tactical  requirements.  Each  section  and  subsection  contains  a  summary  of  our 
approach  to  a  problem  and  our  main  results  and  conclusions.  Some  details  are  provided 
in  Appendices. 

2.1  Beam  Equilibria  and  Matching 

In  a  periodic,  strong  focusing  transport  system  a  matched  beam  is  defined  as  one 
whose  macroscopic  properties  (beam  radii  and  orientation)  repeat  with  the  period  of  the 
focusing  system;  matched  beams  are  typically  defined  to  be  monoenergetic.  While  never 
precisely  realized  in  practice,  the  nature  of  a  matched  solution  for  a  particular  focusing 
system  is  important  to  know  so  that  the  dependence  of  the  beam  size  on  its  current,  energy 
emittance  and  on  the  applied  focusing  strength  can  be  estimated.  One  would  also  like  to 
know  the  form  of  the  matched  solution  so  that  a  beam  may  be  prepared,  as  much  as  is 
practical,  in  a  matched  state  before  launching  into  the  transport  system. 

2.1.1  Self-Consistent  Treatment  of  Space  Charge 

For  a  continuously  twisted  quadrupole,  with  no  added  longitudinal  field,  Gluckstern34 
has  described  K-V  type  equilibria;  the  form  of  these  equilibria  when  a  longitudinal  field  is 
present  has  been  worked  out  by  Chernin.35  The  approach  is  to  identify  constants  of  motion 
for  the  single  particle  equations  of  motion,  assuming  a  uniformly  populated,  elliptical 
beam  and  then  to  find  a  function  of  these  constants-a  solution  of  Vlasov’s  equation-which 
reproduces  the  uniform  beam.  The  correct  function  of  constants  of  motion  turns  out  to  be 
a  delta  function  of  a  particular  argument,  in  close  analogy  to  the  distribution  function  used 
by  Kapchinsky  and  Vladimirsky32  in  their  now  classic  work.  Once  the  distribution  function 
is  found,  it  is  straightforward  to  calculate  beam  radii  as  functions  of  system  parameters, 
assuming  certain  emittances,  or  phase-space  areas,  are  known.  For  non- zero  current  the 
functions  are  implicit,  due  to  the  dependence  of  the  single  particle  betatron  frequencies 
on  the  current;  a  numerical  technique  to  solve  for  the  equilibrium  beam  radii  as  functions 
of  various  system  parameters  was  developed.  Details  of  the  calculation  with  numerical 
examples,  are  given  in  the  paper  of  reference  (36),  included  in  Appendix  A. 

In  order  to  study  the  behavior  of  beams  away  from  equilibrium,  it  is  necessary  to  have 
a  formalism  which  describes  a  beam  in  coupled  focusing  systems,  like  the  stellarator.  It 
is  possible  to  formulate  a  general  envelope  equation  for  a  coupled  system,  including  space 
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charge  for  a  uniform  beai  jy  using  the  notion  of  rms  quantities.  This  was  carried  out 
during  the  present  program  and  is  described  in  the  following  section.  A  full  accounting  is 
given  in  reference  (36). 

2.1.2  General  Formulation  of  the  Envelope  Equation 

When  the  two  transverse  degrees  of  freedom,  x  and  y,  of  particle  motion  are  linearly 
coupled,  the  K-V  envelope  equations  no  longer  apply  to  the  beam.  In  some  simple  cases, 
like  solenoid  focusing,  it  is  possible  to  transform  to  a  rotating  frame  (e.g.  the  Larmor  frame) 
in  which  the  single  particle  equations  of  motion  decouple,  but  in  general  the  problem  is 
most  directly  treated,  and  interpreted,  in  the  laboratory  frame. 

The  derivation  of  the  appropriate  envelope  equation  for  coupled  systems  proceeds  in 
a  different  manner  than  the  classic  derivations  of  the  K-V  equations.32,47  One  begins  with 
the  paraxial  single  particle  equations  of  motion,  written  as 

t/(«)  =  M(s)v(s)  (2-1) 

where  u(s)  is  a  4  x  1  column  vector  with  elements  (x,  x' ,  y,  y')  and  M(s )  is  a  4  x  4 
matrix  describing  the  focusing  system;  s  is  the  path  length  variable  and  a  prime  mark 
denotes  d/ds.  If  the  beam  cross  section  is  a  uniformly  filled  ellipse,  then  the  space  charge 
electrostatic  and  magnetostatic  forces  are  linear  and  a  space  charge  term  may  be  included 
in  M(s);  the  form  of  this  term  is  given  in  reference  (36). 

In  the  paraxial  approximation,  all  particles  in  a  monoenergetic  beam  have  the  same 
value  of  longitudinal  velocity,  i,  therefore  all  particles  in  a  particular  “slice”  of  the  beam 
remain  in  that  slice.  To  proceed  to  an  envelope  equation  which  describes  the  evolution 
of  the  rms  size  of  a  slice  as  it  propagates  through  a  transport  system,  one  defines  the 
quantities 


Eij(s)  =<  vi(s)vj(s)  >  -  <  Vj(s)  ><  Vj(s )  >  (2-2) 

for  i,j  =  1, 2, 3,4  where  the  angle  brackets  denote  an  average  over  particles  in  a  slice.  The 
beam  covariance  matrix  £  satisfies 


£'(s)  =  M(s)£(s)  +  (M(s)£(s))T  (2-3) 

by  virtue  of  Eq.  (2-1);  in  Eq.  (2-3),  T  denotes  transpose.  Equation  (2-3)  is  the  gener¬ 
alization  of  the  K-V  envelope  equations,  the  purely  spatial  components  £n,  £33,  and 
£13  =  E31  defining  the  axes  and  orientation  of  the  rms  beam  ellipse.  Given  an  initial  value 
for  the  beam  matrix,  Eq.  (2-3)  allows  it  to  be  calculated  at  any  point  downstream  from  the 


launching  point.  In  practice  such  an  initial  value  may  be  easily  calculated  by  an  electron 
gun  code,  at  the  entrance  to  the  transport  system,  for  example. 

To  recover  the  K-V  envelope  equations  for  decoupled  2x2  systems,  we  begin  by  noting 
that  M  and  E  are  of  2  x  2  block  diagonal  form  for  decoupled  systems,  assuming  E(0)  is. 
Defining  azz  =  Sn,  ox'x'  =  E22,  and  crZX‘  =  E12  =  E21  the  equations  for  the  x-coordinate 
may  be  written  out  using  Eq.  (2-3)  as: 


a  zx  —  2  <JXX>  (2  —  4) 

<J xxt  —  <tZ'z'  "b  (  kz  -f-  IQiz^zx  (2  —  5) 

K'z'  =2(~kl+ hxz)<7zz‘  (2-6) 

where  k2  =  — M21  contains  the  external  focusing  strength,  I  =  u/(/32-/3),  v  =  Budker’s 
parameter,  and 


Qzx  — 


1/2/ 
(Tzx  {<T 


1/2  4.  al/2\ 
xz  1  Oyy  ) 


(2-7) 


(The  space  charge  term,  derived  in  reference  (36),  has  been  included  in  Eqs.  (2-4,  5,  6).) 
Equations  (2-4)-(2-6)  are  coupled  non-linear  equations  but  they  do  admit  a  simple  constant 
of  motion,  the  un-normalized  emittance  ex  defined  by 


€i  =  1 6[<r,*<7t.«»  -  <t2zz.}  .  (2-8) 

Using  Eq.  (2-4)  to  eliminate  azx<  and  Eq.  (2-8)  to  eliminate  ax>x'  one  may  find  a  single 
equation  for  axx.  If  we  define 


rx  =  2 ali2  ,  (2-9) 

the  factor  2  arising  from  the  relation  between  a  radius  and  an  rms  radius  of  an  ellipse,  the 
resulting  equation  for  rx  is 


r"  +  k 


x'  X 


4/ 

r-z 


(2  -  10) 


which,  along  with  the  analogous  equation  for  ry,  are  the  K-V  envelope  equations. 

When  coupling  between  x  and  y  motion  is  present,  it  is  necessary  in  general  to  integrate 
Eq.  (2-3)  numerically.  A  computer  program,  XYENV,  has  been  written  for  this  purpose. 
The  program  allows  a  completely  general  specification  of  the  paraxial  fields  as  functions 


of  axial  distance.  It  has  been  used  to  study  the  envelope  oscillations  in  the  stellatron  and 
beam  matching  problems  in  the  SLIA.  As  an  example  of  the  latter,  we  consider  the  behavior 
of  a  beam,  initially  matched  in  a  uniform  solenoid,  as  it  propagates  into  a  stellarator  field. 

Figure  2-1  shows  plots  of  the  majr  and  minor  beam  radii  versus  distance  (z)  along 
the  system  for  different  values  of  beam  current.  The  beam  starts  out  at  z  =  —20  cm 
matched  in  a  2.98  kG  longitudinal  field;  this  value  is  chosen  so  that  the  1  MeV  beam  has  a 
20  cm  Larmor  wavelength.  The  beam  emittance  is  taken  to  be  en  =  0.3 IkA  ,  except  for  the 
zero  current  case  for  which  we  have  taken  e  =  0.1  rad-cm.  A  stellarator  field  with  on- axis 
gradient  of  500  G /cm  and  period  9  cm  is  slowly  turned  on  using  a  “tanh”  function,  in  such 
a  way  that  Maxwell’s  equations  are  always  satisfied  in  linear  order;  at  2  =  0  the  stellarator 
field  has  reached  half  strength.  One  sees  in  the  figure  the  development  of  oscillations  in 
the  major  and  minor  beam  radii;  the  amplitude  of  the  oscillations  is  seen  to  be  fairly 
insensitive  to  beam  current  in  this  example.  In  other  cases  one  finds  that  certain  modes 
of  envelope  oscillation  are  linearly  unstable,  leading  to  large  oscillation  amplitudes.  This 
situation  is  analyzed  in  reference  (36).  (See  Appendix  A.) 


2.1.3  Beam  Matching:  The  SPIRAL  Code 

As  mentioned  earlier,  special  equilibria  exist  called  matched  solutions  in  which  macro¬ 
scopic  properties  repeat  periodically  with  the  period  of  the  focusing  system.  For  an  actual 
-  as  opposed  to  ideal  -  magnetic  field  system,  there  exist  coil  windings  which  include  feeds 
and  coil  terminations  as  well  as  winding  errors  and  physical  winding  constraints.  To  design 
a  magnetic  focusing  system  with  these  limitations  that  is  matched  to  the  desired  beam 
properties  of  a  device  requires  certain  numerical  capabilities.  Many  of  these  capabilities 
exist  in  the  code  SPIRAL. 

The  SPIRAL  code  has  been  developed  and  is  used  as  a  magnetic  field  coil  design  tool. 
It  computes  the  magnetic  focusing  fields  from  actual  coil  geometries,  including  feeds  and 
terminations.  SPIRAL  is  a  fully  three-dimensional  code  which  calculates  the  magnetic 
fields  using  a  Boit-Savart  solver  and  has  the  capability  to  track  single-particle  trajectories 
through  magnetostatic  fields.  The  code  has  many  diagnostics,  one  being  the  ability  to  map 
out  the  magnetic  field  structure  in  several  ways,  and  another  being  the  ability  to  track  the 
evolution  of  beam  emittance  and  particle  phase-space  through  a  device.  SPIRAL  can  also 
initiate  several  distribution  functions,  including  the  K-V  distribution  as  well  as  a  simple 
uniform  beam  distribution  with  specified  emittance  aniJ  divergence. 

As  an  example,  consider  the  focusing  system  in  the  SLIA  experiment  in  which  a  beam 
is  to  be  matched  from  a  diode  source,  through  a  compression  region  and  into  a  combination 
solenoid/stellarator  transport  section  (see  Fig.  2-2).  The  stellarator  in  this  system  is  an 
£  —  2  (quadrupole)  winding.  Since  the  stellarator  fields  are  to  be  turned  either  on  or  off 
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Figure  2-2  Magnetic  field  system  configuration. 


in  the  SLIA,  the  system  was  designed  such  that  the  compression  region  was  to  match  the 
beam  (energy  0.6  MeV  and  current  800  A)  into  the  solenoidal  field  only.  The  stellarator 
fields  and  any  compensation  magnets  must  work  independently  and  provide  minimal  beam 
disruption.  In  addition,  it  is  desired  that  the  stellarator  coils  begin  at  roughly  the  same 
axial  position  as  the  solenoidal  windings. 

Using  SPIRAL,  it  was  found  that  shifting  the  stellarator  windings  to  begin  about  8.2 
cm  further  down  axially  and  putting  a  discrete  quadrupole  magnet  centered  at  the  opening 
of  the  stellarator  coils  at  a  relative  angle  of  45°  (see  Table  2-1)  provided  a  reasonable 
match  with  minimal  beam  disturbance.  Plots  of  the  longitudinal  single-particle  trajectories 
are  shown  in  Fig.  2-3  for  the  compression/ solenoidal  configuration  with  and  without 
the  stellarator  fields.  Also,  Fig.  2-4  shows  the  transverse  beam  emittance  for  the  same 
two  cases.  It  can  be  seen  from  the  plots  that  the  quadruople  compensation  coil  and 
stellarator  coil  placement  provide  minimum  disruption  of  the  beam  as  compared  with  the 
compression/ solenoidal  system  alone. 

2.2  Beam  Stability 

Three  fundamental  types  of  stability  have  been  investigated  for  the  stellatron  con¬ 
figuration:  single  particle  stability,  envelope  stability  and  collective  stability.  The  basic 
stability  criterion  for  single  particle  motion  in  a  linear,  error-free  stellarator  field  was  given 
in  reference  29;  two  pass-bands  and  two  stop-bands  were  found  to  exist  in  the  parameter 
space  describing  the  focusing  system.  The  derivation  in  reference  29  includes  the  dc  space 
charge  effect  of  a  uniform,  round  beam.  In  the  sections  which  follow  we  will  assume  that 
the  fundamental  stability  criteria  are  satisfied  and  go  on  to  discuss  the  types  of  instability 
which  can  still  occur  under  these  conditions. 

2.2.1  Orbital  Resonances 

Simple  considerations  balancing  external  focusing  forces  against  the  defocusing  effects 
of  space  charge  lead  to  very  large  estimates  for  the  limiting  current  in  cyclic  accelerators. 
Much  more  stringent  limits  become  apparent,  however,  when  constraints  on  beam  stability 
are  investigated.  One  interesting  potential  source  of  instability  is  the  resonant  excitation  of 
betatron  oscillations  due  to  field  and  focusing  errors.  When  extra  focusing  fields  axe  added 
to  the  basic  weak  focusing  betatron  field,  the  possibilities  for  exciting  disruptive  orbital 
resonances  may  become  quite  numerous.  Avoidance  of  certain  dangerous  resonances  is 
essential  for  beam  confinement.  Since  the  betatron  frequencies  are  in  general  dependent 
both  on  beam  energy  and  current,  avoidance  of  resonances  places  definite  limits  on  these 
accelerator  parameters.40 

When  field  errors  are  present,  a  particle  in  a  stellatron  is  subject  to  a  periodic  per¬ 
turbation  of  its  orbit  as  it  travels  around  the  device.  If  the  period  of  the  perturbation  is 
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system  of  Table  2-1.  a)  Without,  and  b)  with  the  quadrupole/ stellarator 


coils. 


Table  2-1 

Magnetic  Field  Coil  Data 


MAGNETIC  FIELD  COIL  DATA 

Focusing  Coil:  Single-ring  coils  at  axial  position  2,  radius  r,  and  current  I 

FC1:  2  =  22  cm,  r  =  22  cm,  and  I  =  86,000  amp-turns. 

FC2:  2  =  22  cm,  r  =  26  cm,  and  I  =  -80,000  amp-turns. 

FC3:  2  =  5  cm,  r  =  5  cm,  and  I  -6,200  amp-turns. 

FC4:  2  =  26  cm,  r  =  5  cm,  and  I  =  -7,200  amp-turns. 

AGF  -  Axial  Guide  Field  Coils:  Two  stacks  of  62  single-ring  coils  beginning  at  axial 
position  2  =  33.64  cm  and  ending  at  2  =  138.37  cm.  Inner  coil  set  at  r  =  4.585  cm  and 
outer  set  at  r  =  5.518  cm,  both  sets  carrying  current  of  I  =  1934.4  amps. 

Stellarator  Coils:  The  stellarator  coil  set  consists  of  two  sets  of  /  =  2  helical  quadrupole 
coils,  each  consisting  of  two  wires  spread  at  an  angle  a,  carrying  current  I.  The  quadrupole 
set  begins  at  axial  position  2  and  continues  with  a  pitch  length  of  18  cm  with  a  total  length 
L,  and  beginning  at  a  relative  angle  ar  with  respect  to  a  transverse  axis.  All  coils  have 
radius  r  =  4.154  cm  and  current  I  —  3457  amps. 

ST1:  2  =  41.416  cm,  L  =  85.462  cm,  ar  =  45.38°,  a,  =  4.660°. 

2  =  41.812  cm,  L  =  84.670  cm,  ar  =  53.30°,  a,  =  15°. 

Quadrupole  Lenses:  The  quadrupole  lens  is  defined  by  an  axial  beginning  position  2, 
an  axial  length  L,  a  transverse  magnetic  field  gradient  strength  B1,  and  a  relative  angle  to 
the  stellarator  coil  sets  which  is  9  =  45°. 


QL1:  2  =  40.614  cm,  L  —  2  cm,  B'  =  340  G/m. 
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Figure  2-4  RMS  emittance  plots  of  euu>  vs.  t  for  the  compression/solenoidal  coil  system 
of  Table  2-1.  a)  Without,  and  b)  with  the  quadrupole/stellarator  coils. 


not  equal  or  nearly  equal  to  a  period  of  one  of  the  four  modes  of  betatron  oscillation,  then 
the  resulting  driven  oscillation  remains  of  finite,  usually  tolerable  amplitude;  on  the  other 
hand,  if  there  is  equality  or  near  equality  between  the  periods  of  the  driving  force  and  a 
natural  betatron  mode,  transverse  oscillation  amplitudes  can  be  driven  to  large  values  and 
particles  can  be  lost  to  the  walls  or  beam  emittances  can  grow  intolerably. 

It  is  a  fairly  simple  matter  to  design  a  modified  betatron  or  stellatron  so  that  the 
tunes  axe  away  from  dangerous  integer  and  half-integer  values  for  low  current  beams.  The 
dependence  of  the  tunes  on  current  leads  to  limits  on  the  current,  in  the  usual  way,  but 
the  functional  dependence  of  the  tunes  on  current  is  not  simple.  For  example,  in  the  case 
where  only  vertical  and  toroidal  fields  are  present,  the  fast  cyclotron  mode  tune  decreases 
with  increasing  current,  whereas,  the  slow  drift  mode  tune  may  increase  with  current, 
depending  on  parameters.  For  typical  designs,  the  fast  mode  tune  is  actually  only  weakly 
dependent  on  current  while  the  slow  mode  tune,  while  more  sensitive  to  current,  is  so  slow 
(<C  1)  for  all  anticipated  currents  that  it  is  never  near  a  dangerous  resonant  value,  at  least 
when  the  toroidal  field  is  large. 

The  most  basic  resonance  is  the  dipole  resonance  which  occurs  when  a  particle  executes 
an  integer  number  of  betatron  oscillations  around  the  machine.  If  this  integer  is  large, 
the  resonance  is  usually  not  dangerous  and  is  easily  passed  through  during  acceleration. 
(Experimental  evidence  of  this  fact  has  been  reported  by  the  Irvine  group.)  At  low  energies 
the  ratio  of  toroidal  to  vertical  field  is  usually  large  and  the  resonances  encountered  tend 
to  be  of  high  order.  As  acceleration  proceeds,  however,  the  field  ratio  can  drop  to  0(1)  or 
less  and  the  more  dangerous  low  number  resonances  axe  crossed. 

If  we  include  the  most  general  field  and  focusing  errors  allowed  by  Maxwell’s  equations 
the  equations  for  the  betatron  oscillations  in  the  stellatron  become 


(1  +  Ai)x"  + 


-  +  j/cosm0  -  6  4-  ezo  +  ez 


x  = 


^sinm#  +  tzy 


(2-11) 


y  +  (&  +  eg0)y'  +  Ai  -  e2 


(1  +  A1)y,,+ 


—  —  fxcosmO  +  S  —  ez 
£ 


psinm#  +  e*y  +  e8o 


(2-12) 


x  -  (b  +  eg0)x'  +  er 


where  Ai  is  the  momentum  error,  fi  is  a  dimensionless  quadrupole  strength,  m  is  the 
number  of  field  periods  around  the  device  and  b  is  the  ratio  Bg0/BZ0.  ezo,  ezz,  ezy, 
tg0,  and  ero  axe  five  arbitrary,  but  presumed  small,  periodic  functions  of  6  (period  2n) 
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representing  errors  in  Bz ,  dBz/dr,  dBz/dz,  Be,  and  Br  respectively,  and  where  we  have 
written  n  =  \  +  8.  If  we  assume  that  an  error  term  (i.e.,  the  e’s)  behaves  as  ~  cosk#  for 
some  integer  k  and  if  we  treat  the  e’s  as  infinitesimally  small,  then  a  perturbative  treatment 
gives  a  resonance  when  any  of  the  following  conditions  hold: 


m 

~2 


±  u±  =  ±k 


m  , 

m 

2  ±  l/± 

± 

2  ±  "± 

=  i:k 


where  the  linear  tunes  u±  are  given  by 


v±  =  n  +  ^m2  ±  (nm2  +  ft2)1/2 
m  —  m  +  6/(1  +  Aj) 

n  =  -/(l  +  A,)  +  — /(1+ A,)2 

f*  =  M  1  +  Ai) 


(2-13) 

(2  -  14) 


(2  -  15) 
(2  - 16) 

(2  -  17) 
(2  -  18) 


and  where  all  (±)  signs  in  (2-13,  14)  may  be  chosen  independently.  The  tunes  are  real 
when  the  stability  condition 


-  1-2 
n  —  -m 
4 


(2  -  19) 


is  satisfied.  Clearly  (2-13)  represents  integer  resonances  and  (2-14)  (a  slight  generalization 
of)  half-integer  resonances;  (2-14)  gives  the  location  of  the  first  stop-band  for  infinitesimal 
e’s  while  approximate  locations  of  higher  stop-bands  involve  additional  terms  of  ±(m/2  ± 

i/±). 

When  the  error  terms  are  of  finite  size  (still  small,  but  no  longer  infinitesimal)  one 
employs  Floquet’s  theorem  to  locate  the  stop-  and  pass-bands  in  a  standard  way. 

We  present  results  using  Floquet’s  theorem  below  for  three  special  cases;  in  all  of  these 
we  will  set  m  =  12  and  6  =  10. 

Case  1:  All  e’s  vanish,  but  6^0.  This  case  corresponds  simply  to  introducing  an 
asymmetry  in  the  weak  focusing  field  between  radial  and  vertical  directions  (n  |).  In  Fig. 
2-5  we  show  the  stable  and  unstable  regions  of  the  \i  —  6  plane.  The  resonance  conditions 
t/+  —  j/_  =  m,  of  type  (2-14)  with  k  =  0,  occurs  when  y.  =  59.75.  Unstable  regions  of  the 


plane  are  crossed  with  solid  lines,  contours  of  equal  growth  rates.  (The  finite  resolution 
of  the  grid  may  miss  some  growth  rate  contours.)  The  five  disjoint  regions  of  instability 
have  the  following  origin:  The  two  regions  adjacent  to  the  6-axis  are  the  extensions  into 
regions  of  finite  n  of  the  usual  (fx  =  0)  betatron  stability  condition  0  <  n  <  1.  The  regions 
adjacent  to  the  line  /x  =  59.75  is  the  resonance  region  where  v+  —  u-  «  m,  and  the  band 
for  /i  >  95.5  is  the  basic  instability  region  where  (2-19)  is  violated  for  our  parameters. 

On  Fig.  2-5  two  points  have  been  labelled  A  and  B  in  the  stable  and  unstable  regions, 
respectively.  In  Fig.  2-5A  we  show  results  of  a  particle  orbit  integration  for  a  stellatron 
operated  at  point  A;  the  orbit  shown  is  the  projection  in  the  r-z  plane.  The  orbit  is  clearly 
well  behaved.  In  Fig.  2-5B  we  show  results  of  a  similar  orbit  integration  for  operation 
at  B.  The  orbit  is  plotted  versus  time  in  Fig.  2-5C  where  we  can  see  clearly  the  growth 
saturation,  and  turn-over  of  the  oscillation  amplitude  due  to  nonlinearities.  While  the 
saturation  amplitude  is  not  small,  this  resonance  does  not  seem  to  be  destructive  and  only 
appears  to  couple  energy  back  and  forth  between  two  oscillation  modes,  as  in  the  modified 
betatron.49 

Case  2:  ezo  =  e  cos  3$;  all  other  errors  vanish.  This  case  considers  a  periodic  error 
in  the  vertical  field.  We  illustrate,  in  Fig.  2-6,  the  fx  —  e  plane  where  two  regions  of 
instability  are  seen;  there  is  the  fundamental  instability  band  where  (2-19)  is  violated  for 
/i  >  95.5  and  a  second  band  at  fx  fts  89  at  which  one  finds  the  condition  r/+  —  i/_  =  m  +  k 
is  approximately  satisfied;  here  u+  ss  17,  i/_  «  2,  m  =  12,  k  =  3.  Results  of  particle 
orbit  integrations  at  points  A  and  B  are  given  in  Figs.  2-6A  and  2-6B.  In  this  case  we 
must  choose  a  rather  small  value  of  e  for  our  orbit  integrations;  otherwise  the  resulting 
energy  error,  represented  by  the  inhomogeneous  term  of  (2-11),  causes  an  unacceptably 
large  driven  (but  non-resonant)  oscillation.  (Recall  that  all  the  error  terms  in  (2-11,  12)  are 
normalized  to  the  vertical  field.)  The  only  reason  Fig.  2-6  includes  such  large  e  values  is 
so  that  the  stop-band  could  be  resolved  by  the  mesh  we  used  when  applying  the  numerical 
stability  test.  Near  the  fx- axis,  the  stop-band  is  extremely  narrow. 

It  is  interesting  to  note  that  for  n  «  81  there  is  an  integer  resonance  of  the  form 
(2-13),  m/2  —  ia_  =  k,  (i/-  =  k  =  3);  since  such  a  resonance  is  due  to  the  inhomogeneous 
part  of  (2-11),  it  does  not  appear  els  an  unstable  band  in  Fig.  2-6  which,  as  with  Fig.  2-5, 
is  obtained  by  analysis  of  only  the  homogeneous  part  of  (2-11,  12).  An  orbit  integration 
(not  illustrated)  shows  very  strong  disruption  on  this  resonance,  as  expected  of  a  low  order 
integer  resonance. 

Case  3:  eg0  =  e  cos  60;  all  other  errors  vanish.  For  our  final  case  we  consider  a  small 
periodic  bumpiness  in  the  toroidal  field,  els  would  be  produced  by  discrete  toroidal  field 
coils.  For  this  case,  the  stable  and  unstable  regions  of  the  /x  —  e  plsine  are  shown  in  Fig. 
2-7.  Here  we  see  two  unstable  regions,  the  by-now-familiar  band  for  fx  >  95.5,  and  a  second 
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Figure  2-5  The  /x  —  6  plane  for  case  1.  The  contours  are  lines  of  constant  growth  rate. 

Points  A  and  B  refer  to  cases  for  which  orbits  are  shown  in  Figures  1A  and 
IB  respectively. 


Figure  2-5A  A  particle  orbit  projection  in  the  transverse  p! 
1,  point  A.  ( fx  =  40.,  6  —  4.) 


Figure  2-5B  A  particle  orbit  projection  in  the  transverse  p! 
2,  point  B.  (n  =  59.75,  6  =  4.) 


Figure  2-6  The  /x  —  e  plane  for  case  2. 
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near  the  point  fj.  fa  81,  where  one  finds  that  V-  fa  3  and  a  condition  of  the  form  (2-14) 
is  satisfied  for  k  —  6.  Again  we  have  labelled  points  A  and  B  in  the  stable  and  unstable 
regions;  Figs.  2-7A  and  2-7B,  respectively,  illustrate  results  of  particle  orbit  integrations. 
In  this  case,  the  effects  of  non-linearities  are  insufficient  to  prevent  the  particle  orbit  from 
intersecting  the  pipe  wall,  located  in  this  example  4  cm  from  the  axis. 

Many  other  examples  could  be  given,  of  course;  the  above  three  cases  were  chosen 
only  to  illustrate  some  possibilities. 

To  give  some  feeling  for  the  global  resonance  situation  we  show,  in  Fig.  2-8,  the 
contours  in  the  /x  —  b  plane,  for  fixed  m,  on  which  certain  of  the  conditions  (2-13,  14)  are 
satisfied.  Clearly  it  is  difficult  if  not  impossible  to  avoid  all  resonances.  The  stellatron 
experiment  at  UC  Irvine8  has  demonstrated,  however,  that  many  high  order  ( k  >  20) 
resonances  may  be  passed  through  with  virtually  no  beam  disruption.  Plots  of  the  four 
single  particle  tunes  versus  time  for  a  particular  shot  of  the  Irvine  device  are  given  in  Fig. 
•2-9,  clearly  showing  the  crossing  of  integer  values.  A  plot  of  beam  current  versus  time 
for  the  shot  is  given  in  Fig.  2-10  from  which  one  sees  that  the  current  is  lost  gradually, 
over  the  course  of  50  fj, s  during  the  second  half  of  the  shot,  not  in  the  abrupt  manner 
typical  of  an  integer  resonance  crossing.  (An  electron  makes  fa  115  revolutions  around  the 
device,  per  microsecond.)  It  may  be,  however,  that  while  no  single  high  order  resonance  is 
responsible  for  the  current  loss,  the  net  effect  of  crossing  many  of  them  in  rapid  succession 
may  be  to  induce  a  type  of  slow,  diffusive  particle  loss.  The  above  three  cases,  and  other  we 
have  considered,  support  one’s  intuition  that  low  order  (small  k)  resonances  will  be  more 
of  a  problem  than  large  k  resonances  and  may  well  limit  the  current  and/or  the  energy 
attainable  in  modified-betatron  type  devices. 

When  space  charge  is  included  in  the  calculation  of  the  tunes,  one  finds  that  two  tunes 
are  decreased  and  two  are  increased  by  the  effect  of  space  charge.  The  tunes  all  depend 
on  current  through  the  combination  current/6,  from  which  we  may  draw  the  important 
conclusion  that  the  tunes  are  less  sensitive  to  current  the  stronger  is  the  toroidal  field. 
This  is  an  important  feature  of  the  toroidal  field  and  immediately  suggests  the  use  of  a 
large  toroidal  field  for  control  of  the  tune  shifts  for  large  currents.  While  this  is  correct  in 
principle,  a  practical  consideration  may  limit  its  application.  Specifically,  the  toroidal  field 
needs  to  be  large  at  injection  time  when  space  charge  effects  are  large;  during  subsequent 
acceleration,  then  b  (=  Bg/Bz)  must  be  kept  fixed  in  order  to  keep  the  tunes  fixed,  that  is, 
Be  must  be  raised  in  direct  proportion  to  the  electron  momentum.  If  Be  starts  out  large, 
(initial  fields  of  the  order  of  1  kG  are  typical  in  present  experiments),  practical  limits  to  the 
final  beam  energy  become  clear.  Also,  a  very  large  toroidal  field  will  probably  complicate 
the  beam  ejection  problem  at  the  end  of  the  acceleration. 
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Figure  2-7 A  A  particle  orbit  projection  in  the  transverse 
3,  point  A.  ( fi  =  40.,  t  —  0.75) 


Figure  2-8  The  n~b  plane  for  m  =  12,  showing  con 
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Figure  2-9  Tunes  versus  a  dimensionless  time  for  shot  #36772  of  the  UC  Irvine  stella- 
tron;  the  time  interval  shown  is  «  90  /is. 


(V)  ql 


Figure  2-10  Current  versus  time  for  shot  # 


The  resonances  described  so  far  axe  all  consequences  of  linear  field  errors,  that  is,  the 
single  particle  equations  of  motion  are  still  linear  in  the  transverse  displacements  x  and 
y  in  the  presence  of  the  field  errors.  A  much  expanded  class  of  resonances  is  possible  if 
one  considers  the  effect  of  non-linearities  in  the  applied  fields.  A  fully  general  treatment 
of  non-linear  resonances  would  be  very  complicated  but  a  single  interesting  example  was 
studied50  and  can  be  reported  on  here. 

We  consider  the  magnetic  field  due  to  four  wires  carrying  current  (/,  — /,  —I) 
wound  on  a  cylinder  of  radius  r0  in  a  helix  of  period  2n/k.  To  leading  non-linear  order 
the  fields  in  cylindrical  coordinates  are 


Br(r,  9  —  kz )  ~  B0kr  sin  2(0  —  kz ) 


l-f-(fcr)2 


B#(r,  9  —  kz)  ~  Bakr  cos  2(9  —  kz) 
Bz(r,  9  —  kz)  ~  — B0(kr )2  cos  2(0  —  kz) 


1  +  kr )2 


where  B0  —  16 (Ik/c)  kr0K'2(2kr0).  In  normalized  cartesian  coordinates,  2 kx 
y,  2 kz  — ►  2,  and  2 kct  — >  t,  the  equations  of  motion  are 


(2  -  20) 

(2-21) 

(2  -  22) 

x,  2 ky  —> 


Jl  . 

x  =  a<-y 


+  z 


y  = 


(x2  —  y2)  cos  z  +  2  xy  sin  z 
c  ^1  +  ^x2  +  iy2^  cos  z  +  y  (  1  + 

(x2  —  y2 )  cos  z  +  2 xy  sin  z 


iy2)sinz]} 


(2  -  23) 


+  z  y  l  + 


\x2  +  1 2y 2)  cosz  ~x[l  +  smz  } 

z  =  -a  jxjx  ^1  +  y~x2  +  ^y2^  cos 2  +  y  ^1  +  ^y 


(2  -  24) 


-  y  y  ( 1  +  ^x2  +  ^y2^  COS2  -  X  ^1  +  ^x2  j  sin 2  | 


(2  -  25) 


where  a  =  eBo/Arwykc2 .  These  equations  have  a  constant  of  motion,  Pz  +  kPg  where  PZtg 
are  canonical  momenta: 


yx  —  yx  +  2z  +  a 


(x2  —  y2 )  cos  2  +  2xy  sin  z 


i  +  ^(z2  +  y2) 


=  constant  .  (2  —  26) 
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The  linear  equations  obtained  from  Eqs.  (2-23,  24,  25)  by  neglecting  all  non-linear 
terms  have  solutions  which  oscillate  in  z  with  frequencies  \  ±  k±  where 


k+  = 


1  1/2 


J±(«/A) 


(2  -  27) 


and  where  /3Z  =  zjc.  The  basic  stability  condition  for  these  betatron  oscillations  is 


a 


Pz 


1 

<4 


(2  -  28) 


If  we  return  to  the  non-linear  equations  and  treat  the  non-linearities  perturbatively  we 
see  that  a  non-linear  term  like  x3  cos  z  will  oscillate  at  frequencies  \  +  q  and  |  +  q  where 
q  consists  of  a  sum  of  three  of  ±k±.  In  particular,  the  \  +  q  term  will  be  resonant  if 
q  =  k+  —  2k-  =  k-  or 


k+  =  3k- 


(2  -  29) 


which  occurs  when 


a/0,  =  0.2  .  (2-30) 

Integration  of  the  non-linear  equations  shows  a  dramatic  instability,  even  for  a/ 0Z  close 
to  0.2.  As  an  example,  if  we  take  7  =  3  this  gives  a  resonant  value  for  a  of  0.1886.  Figure 
2-11  shows  the  results  of  a  numerical  integration  of  Eqs.  (2-23,  24,  25)  for  a  =  0.175.  The 
solution  to  the  linear  equations  for  this  a  is  of  course  completely  well  behaved. 

The  presence  of  a  longitudinal  magnetic  field,  in  addition  to  the  twisted  quadrupole, 
changes  the  values  of  k±  and  the  stability  condition  but  not  the  resonance  condition  (2-29). 
k±  becomes 


where 


+  2o.  +  2  ±  [«;(«,  +  2 )2  +  4 


(2-31) 


a:  =  -eBz/(m0yc2k)  (2  -  32) 

ot0  =  4 a/0z  (2  -  33) 


and  the  linear  stability  condition  becomes 
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t  versus  time  for  a  =  0.175. 


a*  +  1|  >  |a„ 


(2  -  34) 


the  resonance  condition  becomes 

[\az  +  1|2  -  9][9|a,  4-  1|2  -  1]  +  100a2  =  0  . 
It  follows  that  no  resonance  is  possible  if 


|a2  +  1|  <  1/3  or  [a*  +  1|  >  3  . 


(2  -  35) 


(2  -  36) 


Figure  2-12  shows  the  |a2  +  1|  —  |a0|  plane  and  the  region  of  linear  stability.  The  dotted 
curve  is  a  plot  of  the  non-linear  resonance  condition,  Eq.  (2-35),  which  must  be  avoided 
in  practical  designs. 

2.2.2  Envelope  Stability  of  l  =  0  Systems 

As  alternatives  to  a  twisted  quadrupole  or  t  =  2  stellarator  field,  poloidally  symmetric 
£  =  0  fields  have  been  investigated  for  use  in  high  current  beam  transport.  Two  basic  types 
have  been  considered,  the  bumpy-torus  betatron51  which  uses  toroidal  field  lines  which 
bulge  in  and  out  periodically  and  the  reversing  solenoid  lens  accelerator52  which  uses  the 
focusing  action  of  periodically  reversing  solenoids  which  produce  field  cusps  at  their  ends. 
We  have  considered33  the  stability  of  K-V  beam  envelopes  in  the  presence  of  space  charge 
in  these  two  focusing  systems,  following  the  approach  of  Struckmeier  and  Reiser54  who 
treated  the  stability  problem  for  i  =  0  and  £  =  2  discrete  element  systems. 

The  K-V  envelope  equations  are 


X"  +  kl(s)X  - 

r"  +  kl(s)Y  - 


2 

-  -  -25-  =0 


2  K  e 


X  +  Y  X3 


2  K 


-  -2-  =  0 

■w  r'i  v 


X  +  Y  Y3 


(2  -  37) 
(2  -  38) 


where  kx  y(s)  represent  external  linear  focusing,  K  is  the  beam  perveance,  and  eIiS,  axe 
unnormalized  beam  emittances.  In  a  torus  with  t  =  0  focusing  these  K-V  equations  hold 
in  the  Larmor  frame  with  ez  =  ey  =  e  and 

**  =  *!  3  *2  =  2^1  (‘ +  ij2(i))  (2  -  39) 
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where  b(s )  =  B$0(s)/ BZo  and  we  have  taken  the  betatron  field  index  n  =  1/2.  We  will  take 
b(s)  to  be  periodic  with  period  S.  In  normalized  variables  the  envelope  equations  become 


(2  -  40) 
(2-41) 


where  the  dimensionless  variable  t  =  ns/ S  has  been  defined  so  that  the  focusing  period  is 
7r.  In  Eqs.  (2-40,  41),  a  dot  means  d/dt  and 

g  =  KSj  ne 

(x,y)  =  (X,Y)/(eSI, r)'/2  . 


The  two  cases  for  b(t)  corresponding  to  the  bumpy  torus  (BT)  and  reversing  solenoid 
lens  (RSL)  are  shown  in  Figure  2-13.  In  each  case  one  wishes  to  know,  (1)  under  what 
circumstances  do  matched  beam  solutions,  i.e.  periodic  solutions  with  period  7r  exists  to 
Eqs.  (2-40,  41),  and  (2)  when  are  these  solutions  stable  to  small  perturbations?  The 
second  of  these  questions  is  of  particular  practical  importance  because  one  would  like  to 
know  whether  a  beam  launched  with  nearly  the  correct,  matched  initial  conditions  will 
remain  nearly  matched  as  it  propagates  through  the  focusing  system. 


To  find  matched  beam  solutions  we  set  x  =  y  =  R  in  Eqs.  (2-40,  41)  since  a  cir¬ 
cular  matched  beam  is  expected  in  a  symmetric  focusing  system  for  a  beam  with  equal 
emittances  in  the  two  transverse  planes.  One  then  may  find  initial  conditions  -ft(O)  and 
f?(0)  numerically  such  that  the  resulting  solution  is  periodic  with  period  7r;  such  a  solution 
exists  when  the  phase  shift  per  period  of  single  particle  motion 


H  = 


dt 

E?<lr  ’ 


(2  -  42) 


which  is  the  condition  for  stability  of  betatron  oscillations.  In  the  absence  of  space  charge 
(g  =  0),  fj.  =  (i0;  in  general  p  <  fi0. 

To  investigate  stability  of  the  matched  solutions  one  carries  out  a  straightforward 
linearization  of  Eqs.  (2-40,  41),  writing 


—  R  +  <5r 


y  —  R  -f-  by 


(2  -  43) 
(2  -  44) 
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External  Focusing  Function 


s 


2nr 
_ o 


b2(t+77)  -  b2 ( t ) 


b2(t)  -  b2  (1+  €  cos  2t ) 2  (BT) 


b2(t) 


b2  0  *  t  t-  tQ/2  and  rr-to/2  l  t  *  n 
0  tQ/2  <  t  t r  -  tQ/2  (RSL) 


Figure  2-13  External  focusing  functions  for  bumpy  torus  (BT)  and  reversing  solenoid  lens 
(RSL)  systems. 
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Defining  6±  =  8Z  ±6y  one  finds  that  the  equations  for  6+  and  8-  are  decoupled,  representing 
independent  modes  of  oscillation: 


h  + 


g  3 
R 2  +  R 4 


6+  =  0 


L  + 


(2  -  45) 
(2  -  46) 


The  “+”  mode  is  a  monopole  or  “breather”  mode  and  the  mode  is  a  quadrupcle  or 
elliptical  mode. 

Equations  (2-45,  46)  are  linear  equations  with  periodic  coefficients  to  which  Floquet’s 
theorem  applies.  The  stability  condition  is 


He  <  *  (2  ~  47) 

where  fie  is  the  phase  shift  per  period  of  either  the  <5+  or  8-  equations.  This  is  a  more 
stringent  condition  than  Eq.  (2-42),  the  condition  for  single  particle  stability.  In  fact  we 
have  found  that  Eq.  (2-47)  implies  that 


Ho  <  */2  (2  -  48) 

for  the  bumpy  torus  and  RSL  systems;  this  is  also  the  envelope  stability  condition  found 
by  Struckmeier  and  Reiser.54  When  (2-48)  is  satisfied  the  beam  current  limit  is  set  only 
by  the  aperture  size  or  by  some  other  collective  effect. 

2.2.3  Negative  Mass  Instability 

The  negative  mass  instability  in  the  modified  betatron  was  studied  by  Sprangle  and 
Vomvoridis.26  They  found  that  the  longitudinal  field  coupled  transverse  and  longitudinal 
modes  and  that  the  field  helped  to  reduce  the  growth  rate  of  the  instability.  Godfrey  and 
Hughes54  subsequently  improved  the  dispersion  relation,  and  resolved  a  puzzle  about  the 
existence  of  a  growth  rate  cutoff  above  a  certain  critical  current,  by  including  toroidal 
corrections  to  the  first  order  fields. 

When  strong  focusing  fields  are  added  an  important  feature  of  the  dispersion  relation 
is  changed.40  Since  the  momentum  compaction  factor  a  is  greatly  reduced,  to  a  value  less 
than  one,  by  the  presence  of  the  strong  focusing  fields,  a  transition  energy  is  introduced 
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below  which  there  is  no  negative  mass  effect;  the  beam  consequently  is  stable  for  ener¬ 
gies  below  mc27t.  This  transition  energy  is  derived  in  terms  of  system  parameters  (field 
strengths,  periods,  etc.)  in  reference  40,  attached  here  as  an  appendix. 

Reference  40  derives  the  full  linear  eigenvalue  problem  for  long  wavelength,  low  fre¬ 
quency  collective  modes  in  the  stellatron.  The  eigenvalue  problem  is  then  solved  in  the 
‘smooth  approximation,’  which  assumes  that  the  wavelength  of  the  perturbation  is  many 
stellarator  field  periods  long.  This  results  in  a  dispersion  relation  identical  in  form  to  that 
for  the  modified  betatron,  under  the  replacement  of  a  single  term;  the  dispersion  relation 
has  been  tested  against  a  numerical  solution  of  the  eigenvalue  problem  and  excellent  agree¬ 
ment  has  been  found.  The  net  effect  of  the  stellarator  windings,  we  have  found,  is  to  add 
a  smooth,  transverse  focusing  term  which  can  be  large  compared  to  the  weak-focusing  or 
space  charge  defocusing  terms;  stellarator  fields  enhance  stability  in  this  sense;  as  a  result, 
increasing  the  stellarator  field  strength  reduces  the  negative  mass  growth  rate.  Numerical 
examples  given  in  reference  40  quantify  this  effect. 

2.2.4  Beam-Breakup  Instability  (BBU) 

The  BBU  is  the  most  serious  instability  facing  the  SLIA.  The  magnitude  of  the  prob¬ 
lem  can  be  understood  by  examining  the  asymptotic  formula  for  the  maximum  beam 
displacement  during  the  beam  pulse  in  the  n-th  gap,  as  given  by  V.K.  Neil  et  a/.,55 


;  exp(fn) 

~  /  —  „  -  > 
Vl67!Tn 

where 


(2  -  50) 


f„=  5.49x10-^10]-^  ,  (2-51) 

and  where  u>  =  deflecting  mode  frequency,  a>o  =  27t  (785  MHz),  I  =  beam  current,  and  Bo 
=  solenoidal  magnetic  field. 

This  maximum  amplitude  occurs  at  a  time  tn,  given  by 


uin  =  2TnQ  ,  (2-52) 

where  Q  is  the  quality  factor  for  the  accelerating  gap.  The  beam  pulse  is  assumed  to  be 
longer  than  tn. 

The  coupling  impedance,  Zj.,  is  given  in  the  limit  b^>  w  by56 


Z±\n]  =  ~730^~  ImPi(u;) 


(2  -  53) 
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where  w  is  the  gap  width,  b  is  the  pipe  radius,  and  Pi(w)  is  the  cavity  response  function. 
Substituting  this  formula  into  the  expression  for  r„  yields 


f,  imV/lGeV]fM/mf,(u,) 

E,[MVIcm\Ba[kG)V  ' 


(2  -  54) 


where  Eg  is  the  electric  field  in  the  gap,  Egw  is  the  energy  gain  per  gap,  and  Vj  is  the 
final  beam  energy.  In  the  limit  b  «  w,  Godfrey57  has  shown  that  the  impedance  falls  as 
b~ 1  instead  of  b~2.  The  ATA  gaps  satisfy  the  expression  given  above,  but  the  SLIA  gaps 
will  be  better  represented  by  Godfrey’s  formula. 

Figure  2-14  shows  contours  of  constant  rn  =  10,  corresponding  to  £„  «  1000,  plotted 
in  the  Eg  —  b  plane  for  various  combinations  of  I  and  Vf,  assuming  ImPi(u )  =  -0.7.  There 
is  a  clear  advantage  to  designing  gaps  having  as  large  as  an  Eg  as  possible.  At  present 
LLNL  operates  at  «  200  kV/cm,  while  Sandia  has  operated  their  gaps  up  to  fa  400  kV/cm. 
The  limit,  set  by  emission  levels  over  areas  of  interest,  must  be  experimentally  studied  with 
the  use  of  new  electrode  materials,  coatings  and  conditioning  techniques. 

The  cavity  response  function  has  been  computed  for  the  deflecting  mode  in  ATA  as 
ImPi(uj )  =  -1.3.  For  a  perfectly-coupled  gap  (ie.  with  no  reflections)  the  response  would 
be  roughly  half  the  ATA  value,  ImPi(u> )  fa  -0.7.  LLNL  has  designed  gaps  with  nearly 
optimal  coupling,  although  they  have  not  yet  measured  their  BBU  characteristics. 

The  BBU  can  be  detuned  through  the  use  of  nonlinear  focusing  fields.  Solenoidal 
and  quadrupole  focusing  forces  are  both  linear  in  the  transverse  displacement  from  the 
matched  reference  orbit.  The  use  of  quadrupole  focusing  in  the  SLIA,  therefore,  does  not 
substantially  affect  the  BBU.  It  acts  only  as  an  additional  solenoidal  focusing  field.  In  the 
smooth  approximation  this  additional  focusing  can  be  calculated  to  be58 


Beff  =  B0 


(2  -  55) 


where  tq  is  proportional  to  the  quadrupole  field  strength.  For  e?  <  1,  the  quadrupole 
focusing  yields  only  a  modest  reduction  in  the  BBU  growth  rate  compared  with  the  result 
with  Bq  alone.  By  adding  a  sextupole  field  to  the  transport  system,  however,  the  focusing 
fields  will  be  quadratic  in  the  displacement  from  the  reference  orbit,  and  the  growth  of 
the  BBU  will  be  strongly  affected.46  The  nonlinear  transverse  focusing  causes  the  betatron 
wavelengths  of  the  electrons  to  depend  on  the  amplitude  of  their  displacement  from  the 
reference  orbit.  As  the  BBU  instability  drives  the  beam  electrons  off  axis,  therefore,  their 
response  to  the  BBU  driving  mode  becomes  detuned,  thereby  preventing  them  from  ex¬ 
citing  that  mode  in  accelerating  gaps  downstream.  Ion  focusing  also  defeats  the  BBU  by 
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Figure  2-14  Contours  for  Tn 


introducing  a  spread  in  betatron  wavelength.45  Since  the  SLIA  is  planned  to  operate  with¬ 
out  IFR  guiding,  magnetic  focusing  systems  that  can  spoil  the  BBU  have  been  explored. 
A  twisted  sextupole  winding  appears  to  offer  great  promise  in  this  regard. 

The  twisted  sextupole  field  is  the  equivalent  of  an  /  =  3  stellarator  field,  whose  com¬ 
ponents  may  be  expressed  (to  quadratic  order  in  x  and  y)  as 

27 

Bx  fa  —  t3B0Oi2[2xy  cos(3os)  —  ( x 2  —  y2)  sin(3o;s)]  (2  —  56) 

16 

27 

By  fa  —eaBoa2[(x2  —  y2)  cos(3as)  +  2xy  sin(3as)]  .  (2  —  57) 

These  fields  add  to  the  other  fields  (longitudinal  field,  quadrupole  field,  and  vertical  field 
on  the  bends)  in  the  SLIA  configuration. 

The  current  that  must  be  used  to  excite  a  continuous  sextupole  winding  to  produce 
a  desired  shear  in  the  field  on  axis  depends  strongly  on  the  radius  of  the  sextupole  coil, 
as  shown  in  Figure  2-15.  The  parameter  e3  on  that  figure  measures  the  sextupole  field 
strength.  When  es  =  1,  the  amplitude  of  the  sextupole  field  equals  the  solenoidal  field, 
Bo-  Figure  2-15  is  plotted  for  B0  =  5  kG.  For  the  SLIA,  where  the  vacuum  vessel  has 
a  radius  of  3  cm,  the  current  required  in  the  sextupole  coils  will  be  fa  240  kA  for  e3  = 
6.  While  this  current  is  large,  the  inductance  of  the  sextupole  winding  is  small,  as  is  the 
energy  associated  with  this  field. 

The  effect  of  this  field  is  illustrated  in  Figure  2-16.  This  figure  shows  the  beam  profile 
after  passing  through  800  gaps.  The  beam  pulse  is  30  ns  long,  and  t  =  0  on  the  figure 
corresponds  to  the  time  that  the  beam  head  passes  the  800th  gap.  The  gaps  in  this  example 
are  modeled  by  an  impulse  model  developed  by  V.K.  Neil  et  al.55  The  gap  characteristics 
are  assumed  similar  to  those  in  ATA  (Z±/Q  =  10  fi,  Q  =  4,  a  deflecting  mode  frequency  of 
785  MHz,  and  energy  gain  per  cavity  of  0.3  MeV).  The  final  energy  of  this  accelerator  would 
be  241  MeV  (1  MeV  injected  +  800  gaps  x  0.3  MeV/gap).  The  transport  system  between 
gaps  is  similar  to  that  planned  for  the  SLIA  ( B0  =  5  kG,  quadrupole  field  amplitude  of  eq 
=  0.3,  quadrupole  pitch  length  of  Lq  =  18  cm),  with  the  addition  of  a  continuously-twisted 
sextupole  field  (ea  =  6  and  Ls  =  18  cm).  The  incident  electron  beam  at  the  first  gap  has 
a  current  of  10  kA,  an  energy  of  1  MeV,  and  is  off-axis  by  10-4  cm.  The  beam  is  injected 
without  energy  spread,  and  the  model  allows  no  energy  spread  to  develop  as  the  electrons 
are  accelerated;  energy  spread  should  be  stabilizing  for  the  BBU  since  it  independently 
introduces  a  spread  in  the  betatron  wavelength.  Without  the  sextupole  field,  Figure  2-16 
shows  that  the  amplitude  of  the  beam  displacement  from  the  axis  would  be  fa  108  cm  after 
800  gaps;  it.  the  beam  would  be  catastrophically  unstable.  With  the  sextupole  field,  the 
beam  displacement  is  limited  to  fa  1  cm,  as  seen  on  the  figure.  The  beam  displacement 


22 


Sextupole  Coil  Current 


CD 

uo 

in 

in 

in 

O 

+ 

+ 

+ 

+ 

+ 

+ 

0 

CD 

0 

0 

0 

0 

T— 

00 

CO 

CM 

o 

[V]  juajjno 


Figure  2-15  Sextupole  coil  current  versus  winding  radius. 
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Figure  2-16  Beam  profile  at  gap  800.  (10  kA  beam) 


saturates  at  this  amplitude  before  the  400th  gap,  and  remains  stable  through  the  800th 
gap,  which  was  the  final  gap  in  the  calculation. 

Figure  2-17  shows  the  same  calculation  carried  out  for  a  40  kA  beam.  Again,  the 
sextupole  field  clamps  the  BBU  growth  when  the  displacement  reaches  1  cm  amplitude. 

Other  tests  have  demonstrated  that  the  saturation  of  the  BBU  near  amplitudes  of  1 
cm  is  insensitive  to  the  initial  amplitude.  When  the  initial  amplitude  in  Figure  2-16  is 
increased  from  10-4  cm  to  0.1  cm,  the  saturation  occurs  after  fewer  gaps,  but  the  saturated 
amplitude  is  still  1  cm. 

These  calculations  offer  the  promise  of  defeating  the  BBU  with  high-order  multipole 
fields,  such  as  sextupoles.  The  effect  of  these  fields  on  the  beam  emittance  or  on  beam 
stability  other  than  the  BBU  has  not  yet  been  investigated.  These  studies  are  on-going  at 
the  present  time. 

Shielded  gaps  may  offer  yet  another  means  to  control  the  BBU.  In  the  SLIA  conceptual 
design  shielded  gaps  are  used  to  avoid  the  low-frequency  beam  displacement  in  off-axis 
pipes  due  to  return-current  asymmetries  in  the  feed.  No  measurable  deflection  of  beams 
was  observed  with  shielded  gaps,  within  diagnostic  resolution,  in  the  experiments  of  Wilson 
at  NBS  and  Hasti  at  SNLA.  Miller  et  al.  of  SNLA  have  also  studied  a  version  of  these 
gaps  in  the  RADLAC  development  program.  The  high-frequency  character  of  these  gaps 
is  unknown,  and  must  be  evaluated  in  “cold-test”  cavity  measurements. 

Miller  et  al ,59  have  used  the  SOS  code,  a  3-D  code  developed  at  Mission  Research 
Corporation  (MRC),  to  analyze  the  differences  between  pillbox  cavities  of  the  type  used 
at  ATA  and  the  so-called  “shielded-gap”  configurations  which  are  envisioned  for  the  SLIA. 
Their  model  utilized  perfectly- conducting  structures  everywhere  (including  the  feed  ter¬ 
mination).  and  therefore  does  not  provide  design  data.  The  model  is  useful,  however,  in 
displaying  the  relative  differences  between  pillbox  and  shielded-gap  configurations.  Their 
conclusion  is  that  a  shielded-gap,  made  by  inserting  a  short  length  of  coaxial  transmission 
line  in  a  pillbox  cavity,  can  substantially  reduce  the  coupling  impedance  to  the  BBU  in¬ 
stability,  provided  that  the  coax  is  designed  so  that  the  lowest  TE  mode  cutoff  frequency 
is  greater  than  the  frequency  of  the  TM  cavity  mode  which  drives  the  BBU.  More  recent 
calculations  at  MRC60  and  at  LLNL61  have  found  that  shielded-gap  configurations  can  in¬ 
troduce  new,  weakly-damped  modes  which  increase  the  BBU  coupling  impedance.  These 
calculations  have  so  far  not  found  any  shielded-gap  design  that  substantially  improves  the 
BBU  coupling  impedance  over  that  obtainable  with  a  simple  radial  line.  This  subject  is 
still  under  study  at  SAIC,  MRC,  and  LLNL. 
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Figure  2-17  Beam  profile  at  gap  800.  (40  kA  beam) 


2.3  The  Spiral  Line  Induction  Accelerator  (SLIA) 

Our  investigations  of  the  properties  of  high  current,  strong  focused  accelerators  led  to 
the  joint  proposal,  in  1985,  of  what  is  now  called  the  SLIA,  by  SAIC  and  Pulse  Sciences, 
Inc.  The  SLIA,  illustrated  in  Figure  2-18,  consists  of  a  continuous  beam  line  which  is 
recirculated  through  an  induction  module  multiple  times;  the  number  of  recirculations 
depei  T  on  the  sizes  of  the  beam  pipe  and  the  module  cores. 

The  SLIA  as  proposed  has  several  advantages  over  closed  orbit  machines.  In  particular, 
injection  and  ejection  become  much  easier  than  for  a  closed  orbit  device.  In  addition,  the 
open-ended  design  of  the  SLIA  allows  very  long  pulse  lengths  to  be  envisioned;  the  pulse 
length  is  set  by  the  injector  pulse  length  rather  than  by  the  circulation  time,  as  in  a  closed 
orbit  machine.  The  magnetic  fields  used  for  focusing  and  bending  in  the  SLIA  beam  line 
are  effectively  static  in  time,  simplifying  their  design  and  allowing  low  field  errors  to  be 
achieved.  It  may  even  be  possible  to  use  permanent  magnets  for  some  of  the  fields  in 
the  transport  line.  These  advantages  come  at  the  cost  of  added  mechanical  and  electrical 
complexity,  however.  The  beam  transport  line  design  is  complicated  by  the  simultaneous 
presence  of  many  coils  and  the  induction  module  design  is  quite  novel,  requiring  careful 
considerations  to  reduce  asymmetries  in  the  accelerating  fields  applied  to  the  different  gaps. 
Nonetheless,  based  on  a  1985  proposal  by  SAIC  and  PSI,  a  SLIA  experiment  was  funded 
and  a  transport  line  was  built  to  study  the  properties  of  a  strong  focused  recirculator.  This 
experiment  is  now  planned  to  lead  to  a  2  module,  3-turn  machine  producing  an  8.5  MeV, 
10  kA  beam. 

The  research  funded  by  ONR  under  this  contract  contributed  and  continues  to  con¬ 
tribute  directly  to  the  evolution  of  the  SLIA  by  its  concentration  on  fundamental  issues  of 
beam  transport  and  stability  in  strong  focused  systems.  The  original  work  by  Roberson, 
Mondelli  and  Chernin  on  momentum  compaction  in  the  stellatron  has  been  basic  to  the 
design  of  the  recirculator.  The  work  performed  on  beam  stability  and  envelope  matching 
has  also  facilitated  the  design  of  the  beam  line.  One  of  the  most  fundamental  physics  is¬ 
sues  for  the  SLIA,  the  beam  break-up  instability,  has  been  extensively  studied  under  ONR 
sponsorship,  resulting  in  some  promising  proposals  for  stabilization  by  non-linear  fields 
and  by  energy  spread.  Work  on  orbital  resonances,  while  most  pertinent  to  closed  orbit 
devices  in  which  a  beam  circulates  many  times,  may  also  find  application  in  the  SLIA,  or 
certainly  in  the  induction  synchrotron  accelerator  (ISA)  for  which  it  has  been  proposed  to 
use  the  SLIA  as  an  injector. 

Though  a  number  of  important  physics  and  design  issues  remain  to  be  resolved,  the 
technology  represented  by  the  SLIA  appears  to  be  a  promising  direction  for  future  research 
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FIXED  VERTICAL  FIELDS 


Figure  2-18  Schematic  of  the  Spiral  Line  Induction  Accelerator  (SLIA):  a  two  pass  ex- 


in  high  power  electron  accelerators.  ONR  sponsorship  of  the  basic  physics  and  scaling  of 
strong  focused,  compact  accelerators  has  laid  important  groundwork  for  on-going  research 
in  an  area  of  continuing  interest  to  the  Navy. 
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3.0  Conclusions  and  Recommendations 

Our  studies  of  compact,  strong  focused,  high  current  electron  accelerators  have  ad¬ 
dressed  many  aspects  of  beam  behavior  in  these  devices.  We  have  found  that  systems 
employing  l  =  2  stellarator  windings  can  be  designed  to  be  stable  and  robust,  with  excel¬ 
lent  energy  bandwidth.  Experimental  programs  at  NRL,  UC  Irvine,  and  Pulse  Sciences, 
Inc.  continue  to  investigate  the  properties  of  these  focusing  systems. 

The  criteria  for  stability,  both  due  to  single  particle  resonances  and  to  collective  effects, 
are  reasonably  well  understood  at  this  point,  except,  perhaps,  for  high  frequency  phenom¬ 
ena  like  the  short  wavelength  negative  mass  instability  and  the  electromagnetic  three-  wave 
instability.  It  appears  that  current  limits  in  closed  devices  like  the  betatron/stellatron,  axe 
due  to  tune  shifts  and  to  the  negative  mass  instability,  while  in  induction  linacs,  including 
SLIA  configurations,  current  is  limited  by  the  beam  break-up  instability.  Energy  limits  are 
similarly  constrained  by  tune  shifts  in  closed  devices  if  one  realizes  the  condition  required 
to  avoid  resonanaces  is  to  raise  all  fields  proportionally,  in  particular,  BgaBz;  in  the  SLIA, 
energy  limits  will  be  dictated  by  induction  module  technology  and  by  the  number  of  re¬ 
circulations  allowed  by  pipe  and  core  radii.  Finally,  beam  pulse  length  limits  in  a  closed 
device  are  clearly  bounded  by  the  circumference/c  while  in  a  SLIA  or  induction  linac, 
output  pulse  length  is  limited  only  by  injector  pulse  length  and  the  number  of  volt-seconds 
that  can  be  supplied  to  the  cores  (for  a  continuous  beam).  These  conclusions,  and  the 
analyses  that  led  to  them,  represent  the  main  results  of  our  work  on  compact  electron 
accelerators. 

There  remain  several  important  issues  in  these  devices  which  require  further  study. 
Prominent  among  these  is  the  stabilization  of  the  beam  break-up  instability  in  the  SLIA. 
Many  lines  of  investigation  should  be  followed,  especially  the  potential  role  of  energy 
spread,  including  the  variation  in  recirculation  time  with  energy,  and  the  application  of 
non-linear  fields.  Also,  the  use  of  feed-forward  stabilization,  whereby  the  kick  received  by 
the  beam  is  sensed  as  it  exits  a  module  and  a  signal  is  sent  across  a  chord  of  the  device  to 
a  point  at  which  a  correcting  kick  is  applied,  may  have  practical  application  in  the  SLIA. 
Virtually  no  work  has  been  done  on  this  question. 

The  related  issues  of  beam  matching  and  beam  quality  in  these  devices  also  deserve 
more  thorough  treatments.  Techniques  for  matching  a  beam  on  and  off  the  stellarator 
field  in  order  to  minimize  the  resulting  envelope  oscillations,  are  presently  imperfect  and 
much  work  needs  to  be  done  along  these  lines.  Very  little  also  is  presently  known  about 
required  field  tolerances  and  the  effect  of  field  errors  on  beam  behavior  in  general  and  on 
beam  quality  in  particular. 


Though  significant  progress  has  been  made  in  the  understanding  of  compact  electron 
accelerators  further  work,  both  theoretical  and  experimental,  is  essential  to  ensure  that 
practical,  workable  designs  can  be  developed. 
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The  magnetic  field  configurations  consisting  of  the  combination  of  a  weak-focusing  betatron  field, 
toroidal  field  and  either  1  =  0  or  1  =  2  stellarator  windings  are  assessed  for  their  potential  as  focusing 
fields  for  high-current  cyclic  electron  accelerators.  These  accelerators,  named  “stellatrons,”  are  shown 
to  have  unproved  tolerance  to  mismatch  between  the  average  beam  energy  and  the  equilibrium  beam 
energy  that  matches  the  vertical  magnetic  field  compared  to  devices  without  the  stellarator  fields.  Both 
analytical  calculations  in  the  paraxial  approximation  and  numerical  particle-orbit  calculations  are 
presented  to  substantiate  this  finding.  The  problems  of  orbital  resonances  and  of  injection  into  these 
devices  are  discussed.  Earlier  work  in  the  field,  much  of  it  unpublished,  is  discussed  and  compared 
with  the  stellatron  concept. 


I.  INTRODUCTION 

Electron  accelerators  built  to  carry  out  nuclear  and  high-energy  physics  experi¬ 
ments  typically  carry  average  currents  of  less  than  an  ampere  at  kinetic  energies 
as  high  as  50  GeV.  During  the  past  twenty  years,  significant  progress  has  been 
made  in  high-current,  low-energy  beam  generation.  These  intense-beam  ac¬ 
celerators  are  generally  motivated  by  applications,  mostly  for  radiation  source 
development  and  for  controlled  nuclear  fusion  research.  Megampere  currents  at 
megavolt  energies  are  now  routinely  obtained  using  pulselines  charged  by  Marx 
capacitor  banks,  connected  to  cold  cathode  (field-emission)  diodes.1 

Multikiloampere  electron  currents  at  tens  of  megavolts  have  been  obtained 
using  linear  induction  accelerators.2-6  Recently,  there  has  been  considerable 
interest  in  extending  the  current-carrying  capabilities  of  cyclic  induction  ac¬ 
celerators  such  as  variations  of  the  betatron  and  linear  induction  modules  in  cyclic 
(racetrack  or  other)  configurations. 
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Conventional  betatrons7  8  are  current  limited  due  to  the  defocusing  effects  of 
space  charge  at  injection.  Overcoming  the  space-charge  limit  requires  high-energy 
injection.  By  injecting  a  beam  at  high  energy,  approximately  200  A  of  circulating 
current  has  been  obtained  in  a  small  betatron  in  which  the  beam  was  ultimately 
accelerated  to  100  MeV.9  In  another  experiment  using  a  4-MeV  induction  linac, 
approximately  500  A  of  circulating  current  has  been  confined  in  a  conventional 
betatron  configuration.10  Devices  designed  to  improve  the  current-carrying  capa¬ 
bility  of  betatrons  have  included  fixed-field  alternating-gradient  (FFAG)  beta¬ 
trons11  and  plasma  betatrons.12  The  FFAG  uses  alternating  gradient  strong 
focusing  fields1314  in  addition  to  a  vertical  betatron  field.  In  the  FFAG,  the 
magnetic  fields  remain  constant  in  time,  while  the  flux  linking  the  particle  orbit  is 
changed  to  produce  an  inductive  electric  field.  To  avoid  single-particle  resonances 
as  the  energy  is  increased,  the  equilibrium  radius  of  the  orbit  is  allowed  to  vary 
with  energy,  keeping  the  orbits  “self-similar”.  The  currents  achieved  in  pro¬ 
totypes11  have  been  modest,  limited  by  the  tune  shift  due  to  space-charge  effects. 

The  plasma  betatron  employs  a  toroidal  magnetic  field  in  addition  to  the 
betatron  field.  A  plasma  is  injected  or  created  in  the  device  and  the  applied 
inductive  electric  field  causes  a  portion  of  the  electron  distribution  to  gain  more 
energy  from  the  electric  field  between  collisions  than  it  loses  during  a  collision,  a 
phenomenon  called  electron  runaway.54  In  magnetic-fusion  devices,  runaway 
currents  of  hundreds  of  amperes  at  energies  of  several  MeV  have  been 
obtained.1516,39 

The  HIP  AC 17  device,  using  only  a  toroidal  magnetic  field  with  a  cross-field 
injection  scheme  known  as  inductive  charging,18  has  been  investigated  in  connec¬ 
tion  with  an  ion-acceleration  concept.  Average  electron  densities  of  4x  lO’cm-3, 
have  been  achieved.  The  trapped  electrons  are  not  accelerated  in  this  device,  but 
are  used  instead  to  create  a  strong  potential  well  for  ions. 

In  another  series  of  experiments,19  average  electron  densities  of  10locm-3  have 
been  achieved  using  inductive  charging  in  toroidal  device.  A  total  charge  of  100 
microcoulombs  has  been  trapped,  which  if  accelerated  to  relativistic  velocities 
would  result  in  a  current  of  10  kA.  With  a  time-independent  vertical  field  and  a 
transformer  used  to  produce  an  inductive  electric  field,  only  small  electron  current 
(<50  A)  is  obtained.  The  results  indicate  that  the  electric  cloud  is  not  accelerated, 
but  remains  trapped  in  the  torus. 

Current  interest  has  been  focused  on  high-current  nonneutral  electron- 
beam  acceleration.20  Recently,  a  modified  betatron  configuration  has  been  sug¬ 
gested21"23  that  employs  a  conventional  weak-focusing  betatron  field  and  toroidal 
magnetic  field.  The  principal  advantage  offered  by  the  modified  betatron  is  that 
the  toroidal  magnetic  field  reduces  the  required  injection  energy  by  containing  the 
space-charge  defocusing  forces.  Extensive  analysis  of  this  configuration  has  been 
carried  out  both  analytically  and  numerically.24-34 

Recent  experiments  employing  inductive  charging  in  a  modified  betatron  con¬ 
figuration  have  achieved  beam  currents  of  approximately  200  A  and  energies  of 
approximately  !  Me'.’  vsing  ?Q-kV  injection  vcltage.35  In  an  elongated  or 
“stretched"  modified  betatron,  beam  currents  of  approximately  50  A  at  energies 
of  1  MeV  have  been  achieved  using  50-kV  injection.34 
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Another  approach  to  high-current  cyclic  accelerators  is  the  use  of  linear 
induction  modules  in  a  cyclic  configuration.  One  such  scheme3  has  used  a 
long-pulse  linear  induction  module  with  beam  recirculation  through  isolated  beam 
paths  for  eagh  transit  through  the  induction  module.36  This  configuration  is 
essentially  a  folded  induction  linac.  Yet  another  approach  is  to  use  a  long-pulse 
induction  module  in  a  racetrack  geometry.  High  currents  can  be  handled  by 
adding  a  toroidal  magnetic  field.  This  field  introduces  particle  drifts  in  the  bends, 
however,  which  can  be  averaged  out  by  the  addition  of  stellarator  windings.37  For 
reasonable  parameters,  the  stellarator  windings  can  contain  particles  with  energies 
up  to  approximately  l  MeV  per  kilogauss  of  stellerator  field  on  a  1-m  radius  of 
curvature  bend,38  based  on  single-particle  numerical  orbit  integrations. 

If  a  time-dependent  vertical  magnetic  field  is  added  to  the  bends  in  order  to 
guide  the  beam,  then  the  stellerator  windings  provide  bandwidth  for  an  energy 
mismatch  between  the  beam  energy  and  the  matched  energy  in  the  vertical  field. 
The  matched  energy  in  the  vertical  magnetic  field  is  approximately  30  MeV/ 
kG-m.38 

The  bandwidth  can  be  important,  not  only  because  it  allows  one  to  handle 
beams  with  a  significant  variation  in  kinetic  energy,  but  also  because  it  reduces  the 
sensitivity  of  the  system  to  abrupt  changes  such  as  occur  at  accelerating  gaps. 

This  paper  describes  an  analysis  of  a  configuration  consisting  of  a  combination 
of  stellarator  and  betatron  fields,  called  the  stellatron.40'*3  The  motivation,  in 
part,  has  been  to  increase  the  allowed  mismatch  between  the  beam  energy  and  the 
betatron  field.  In  both  the  conventional  and  modified  betatron,  the  allowed 
mismatch  Ay/y0  is  given  by 


— s|(l-n)-. 

To  ro 


where  n  is  the  betatron  field  index  and  o/r0  is  the  inverse  aspect  ratio  of  the 
accelerator  chamber.  Typically  a/r0  =  0.1  and  so  Ay/y0s2.5%  is  the  allowed 
bandwidth  for  mismatch  in  the  conventional  and  modified  betatron  for  n  =  1/2. 
Since  the  vertical  field  at  injection  in  several  planned  experiments  is  generally  less 
than  100  G,  field  errors  of  a  few  gauss  can  cause  loss  of  the  beam. 

By  adding  a  stellarator  field  to  a  cyclic  accelerator,  a  strong-focusing  system  is 
obtained  that  can  sustain  high  currents  and  large  mismatch  between  particle 
energy  and  vertical  field.  Stellarator  fields  are  most  simply  characterized  by  two 
integers,  1  and  m,  these  being  respectively  the  number  of  field  periods  in  the 
poloidal  and  toroidal  directions  in  the  device  [see  Eq.  (3)  and  Fig.  3 A  below). 
Two  special  cases  have  been  treated  in  some  detail.  The  1  =  2  stellatron4041  is 
shown  in  Fig.  1.  The  stellarator  field  consists  of  a  toroidal  field  plus  a  continuously 
twisted  quadrupole  field.  The  twisted  quadrupole  configuration  is  analogous  to  the 
alternating-gradient  strong-focusing  fields  that  are  routinely  utilized  in  modern 
synchrotrons. 

The  l  =  0  stellatron  shown  in  Fig.  2  is  similar  to  a  bumpy-torus  fusion  device 
with  the  addition  of  the  betatron  field.  Here,  the  varying  radial-field  component 
leads  to  strong  focusing.  The  I  =  1  configuration  has  no  transverse  field  gradient  at 
the  beam  axis,  and  therefore  is  only  weakly  focused.  Stellatrons  with  is 3  are 
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£-2  SYSTEM  -  "STEUATRON 
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FIGURE  1  The  I  =  2  steUatron. 

high-shear  devices  with  low  fields  near  the  magnetic  axis.  The  nonlinear  proper¬ 
ties  of  these  configurations  may  be  useful  in  conjunction  with  1  =  0  or  1  =  2 
focusing  fields  to  provide  detuning  of  orbital  resonances. 

In  Section  II.  the  stability  properties  of  the  1  =  0  and  1  =  2  configurations  are 
derived,  including  the  effects  of  the  self  fields  of  the  beam.  The  particle  tunes  are 
calculated.  In  Section  III  particle  orbits  for  the  conventional  and  modified 
betatrons  and  the  1  =  0  and  2  stellatrons  are  calculated  numerically  and  compared. 
The  energy-mismatch  bandwidth  of  the  steUatron  is  calculated  as  a  function  of  the 
focusing  strength.  The  sensitivity  of  the  particle  orbits  to  the  betatron  field  index 
in  the  steUatron  is  also  examined. 

In  Section  IV,  integer  orbital  resonances  in  the  stellatrons  are  discussed.  Integer 
orbital  resonances  occur  when  the  betatron  “tune”  (the  number  of  betatron 
wavelengths  around  the  circumference)  is  an  integer.  In  Section  V,  various 

I-  0  SYSTEM  -  -  "BUMPY  TORUS" 


FIGURE  2  The  l  =  0  steUatron. 
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high-current  injection  schemes  are  described  and  the  potential  advantages  offered 
9  by  the  stellatron  are  examined.  Section  VI  is  a  discussion  of  the  available 

literature  (much  of  it  unpublished),  and  Section  VII  presents  a  summary  of  the 
results.  t 


II.  PARTICLE  ORBITS  IN  THE  PARAXIAL  APPROXIMATION 

A  linearized  or  paraxial  analysis  of  particle  orbits,  valid  for  a  particle  that  is 
“near”  a  circular  design  orbit  with  “nearly”  the  correct  energy  to  be  matched  on 
that  orbit,  has  been  carried  out  to  gain  some  quantitative  understanding  of 
particle  behavior  in  the  stellatron.  Orbits  of  particles  not  satisfying  these  condi¬ 
tions  must  in  general  be  found  numerically,  as  discussed  in  the  following  section. 
The  linearized  analysis,  however,  yields  important  information  about  the  fre¬ 
quency  and  the  stability  of  the  particle  oscillations  about  the  design  orbit. 

The  analysis  employs  a  general  magnetic-field  configuration  which  we  describe 
in  the  coordinates  of  Fig.  3A.  The  applied  fields  consist,  in  part,  of  a  vertical  and 
radial  "betatron"  field  of  the  form 

B6  -  -  —  8,0r  +  Bsof  1  -  — (1) 
r0  L  r0  J 

where  (r,  z)  =  (r0,  0)  is  the  location  of  the  design  orbit  (assumed  circular),  Bl0  is 
the  vertical  field  at  the  design  orbit,  and  n  is  the  betatron  field  index.  A  toroidal 
field, 

“-[‘-T*]4  (2) 

is  superimposed  on  (1),  as  well  as  a  “stellarator”  or  multipole  field  which  is 
written  in  the  cylindrical  approximation  as  the  negative  gradient  of  the  scalar 
potential 

<bs(p,  <f>,  s)  =  ~  -t*  2'/| (fcp)  sin  ( l<t>  +  ks),  (3) 

k 

where  k  =  m/r„  and  are  constants  and  l  and  m,  referred  to  as  the  poloidal  and 
toroidal  field  numbers  respectively,  are  taken  as  integers.  In  addition,  s  is  defined 
to  be  -ro0,  so  that  (p.  <t>,  s)  is  a  righthanded  system.  The  axis  of  the  stellarator 
field  is  assumed  to  be  aligned  with  the  symmetry  plane  of  the  betatron  field,  z  =  0. 

Each  of  the  fields  ( 1 )— ( 3)  has  its  own  special  purpose.  The  betatron  field  (1),  of 
course,  acts  simply  to  cancel  the  centrifugal  force  experienced  by  a  circulating 
particle.  In  a  stellatron,  this  field  rises  with  the  particle  energy  so  that  the  betatron 
condition  is  (at  least  approximately)  satisfied.  The  weak-focusing  nature  of  (1), 
when  n  is  between  o  and  l,  which  is  crucial  to  the  successful  operation  of 
betatrons,  is  of  secondary  importance  to  the  stellatron.  Stable  orbits  still  occur  in 
the  stellatron  when  n  is  outside  of  the  interval  (0,  1).  The  stellatron  configuration 
therefore  has  the  virtue  of  being  insensitive  to  errors  in  the  vertical  field  or  its 
gradient.  This  insensitivity  has  beneficial  consequences  in  practical  designs. 
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FIGURE  3  Gcumetrv  and  coordinate  systems.  The  origin  in  Fig.  3A  is  at  the  center  of  a  torus  of 
major  radius  r„.  Fig.  3B  illustrates  a  cut  along  the  minor  cross  section,  showing  the  beam  and 
individual  particle  positions. 


A  toroidal  field  (2)  is  added  to  the  conventional  betatron  in  the  hope  of 
increasing  the  contained  current.  Indeed,  it  has  been  found  that  with  only  fields 
(1)  and  (2)  the  beam  densities  which  can  be  contained  in  equilibrium  increase  by  a 
factor  (BJ2B:)2.  for  large  B„,  over  those  in  a  conventional  betatron,  (1).  Addi¬ 
tionally,  the  toroidal  field  acts  to  control  the  tune  shift  due  to  the  effects  of  space 
charge.  This  feature  is  important  when  orbital  resonances  must  be  avoided. 

The  stellarator  field  (3)  is  added  to  provide  the  beneficial  effects  of  alternating- 
gradient  strong  focusing.  For  l  =  0  or  2,  Eq.  (3)  gives  an  alternating  gradient  field 
component  at  the  design  orbit.  This  field  greatly  improves  the  momentum 
compaction  over  that  of  a  configuration  consisting  of  (1)  and  (2)  only,  and  so 
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greatly  relaxes  certain  design  tolerances  of  the  machine.  Windings  of  other 
(-number  may  also  be  useful;  in  particular  fields  with  i  2:3  may  help  in  controlling 
orbital  resonances,  but  the  analysis  must  be  done  numerically.  The  l  =  1  field 
produces  an  oscillating  magnetic  field  (not  a  field  gradient)  on  the  axis;  such  a 
configuration  would  not  possess  a  circular  design  orbit  and  is  not  expected  to 
possess  any  noteworthy  focusing  properties.  This  paper  is  concerned  only  with 
stellatrons  having  l  =  0  (a  "bumpy  torus")  and  1  =  2  (a  “twisted  quadrupole"). 

Since  the  stellatron  is  intended  to  carry  large  currents,  it  becomes  important  to 
consider  the  effects  of  the  self  fields  of  the  beam.  In  general,  this  is  an  extremely 
hard  problem.  To  make  progress  analytically,  a  simple  model  for  these  fields  is 
assumed,  consisting  of  a  circular  beam  cross  section  (minor  radius  rb)  and  uniform 
density  with  center  located  at  (rn  +  Ar,  Az)  in  a  perfectly  conducting  cylindrical 
chamber  of  radius  a  (Fig.  3B).  The  self  fields  in  the  stellatron  at  the  particle 
location  (r0  +  Ar  +  Sr,  Az  +  Sz)  are  given  by 


Er  =  -2t Tn0e{ 

>4-) 

(4a) 

Ee  =  0 

(4b) 

£z  =  -27rn0e( 

>2+^42) 

(4c) 

Br  =  0o  Et 

(4d) 

B9  =  0 

(4e) 

Bz  —  ~&oEr 

(4f) 

in  cgs  units,  where  n„  is  the  beam  density,  -e  is  the  electron  charge,  and  0O  =  u0/c 
is  the  velocity  of  a  particle  on  the  design  orbit  normalized  to  the  speed  of  light.  In 
Eqs.  (4a-f),  the  first  term  in  each  set  of  parentheses  is  due  to  the  self  field  of  the 
beam,  while  the  second  term  is  due  to  the  image  of  the  beam  in  the  perfectly 
conducting  wall.  Toroidal  corrections  to  Eqs.  (4)  may  be  important  in  some  cases 
but  are  ignored  here. 


( =  2  Stellatron 

The  case  l  =  2  is  considered  first  in  detail  since  this  case,  rather  than  I  =  0, 
corresponds  most  closely  to  a  conventional  quadrupole  strong-focusing  system. 
The  /  =  2  stellarator  field  near  the  axis  is,  from  (3) 

B1  =  z  |  cos  md  -  rt  sin  md  ]f  +  i  rt  cos  md  +  z ,  sin  md  |z},  (5) 

where  r,  =  r-r„,  z,  =  z.  Taking  the  fields  (1),  (2),  (4),  (5)  and  using  them  in  the 
equations  of  motion  gi\es  equations  correct  to  first  order  in  the  small  quantities 
r,/rn,  A r/r„,  z,/r,„  and  Az.'r,, 

1  to2  /  r2  \ 

r,  +ft?„(l-  n  +  a  cos  --  — r  (r,  -  ArH — \  Ar) 

2  To  \  a  / 

=  ftl0Vntffn  —  ( sin  m^z.+fieoz,  (6a) 
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z,+flz0(n-n  cos  m0)zt 


e 


2yg\  '  a‘  I 

i  =  o. 


-(^O^osin  md)rl-nm)rl 


(6b) 

(6c) 


where  ft.0  =  eB:Q/mty0c,  n  =  kr0Bl0l Bz0,  ojI  =  4-im0e2lmey0,  mt  is  the  electron 
mass,  and  u#1  =  r,Oo  +  ro0t.  A  subscript,  0,  refers  to  quantities  evaluated  on  the 
design  orbit  where  vn/r0  =  Q,0  and  dots  denotes  time  derivatives. 

As  they  stand,  Eqs  (6a-c)  are  not  straightforward  to  solve  since  they  involve 
both  the  motion  of  the  particle  about  the  beam  center  and  the  motion  of  the  beam 
center  about  the  design  orbit.  Since  r,  =  Ar  +  Sr  and  z,  =  Az  +  8z,  separate  solu¬ 
tions  are  required  for  (rt,  zt)  and  (A r,  Az).  A  self-consistent  set  of  equations  for 
both  beam-center  motion  and  for  motion  of  particles  about  the  beam  center  may 
be  found,  however,  by  averaging  (6a-b)  over  a  distribution  of  particle  initial 
conditions.  Details  of  this  averaging  procedure  are  given  in  the  Appendix.  The 
result  is  that  motion  of  the  beam  center  is  governed  by  the  equations 


Ar+  ( 1  -  n  n,  +  ^  cos  mO'^Ar  =  C-r  -^-(p.  sin  md')Az  +  b 
dd  \  a~  /  r0flro  y o 

Az  +  [n  -  -p  n,  -  fi.  cos  md'  jAz  =  — (fi  sin  m0') Ar  -  b  — , 


d  Az 
d0' 


(7a) 

(7b) 


where  we  have  changed  independent  variables  from  (0,  t)  to  (0',  t')a(0,  t-6/Slz0), 
i.e..  Eqs.  (7a-b)  describe  the  movement  of  the  beam  center  as  seen  by  a  reference 
particle  moving  with  angular  velocity  Ol0-  In  Eq.  (7),  the  notation  n, » 
wh(2yl(ll0),  b  ■  B90/Bz0,  and  m0oyoPi  is  introduced.  In  Eq.  (7a),  the  qutuitity 
(y,>  is  just  the  ensemble  average  of  the  difference  in  energy  between  a  particle 
and  the  reference  particle,  i.e.  - 7o)-  (fi)  is  independent  of  time  and  for 

a  matched  beam,  <>,>  =  0.  Once  Eq.  (7)  is  solved,  the  solution  may  be  inserted 
into  Eqs.  (6a,  b)  which  may  then  be  solved  for  the  location  of  any  individual 
particle.  The  method  used  for  solving  Eq.  (7)  is  described  below,  and  may  be  used 
on  Eq.  (6)  as  well. 

Equations  (7a,  b)  may  be  completely  solved  in  the  special  case  n  =  1/2,  a  fact 
which  traces  its  origin  to  the  poloidal  symmeti  y  this  value  of  n  imposes  on  the 
restoring  forces  experienced  by  the  beam.  Although  the  solution  is  special  in  this 
sense,  it  is  expected  to  share  most  of  its  features  with  the  solution  for  arbitrary  n 
(which  must  be  obtained  numerically)  since  the  particle  focusing  should  be 
dominated  by  the  quadrupole  strength  (a  and  only  be  weakly  affected  by  the 
precise  value  of  n.  For  n  =  1/2,  then,  define  the  complex  variable  4  =  (Ar  +  iAz)/r0 
for  which  Eqs.  (7a,  b)  give  the  equation 

€"  +  ib f  +  {\  -  P  +  =  <A>  (8) 

where  <A>  is  the  momentum  mismatch  (yi)/((3oyo)-  In  Eq-  (8),  primes  denote  dldO' 
and  an  asterisk  denotes  complex  conjugate.  The  further  change  of  variable  |  = 
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ipe‘m<> 72  yields  an  equation  with  constant  coefficients, 

(i  2  mb  yy\ 

— 2 - —  +  =  Ae',m®72,  (9) 

C 

the  solution  to  which  consists  of  a  particular  solution  plus  a  sum  of  exponentials  of 
the  form  4>oe7C"'-a  where  i //0  is  a  constant  and  v*  are  given  by 

v±  =  [n+km2±(nm2  +  p,2)1'2!1'2,  (10) 

where  n  =  1/2  -  r2/a2n ,  +  b2l 4,  m  =  m+b.  A  particular  solution  to  Eq.  (9)  is 

4f„  =  Ae‘m9V2  +  Be-‘me72,  (11) 

where 


A  = 


B  = 


fi(A> 


_ <A> _ 

±-^  +  *2[m2  +  m6-±+^n,]  ‘ 


(12) 

(13) 


Referring  back  to  the  definitions  of  £  and  i>,  it  is  evident  that  B  gives  the 
“zero-frequency”  part  of  the  radial  shift  due  to  a  beam  with  momentum  mismatch 
(A),  that  is,  B/(A >  is  the  momentum  compaction  factor,  which  may  be  made  small 
by  choosing  m-  large.  Setting  p.  =  0  in  Eq.  (13)  recovers  the  usual  result  for  a 
conventional  betatron  or  modified  betatron.  (Addition  of  a  toroidal  field  to  a 
conventional  betatron  does  not  affect  momentum  compaction.) 

The  small  oscillations  of  the  beam  center  are  stable  ( v±  are  real)  if  and  only  if 
three  conditions  are  satisfied 


n+km2>  0 

(14a) 

(n  -\m2)2- vl2>0 

(14b) 

nm2  +  p.2>  0. 

(14c) 

These  conditions  are  summarized  in  Fig.  4,  where  the  stable  region  of  parameter 
space  is  illustrated  in  terms  of  the  auxiliary  variables  u  =  4n/m2  and  v  =  ]/x]/m 
Typically  one  operates  in  the  stable  region  with  u<l.  The  orbital  stability 
requirement  places  strict  limits  on  the  size  of  n  that  can  be  tolerated,  i.e.,  strict 
limits  on  the  achievable  momentum  compaction  for  given  b  and  m.  Any  experi¬ 
ment  must  be  designed  so  that  u  and  v  remain  in  the  stable  region  for  all  time:  a 
sample  experimental  trajectory  is  illustrated  in  the  figure.  For  a  low  current  beam 
n  is  positive  and  Eq.  (14b)  gives  the  only  nontrivial  condition 

\km2  +  mb- l|>|2n|  (15) 

The  orbital  stability  criteria  for  individual  particle  motion  within  the  beam  are 
of  identical  form  to  Eqs.  (14a-c)  with  the  replacement,  in  the  definition  of  n,  of 


i 
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FIGURE  4  The  l  =  2  stability  plane,  u  and  u  are  defined  in  the  tex'  The  arrow  indicates  a  sample 
experimental  trajectory. 


(r,Ja)2n,  — »  n,.  That  this  is  so  may  be  seen  by  an  analysis  of  Eqs.  (6a,  b)  exactly 
parallel  to  that  carried  out  for  Eqs.  (7a,  b)  considering  that  now,  in  Eqs.  (6a,  b),  A r 
and  A z  are  known  quantities.  With  this  substitution  in  Eq.  (13),  the  zero- 
frequency  term  in  the  individual  particle  motion  is  obtained  as 

a  *-<*>  (16y 

p  ^-n,  +  fx2(m2  +  mb-|  +  n,)  1  ’ 

where  A  =  yJyQPl  is  proportional  to  the  energy  difference  between  the  individual 
electron  energy  and  the  matched  beam  energy.  With  p.  =  0,  the  results  for  a 
conventional  betatron  and  a  modified  betatron  are  obtained.  This  solution  with 
pL  =  0  suffers  a  loss  of  momentum  compaction,  i.e.,  Bp  — ►  =»,  when  n,  passes 
through  1/2.  In  a  conventional  betatron,  n,  s  1/2  is  required  for  orbital  stability. 
The  modified  betatron,  however,  would  typically  use  n, » 1/2  at  injection,  and 
would  then  pass  through  n,  =  1/2  during  acceleration,  since  n,  «  y_J.  This  transi¬ 
tion  at  n,  =  1/2  corresponds  to  a  momentary  loss  of  equilibrium  in  the  modified 
betatron  as  the  configuration  switches  from  a  diamagnetic  to  a  paramagnetic 
equilibrium.2'’ 27 

This  transition  can  be  avoided  in  the  stellatron  by  using  large  |jjl|.  The  stability 
condition  ( 14b),  however,  sets  an  upper  limit  to  |jx|  for  a  particular  choice  of  m,  b 
and  n,.  The  requirement  on  (i!  may  be  expressed  as 

(ns  —  l/2)(m2  +  mb  +  n,  —  1/2)  <  p2  <(n,  —  1/2+  m2/4+  mb/2)2.  (17) 


For  the  condition  (17)  to  be  meaningful,  n,  must  satisfy 


n. 


1  1 
2<2 


(18) 


which  is  clearly  compatible  with  n,  »  1/2.  The  l  =  2  stellatron,  therefore,  can  avoid 
a  loss  of  momentum  compaction  during  acceleration,  while  retaining  the  ability  to 
use  large  n,  at  injection. 
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1  =  0  Stellatron 

The  I  =  0  stellatron  has  fields  given  by 

*  sin  mfl 

(19a) 

BJ^B.qCOs  m8 

(19b) 

B’t^zBtomy  sin  md 

(19c) 

in  a  cylindrical  (r.  8,  z)  coordinate  system  (Fig.  3),  where  x  =  rjr0  and  y  =  z,/r0. 
Using  the  procedure  described  above  for  linearizing  the  equations  of  motion  in 
the  displacement  from  a  reference  orbit  and  using  the  paraxial  approximation  for 
the  electron  motion  with  n  =  1/2,  the  equations  of  motion  for  single  elect¬ 
rons  including  beam  self-fields  may  be  expressed  as 

2(2  — 4n,  +  tr(l  +  e  cos  20m)2Ii/<  =  — j^pexp  [-^-(20m  +  e  sin20m)  1  (20) 
m  m  P0  L2m  J 

where  e  =  Bs0/B»o ,  =  mdl 2,  P0  is  the  momentum  of  a  particle  which  is  matched 

on  the  reference  orbit,  SP  is  the  “momentum  error”,  and 

<1*  =  (x  +  iy )  exp  ^  |b(l  +  e  cos  m8)  dd j.  (21) 

In  Eq.  (20),  n,  will  in  general  depend  on  6  in  a  way  governed  by  an  envelope 
equation,  thereby  making  the  inclusion  of  self-field  effects  more  difficult  than  for 
the  1  =  2  case  in  which  a  constant  beam  radius  can  be  maintained.57  Consequently, 
only  single-particle  orbits  will  be  considered. 

Equation  (20)  is  a  Hill  equation,  for  which  there  are  theorems  concerning 
characteristic  frequencies  and  stability.44  The  Floquet  solutions  to  this  equation 
for  n,  =  0,  are  of  the  general  form 

oo 

4,  =  e"*m  I  Cnei2n\  (22) 

rt  -  -00 

which  display  an  infinite  set  of  natural  oscillations  .at  the  frequencies,  v±2n.  The 
value  of  v  can  be  determined  numerically  for  a  given  set  of  accelerator  parame¬ 
ters,  but  the  dependence  of  v  on  these  parameters  is  not  known  analytically. 

It  is  possible  to  determine  and  display  the  regions  of  parameter  space  where  v 
is  real,  i.e.,  where  the  motion  is  stable.  Starting  with  the  homogeneous  equation, 

<T  +  -4l2  +  f>2(l  +  e  cos 26m)2 1^  =  0,  (23) 

m 

the  solution  1 Ii(tt),  which  satisfies  «A(0)  =  1  and  >1/(0)  =  0  is  constructed  numerically. 
The  stability  condition44  is  |<Mrr)|  <  1,  which  is  illustrated  in  Fig.  5  for  a  particular 
example.  The  intersections  of  the  unstable  regions  with  the  abscissa  are  given  by 

— 2  (b2  +  2)  =  q2;  q  =  l,2,...  (24) 

m 

which  correspond  to  resonances  between  the  “focusing  frequency”  which  a 
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/- 0 


B  $/m  Bx 

FIGURE  5  The  1  =  0  stability  plane.  The  shaded  regions  are  unstable. 


particle  experiences  and  its  cyclotron  frequency  in  the  toroidal  field.  Under 
resonant  conditions,  there  is  a  transfer  of  energy  from  longitudinal  to  transverse 
motion  and  the  adiabatic  invariance  of  the  magnetic  moment  in  the  toroidal  field 
is  destroyed.  During  acceleration,  if  B,  is  increased  and  B#  is  held  fixed,  the 
operating  point  on  Fig.  5  will  move  from  right  to  left.  The  accelerator  should, 
operate  in  the  left-most  stable  band  to  avoid  crossing  unstable  bands.  At 
injection,  therefore,  m  >  b  is  required,  which  implies  the  need  for  a  large  number 
of  field  periods.  When  space-charge  effects  are  taken  into  account,  one  finds  the 
stronger  requirement,  m  >2b,  to  avoid  envelope  instabilities.64 

The  momentum  compaction  factor  for  the  1  =  0  configuration  may  be  calcu¬ 
lated,  for  small  values  of  e  and  b/m,  directly  from  Eq.  (20).  For  n,  =0,  the  result 
is 


5x  -Jl  (- 

■by  i  i 

(bp/po)  L  ' 

2/  l-(fe/m)2J 

which  illustrates  the  helpful  effect  of  the  alternating  gradient,  as  measured  by  e, 
on  the  tolerance  of  the  1  =  0  device  to  momentum  mismatch.  When  e  or  b/m  is 
large,  numerical  integration  of  the  orbit  equations  is  required. 


III.  NUMERICAL  CALCULATION  OF  SINGLE-PARTICLE 
ORBITS 

A  computer  code,  which  integrates  the  single-particle  equations  of  motion,  has 
been  utilized  to  study  nonlinear  aspects  of  the  stellatron  configuration.  The  code 
solves  the  fully  relativistic  dynamical  equations,  without  utilizing  the  paraxial 
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approximation  for  the  electron  motion  or  an  expansion  in  the  particle  displace¬ 
ment  from  a  reference  orbit.  The  applied  fields  in  the  code  include  toroidal 
corrections  to  first  order  in  the  inverse  aspect  ratio.  In  addition,  the  code  does  not 
assume  the  fi^Jd  index  to  be  1/2. 

The  total  magnetic  field  used  by  the  code  is  the  superposition  of  a  conventional 
betatron  field  Bb  and  a  stellarator  field  B,.  The  conventional  betatron  field  is 
given  in  Eq.  (1).  The  stellarator  field  may  be  expressed  as 

B,  =  -Vd>,  (26) 


where  <t>,  is  the  magnetic  scalar  potential  for  the  stellarator  field.  To  first  order  in 
the  inverse  aspect  ratio,  4>,  has  been  given  by  Danilkin45  as 

<f>,  s)  =  -B90|  s  +  —  l,(x)  sin  [  l<t>  -  las  | 
l  a 

-  ~2  (x2/',(x)  -  x(l-hl)f,(x)lsin  1(1 las  I 
4/o 


[  JT/;(x)  -  x(  1  -  /)/, (x)  sin  | (J  -  l)<f>  -  las  ] 


K£( 

Ki  +  i(x0)  —  x0IC!+i(xo) 

4/a2 

L  k;(x0) 

KB, 

Ki-i(xo)  ~  x0Ki-i(xo) 

4/a2 

L  K((x  0) 

where  s  =  ro0,  x  =  lap,  a  =  2-rr/L,  L  is  the  helix  pitch  length,  k  =  l/r0,  x0  =  lap,., 
P  =  pc  is  the  location  of  (assumed)  wall  surface  currents,  and  I,  K  are  the  modified 
Bessel  functions.  The  coordinates,  (p,  <f>,  s)  are  defined  in  Fig.  3.  More  generally, 
the  field  may  be  expressed  as  a  superposition  of  fields  of  different  /-number.  This 
feature  has  not  been  included  in  the  numerical  model.  In  addition,  beam  image 
forces  may  be  included  in  the  model,  but  have  not  been  utilized  for  any  of  the 
results  presented  in  this  section.  The  wave  number,  k  =  m/r0,  in  the  preceding 
analytical  model  is  equivalent  to  -la. 

Figure  6  shows  single-particle  orbits  obtained  with  this  code.  The  figure  is  the 
projection  of  the  orbits  onto  the  minor  cross-section.  The  torus  has  a  lm  major 
radius  and  0.1m  minor  radius.  The  vertical  magnetic  field  is  fixed  at  Bl0  =  118  G, 
corresponding  to  a  matched  particle  (at  the  minor  axis)  with  y„  =  7.  The  individual 
frames  of  Fig.  6  show  a  comparison  of  particle  orbits  in  a  conventional  betatron,  a 
modified  betatron,  an  l  =  0  stellatror  and  an  l  =  2  stellatron.  Each  frame  shows 
two  particle  trajectories,  one  with  positive  mismatch  and  one  with  negative 
mismatch.  All  the  trajectories  are  initialized  on  the  minor  axis  with  momentum 
parallel  to  the  minor  axis,  and  with  various  amounts  of  energy  mismatch,  Ay/y0, 
as  shown  under  each  frame.  The  conventional  and  modified  betatron  frames 
demonstrate  that  these  devices  cannot  tolerate  mismatch,  Ay/y0  —  ±3%.  In  fact, 
using  Eq.  (12),  we  can  estimate  the  tolerance  of  these  devices  to  energy  mismatch. 
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PARTICLE  ORBITS 


±3%  MISMATCH  ±3%  MISMATCH 


-0.1  ON  0.1 
i-  0 "BUMPY TONUS* 


±20%  MISMATCH 


1-  2  "STELLATRON* 
±20%  MISMATCH 


FIGURE  6  Particle  orbits  in  conventional  betatron,  modified  betatron,  I  =  0  stellation,  and  1  =  2. 
stellatron. 

For  n,(rb/a)2«  1/2,  the  allowed  energy  mismatch  is  given  by 


Ay  la 

—  -s - . 

yo  4r0 


(28) 


For  a/r0  =  0.1,  we  obtain  Ay/y0s2.5%,  which  is  consistent  with  the  figure.  The 
l  =  0  and  /  =  2  stellatron  configurations,  on  the  other  hand,  can  retain  particles 
with  more  than  Ay/y0  =  ±20%  mismatch,  as  shown  on  the  lower  two  frames  of 
Fig.  6. 

Ay 
r 
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RGURE  7  Energy  bandwidth  as  a  function  of  focusing  strength  e  for  the  I  =  0  stellatron.  (n  =  1/2). 
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FIGURE  8  Energy  bandwidth  as  a  function  of  focusing  strength  e  for  the  1  =  2  steilatron.  (n  =  1/2). 


The  tolerance  of  the  accelerator  to  mismatch  between  the  average  beam  energy 
and  the  vertical  magnetic  held  can  be  described  as  a  bandwidth  of  allowed 
mismatch.  For  larger  mismatch,  the  beam  excursion  from  its  reference  orbit  does 
not  ht  inside  the  accelerator  aperture.  The  bandwidth,  therefore,  is  essentially  a 
measure  of  the  momentum  compaction  of  the  accelerator.  Since  a  calculation  of 
the  bandwidth  involves  particles  with  large  Ay/y0  that  make  large  excursions  from 
the  reference  orbit,  the  analytical  formalism  of  the  preceding  section  is  inade¬ 
quate. 

By  launching  particles  as  shown  in  Fig.  6,  with  various  Ay /y0,  the  largest  |Ay/y0| 
for  which  the  particle  is  contained  can  be  found.  This  bandwidth  can  be  displayed 
as  a  function  of  any  parameter  of  the  accelerator.  Figures  7  and  8,  for  example, 
show  the  bandwidth  for  the  l  =  0  and  1  =  2  stellatrons  as  the  amplitude  of  the 


n 

FIGURE  9  Energy  bandwidth  as  a  function  of  vertical  field  index  n  at  e  =  0.3  for  the  1  =  0  steilatron. 
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FIGURE  10  Energy  bandwidth  as  a  function  of  vertical  field  index  n  at  e  **  0.3  for  the  1 ” 2 
Stella  tron. 

strong  focusing  field  is  varied.  In  both  devices,  the  bandwidth  is  ±2.5%  in  the 
limit  where  e  =  0,  and  exceeds  20%  for  e  >  0.3.  For  l  =  2,  the  quantity,  ji,  in  the 
preceding  section  is  given  by  n  =  emb/2,  where  b  =  B90JBz0. 

Figures  9  and  10  show  the  bandwidth  for  the  1  =  0  and  1  =  2  stellatrons  at 
e  =  0.3  as  the  vertical  field  index  is  varied.  The  results  here  demonstrate  the, 
relative  insensitivity  of  the  stellatron  configurations  to  the  field  index.  The  need  to 
maintain  good  uniformity  in  the  field  index  for  a  conventional  or  modified 
betatron  drives  the  stored  energy  and  cost  for  the  vertical  field  coils  (or  leads  to 
the  use  of  iron)  in  those  accelerators.  The  stellatron  configuration  greatly  relaxes 
the  need  for  such  uniformity. 


IV.  ORBITAL  RESONANCES  IN  THE  STELLATRON 

Until  now,  the  applied  fields  in  the  stellatron  have  been  assumed  to  be  “perfect" 
in  the  sense  that  only  the  fields  described  analytically  above  in  Eqs.  (1-3)  are 
present.  In  general,  of  course,  a  physical  magnet  system  will  have  some  imperfec¬ 
tions  leading  to  small  field  and  focusing  errors,  which  can  be  represented  by 
additional  terms  in  the  paraxial  equations.  Generally  these  small  errors  lead  to  a 
small  response  but  under  certain  conditions — when  some  rational  number  of 
particle  wavelengths  fit  around  the  machine — the  particle  orbits  can  be  violently 
disrupted.  In  this  section,  these  resonance  conditions  are  considered  for  the  l  =  0 
and  l  =  2  stellatrons,  and  possible  ways  either  to  avoid  them  or  to  minimize  their 
effects  are  discussed.  It  must  be  said  however,  that  there  is  no  reason  why  such 
resonances  should  be  of  any  less  concern  in  the  modified  betatron  and  stellatron 
class  of  devices  than  they  are  in  more  conventional  cyclic  accelerators  in  which 
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they  are  of  course,  of  crucial  importance  in  establishing  an  operating  point. 
Though  some  of  the  high  mode-number  integer  resonances  might  be  successfully 
passed  through,  there  is  considerable  reason  to  believe  that  major  beam  disrup¬ 
tions  may  occtir  at  low-order  resonance  crossings,  if  such  are  attempted.  This 
point  is  not  often  sufficiently  stressed  in  the  literature  promoting  the  use  of 
devices  similar  to  some  of  those  employed  in  magnetic-confinement  fusion 
research  to  accelerate  large  currents  to  high  energies. 

This  section  will  focus  on  the  integer  resonances,  driven  by  the  Fourier 
components  of  an  error  in  the  vertical  field,  though  others  may  also  be  important. 
The  tune  shift  due  to  space  charge  will  also  be  neglected  at  first,  to  consider  the 
resonant  response  of  single  particles. 

As  in  the  previous  section,  the  l  =  2  case,  for  which  the  tunes  may  be  explicitly 
calculated  (Eq.  (10)1,  is  considered  first.  For  single-particle  motion,  the  four 
possible  “betatron”  frequencies  are  ml 2  ±  v*  where  v±  satisfies  the  biquadratic, 

u*-(\m2+mb  +  b2+  D^  +  HItn2-)-  mb  -  l)2-p2  =  0.  (29) 

The  most  general  resonance  condition  is  written 


Im 

\  Im  N 

1  /m . 

\  Im  \ 

(r- 

M  1~V-) 

where  n, . n4,  p  are  all  integers;  this  condition  applies  to  all  types  of  resonant 

driving  terms.  Integer  resonances  occur  when  all  n,  except  one,  nk,  vanish  and 
plnk  =  n  is  itself  an  integer,  the  Fourier  mode  number  of  the  field  error.  When  an 
integer  resonance  occurs,  the  displacements  which  a  particle  undergoes  on  succes¬ 
sive  encounters  with  a  field  “bump”  add  in  phase  and  the  result  can  be  loss  of  the 
particle  from  the  machine. 

Figure  11  shows  the  integer  resonance  contours  1,  2, ... ,  10  in  the  b-p  plane 
for  m  =  20.  (The  n  and  m  -  n  contours  coincide.)  The  resonance  contours  are  all 
hyperbolae,  from  Eq.  (29);  the  stability  boundary  coincides  with  the  degenerate 
case  n  =  ml 2.  Ideally,  an  accelerator  should  be  designed  to  avoid  all  resonance 
contours,  a  design  which  may  clearly  be  realized  by  holding  b  and  p  fixed  during 
acceleration  (neglecting  the  tune  shift  due  to  space  charge,  which  may,  in  fact, 
limit  the  current  in  this  as  in  other  cyclic  devices).  If  b  and  p  are  allowed  to 
change  during  acceleration,  it  may  be  possible  to  pass  through  at  least  some  of  the 
integer  resonances,  though  this  is  only  speculation  based  on  a  few  computer 
studies.  The  possibility  of  passing  through  resonances  in  the  stellatron  is  discussed 
further  below. 

The  resonance  situation  in  the  bumpy  torus  is  somewhat  less  straightforward  to 
analyze.  In  the  presence  of  an  error  in  the  vertical  field  having  Fourier  mode 
number  fc,  a  term  of  the  form 

(const)  •  e\p\2ik0Jm  +(ib/2m)(2flm  +  e  sin  20m)]  (31) 

is  added  to  the  right-hand  side  of  Eq.  (23),  where  the  multiplicative  constant  is 
proportional  to  the  magnitude  of  the  field  error.  One  may  employ  the  Green's 
function  constructed  from  the  Roquet  solutions  to  Eq.  (23)  to  deduce  that  a 
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L--2  PLANE 


FIGURE  11  Integer-resonance  contours  in  the  b/m  -  n  plane  for  the  / »  2  ste. 


secular  term,  proportional  to  flm,  will  occur  in  the  particular  solution  if 

cos  7r(~~~)  ~  ^i(w),  (32) 

where  is  the  solution  to  the  homogeneous  equation  satisfying  i/»,(0)  =  l, 
i/r((0)  =  0.  Contours  along  which  Eq.  (32)  is  satisfied  may  be  generated  numeri¬ 
cally  for  each  k  in  the  blm  -e  plane;  an  example  is  shown  in  Fig.  12.  Note  the 
behavior  of  the  k  =  0  resonance  line,  which  crosses  other  fc-lines  at  b/m  =  2-  klm 
in  the  second  stable  band.  A  k  =  0  field  error  is  equivalent  to  a  momentum  error, 
so  one  expects  momentum  compaction  to  be  poor  near  the  it  =  0  contour.  It  is 
unlikely,  however,  that  one  would  try  to  operate  in  the  second  stable  band. 

It  is  clear  that  the  “stable”  regions  of  Fig.  12  are  in  fact  criss-crossed  with 
potentially  dangerous  resonance  lines  which  must  either  be  avoided  or  crossed 
somehow  during  the  course  of  an  experiment.  The  operating  point  of  an  ac¬ 
celerator  may  be  held  fixed  in  the  plane,  of  course,  by  raising  all  fields  in 
synchronism,  i.e..  by  keeping  e  =  8BJBe0  and  b  =  Be0IBz0  fixed  for  all  time.  It  is 
worth  considering  the  feasibility  of  crossing  these  resonances,  however,  since 
increasing  all  the  fields  together  has  definite  large  costs  in  terms  of  field  energy.  It 
might  be  desirable,  for  instance,  to  accelerate  while  keeping  B90  and  8B„0  fixed  in 
time,  allowing  the  tunes  to  change  as  b  does.  The  immediate  question  is,  “can 
resonances  be  passed  through  without  major  disruption  of  the  beam?” 

Since  resonances  occur  when  the  number  of  betatron  wavelengths  completed  in 
one  major  revolution  (i.e.,  the  “tune”)  is  an  integer,  half-integer,  or  other  rational 
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FIGURE  12  Integer-resonance  contours  in  the  b/m  -  e  plane  for  the  I  =  0  steUatron. 


number,  a  possible  solution  may  consisting  of  setting  constant  tune  by  allowing 
the  particle  reference  orbit  to  have  a  variable  major  radius  (i.e.,  to  be  noncircu¬ 
lar).  The  FFAG  accelerator  works  in  this  fashion,  and  achieves  constant  tunes  by 
forcing  “self-similar”  orbits  during  acceleration.  A  conventional  betatron  is  the 
special  case  where  the  self-similar  orbits  also  have  constant  radius.  With  addition 
of  a  toroidal  field  which  varies  as  r”\  however,  there  is  no  nonsynchronous  field 
solution  with  constant  tune,  as  shown  below. 

The  single-particle  tunes  (n,  =  0)  in  a  steUatron  depend  on  the  poloidal  and 
toroidal  field  numbers  (l  and  m),  which  are  assumed  fixed,  as  well  as  on  the  field 
ratios,  b  =  B80/B.0  and  e  -  Bs0/Beo,  at  the  reference  orbit.  In  a  modified  betatron, 
the  tunes  depend  only  on  b.  The  B„  and  Bz  fields  have  the  assumed  dependence 


B9  =  B80(f)  — 

r 

(33a) 

/r-\" 

Bz=Bl0(t)\^)  . 

(33b) 

At  t  =  0,  the  particle  is  assumed  to  be  matched  at  r  =  r0;  hence  (for  a  relativistic 
particle) 

_  me2  y(0) 
f°~  e  BlO(0) 

During  acceleration,  the  particle  is  allowed  to  move  its  matched  radius  to  r(f), 
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satisfying 


or  solving  for  r(t). 


Kf)  = 


me2  y(t) 
e  BMO.t) 
me2  y(t)  j~r(t)1" 

e  Bl0( 0  l  r0  J  ’ 


r(t)  =  r0 


t 


me 2  y(f) 
er0  Bz0(t) 


} 


wo— it> 


If  b  -  constant  at  r  =  r(t)  is  required  to  fix  the  tune, 

■  B.0(r)/r0y- 

Bl0(t)  \r(t)J 


erp  B9o(t) 
me2  -y(t) 


=  constant, 


(35) 


(36) 


(37) 


which  implies  that  B„0  must  increase  synchronously  with  y  during  acceleration.  In 
the  stellatron,  the  requirement  for  fixed  tunes  is  that  both  b  and  e  be  fixed 
separately.  This  demonstration  that  b  cannot  be  fixed  if  the  toroidal  field  is 
asynchronous  with  y,  therefore,  applies  to  both  the  stellatron  and  the  modified 
betatron.  It  is  therefore  not  possible  to  fix  the  tunes  in  a  modified  betatron  or 
stellatron  with  an  asynchronous  toroidal  field  that  varies  as  r-1. 

Another  possibility  for  controlling  the  growth  of  single-particle  resonances  may 
be  to  use  nonlinear  focusing  fields  to  de-tune  the  resonance  by  making  the  tunes  a« 
function  of  the  amplitude  of  the  particle  displacement.  The  Bessel  functions  which 
describe  the  l  =  0  or  1  =  2  transverse  fields,  of  course,  already  provide  such  a 
nonlinearity.  The  use  of  l  =  3  fields  in  conjunction  with  a  modified  betatron  or 
l  -  0,  2  stellatron  provides  a  high-shear  nonlinear  field  which  may  be  of  benefit  for 
the  control  of  resonances. 

Rapid  acceleration  of  the  beam  may  also  provide  a  means  of  avoiding  damage 
due  to  resonances  by  rapidly  accelerating  through  them.  Since  the  betatron 
wavelength  depends  on  the  particle  energy,  a  stellatron  with  sufficient  energy  gain 
per  pass  will  traverse  a  resonance  before  the  beam  can  respond.  This  possibility 
has  been  addressed  in  the  numerical  calculation  shown  in  Fig.  13.  The  first  frame 
shows  a  stable  electron  orbit  with  y  -  49  and  with  no  field  error.  In  the  next 
frame,  a  2%  integer  field  error  (n  =  1)  is  superimposed  on  the  stellatron  fields, 
and  the  unaccelerated  electron  trajectory  is  rapidly  lost  to  the  wall.  The  lower  two 
frames  show  the  effects  of  high-gain  acceleration  through  the  resonance.  At  a  gain 
of  3  MeV/pass,  the  particle  motion  is  essentially  contained. 

This  result  suggests  that  a  high-gain  stellatron  accelerator,  with  rapidly  varying 
fields,  may  be  desirable.  Such  a  device  will  not  only  be  beneficial  for  the  control  of 
single-particle  resonances,  but  also  will  have  a  shorter  acceleration  time  than 
conventional  betatrons,  and  will  therefore  be  less  sensitive  to  several  collective 
instabilities.  Resistive-wall  instabilities,  for  example,  are  ineffective  if  the  acceler¬ 
ation  time  is  ^lOO^sec. 
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RESONANT  ORBITS  IN  THE  8TELLATRON 

FIGURE  13  Effect  of  acceleration  through  a  resonance  on  transverse  particle  motion. 


To  allow  field  penetration  in  a  fast  stellatron,  the  vessel  will  have  to  be  slotted, 
and  acceleration  will  then  occur  primarily  at  the  slots  or  gaps.  The  particles  will 
therefore  be  unavoidably  mismatched  as  they  move  between  slots.  The  stellatron 
focusing  can  tolerate  this  mismatch,  and  is  therefore  compatible  with  this  type  of 
device. 


V.  INJECTION 


Injection  into  a  high-current  cyclic  accelerator  is  a  challenge.  The  toroidal  field 
lines  that  contain  the  space  charge  of  the  beam  must  be  crossed  or  perturbed  to 
get  the  beam  in.  A  number  of  experiments  have  been  carried  out  in  similar 
geometries  using  inductive  charging,18  19  35  magnetic  diverters,46  47  and  drift  injec¬ 
tion.48  A  self-synchronous  scheme,  using  the  current  in  the  cathode  shank  to 
make  a  magnetic  diverter,  has  been  used  to  inject  50  percent  of  a  500-kV  20-kA 
50- (j. sec  beam  into  a  racetrack  torus.47  These  experiments  by  Gilad  et  al.  and 
Benford  et  al.  involve  injecting  the  beam  into  neutral  gas.  A  450-keV  16-kA 
25-fj-sec  beam  has  been  injected  into  a  toroidal  magnetic  field  in  a  hard  vacuum49 
to  obtain  a  trapped  beam  current  of  300  A  and  a  quiescent  equilibrium  for 
20  M-sec,  which  is  approximately  3000  revolutions  around  the  torus. 

Recent  experiments5"  have  been  carried  out  in  which  the  torus  is  filled  using  a 
plasma  gun.  The  spatial  plasma-density  profile  is  controlled  by  adjusting  the 
direction  of  plasma  injection  or  the  timing  between  the  gun  firing  and  the 


100  c.  W.  ROBERSON,  A.  MONO  EL  LI  AND  D.  CHERNIN 

application  of  the  pulse  to  the  cathode.  The  plasma  makes  contact  with  both  the 
cathode  and  the  chamber  wall.  When  the  cathode  is  pulsed  negatively,  a  double- 
layer  plasma  sheath  appears,  and  the  electrons  of  the  cathode  plasma  are 
acceleratedc  in  the  sheath  and  ejected  into  the  plasma.  With  a  1.2-MV  80-pisec 
pulse,  a  current  of  50  to  80  kA  can  be  injected  into  the  plasma.  When  the  ring 
radius  is  held  constant,  the  40-kA  ring  decays  to  approximately  half  of  its  value  in 
400  pi  sec.  When  the  ring  is  compressed,  the  25-kA  ring  current  increases  to 
50  kA  and  is  contained  for  about  4  msec.57 

Recently  an  axial  pinch  scheme  has  been  proposed  for  injecting  electrons  across 
magnetic  field  lines  into  the  NRL  modified  betatron.” 

None  of  the  prior  experience  has  established  a  widely  accepted  solution  to  the 
injection  problem.  There  are,  nevertheless,  several  features  of  the  stellatron  field 
configuration  that  may  be  important  from  the  standpoint  of  designing  an  injection 
scheme.  The  present  discussion  will  be  limited  to  the  l  =  2  stellatron. 

The  stellatron  field  has  a  separatrix  within  which  the  magnetic  surfaces  are 
closed  and  beyond  which  field  lines  run  to  the  outside  world.  Injection  along 
magnetic  field  line,  can  be  used  to  place  a  beam  just  outside  the  separatrix.  A 
fast-rising  coil  may  be  used  locally  to  slip  the  separatrix  over  the  beam,  thus 
trapping  it.  Since  the  separatrix  is  a  null  field,  it  can  be  moved  through  the  beam. 

The  rotational  transform  of  the  particle  orbits  due  to  the  helical  fields  can  be 
utilized  to  move  the  beam  electrons  away  from  an  internal  injector.  Electrons  can 
be  forced  to  miss  the  injector  for  many  revolutions  of  the  accelerator.  Figure  14 
shows  a  beam  injected  at  Ar  =  8  cm,  Az  =  0  moving  in  the  toroidal  direction  with 
no  transverse  motion.  The  beam  is  injected  with  I  =  10  kA,  7  =  7  into  fields. 


HGURE  14  Location  of  beam  center  after  injection. 
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B90  =  5  kG,  e  =0.5,  Bz0  =  118  G,  n  =  0.5.  The  analytical  formalism  developed  in 
Section  2  is  used  to  track  the  beam  center  around  the  toms.  The  points  marked 
2 it,  4i r,  6 it,  ...  on  the  figure  represent  the  location  of  the  beam  center  after 
completing  1, 2, 3, . . .  transits  of  the  toms.  Approximately  nine  revolutions 
(200  jisec)  are  required  to  bring  the  beam  electrons  back  to  the  vicinity  of  the 
injector.  The  open  circles,  filled  circles  and  x’s  on  the  figure  represent  the 
injection  of  beams  having  Ay/y0  =  0.1%  mismatch,  and  ±10%  mismatch.  The 
injection  dynamics  are  very  insensitive  to  this  mismatch.  During  the  approxi¬ 
mately  200ns  required  for  the  beam  to  return  to  the  injector,  the  focusing  fields 
must  be  changed  to  trap  the  electrons.  A  small  change  of  the  helical  field  (e.g., 
variation  of  e)  during  injection  may  be  effective  in  trapping  the  beam. 


VI.  DISCUSSION  OF  EARLIER  WORK 

During  the  course  of  this  work,  a  number  of  people  have  pointed  out  related 
unpublished  work.  In  particular,  a  report  on  a  plasma  betatron  that  has  tested  a 
field  configuration  similar  to  the  1  =  2  stellatron  has  been  recently  brought  to  our 
attention.53  The  device  consists  of  a  betatron  field  together  with  toroidal,  helical, 
and  multipole  field  capability  to  study  plasma  containment  configurations.  The 
addition  of  a  betatron  field  to  the  system  is  motivated  by  the  production  of 
runaway  electrons  observ  d  in  the  ORMAK  tokamak.15  Acceleration  is  obtained 
by  changing  the  flux  in  a  central  solenoid.  The  field  at  the  equilibrium  orbit  from 
this  solenoid  is  nulled  to  zero  by  compensating  coils.  A  separate  set  of  programm¬ 
able  vertical-field  coils  is  used  to  provide  the  equilibrium.  One  of  the  reasons 
given  for  driving  the  vertical  field  coils  separately  is  that  the  2-to-l  flux  rule  is  not 
valid  for  currents  comparable  to  the  Alfven  current,53  due  to  self-field  effects. 

The  equilibrium  beam  production  in  a  plasma  betatron  is  complex.  The 
energetic  beam  is  composed  f  runaway54  electrons.  The  current  due  to  those 
electrons  that  do  not  run  away  is  generally  an  order  of  magnitude  higher  than  the 
runaways.  The  vertical  field  is  set  to  provide  an  equilibrium  for  the  plasma 
currents.  The  interpretation  of  the  limiting  energy  of  the  plasma  betatron  in  terms 
of  the  work  presented  in  this  paper  is  as  follows.  The  plasma  currents  provide  a 
rotational  tranform  that  contains  the  high-energy  particles.  This  transform  can  be 
converted  into  an  effective  energy  bandwidth.  As  the  energy  of  the  particles 
increases,  the  mismatch  between  the  parallel  energy  and  vertical  field  increases. 
Since  the  vertical  field  must  be  programmed  to  maintain  the  plasma-current 
equilibrium,  this  mismatch  will  grow  until  the  beam  orbit  size  exceeds  the 
dimensions  of  the  chamber.  Hence,  the  maximum  energy  in  a  plasma  betatron  will 
be  the  bandwidth  provided  by  the  plasma  current  and  the  energy  of  the  equilib¬ 
rium  orbit  due  to  the  vertical  field.  Runaway  electrons  with  energies  of  a  few 
megavolts  are  not  unusual  in  tokamaks,  and  a  1.5-MeV  beam  has  been  produced 
in  a  plasma  betatron."  In  the  Livermore  plasma-betatron  experiment,  the  beam 
energy  is  several  hundred  kilovolts  and  the  beam  current  is  about  10  A.  The 
plasma  cu>rent  is  H>kA. 

The  complexitv  of  forming  a  beam  from  runaways  makes  the  interpretation  of 
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the  Livermore  results  ambiguous.  With  the  helical  field,  the  observed  runaway 
production  is  approximately  10  times  more  intense  and  occurs  at  a  much  lower 
pressure,  indicating  that  the  helical  field  configuration  confines  single  particles, 
but  that  an  azimuthal  current  is  necessary  to  continue  the  runaways  in  the  absence 
of  a  helical  field.  Without  the  helical  field,  the  x-ray  burst  from  the  beam  hitting 
the  wall  occurs  at  approximately  200  nsec.  With  the  helical  field,  the  x-ray  burst 
occuis  at  800  fisec.  The  observed  energy  is  only  200  keV,  however,  which  is  omy 
a  fraction  of  what  should  be  obtained  if  acceleration  occurs  throughout  the  pulse. 

R.  Moir  has  sent  results  from  interesting  work  that  used  an  E  x  B  injector  in  a 
racetrack  geometry.56  The  trajectories  of  an  1  =  2  helical  field  combined  with  a 
vertical  field  on  one  of  the  bends  of  a  racetrack  geometry  is  studied  at  low 
energies.  Resonant  and  nonresonant  diffusion  of  very  low  current  electron  beams 
has  been  studied.  The  beam  makes  100  transits  when  the  stellerator  winding  is  on, 
which  is  about  a  factor  of  5  better  than  without  the  stellarator  windings. 

The  stability  plane  for  continuously  rotated  magnetic  quadrupoles  in  a  linear 
transport  system  has  been  obtained  previously.57  A  stability  plane  similar  to  Fig.  4 
results. 

In  a  series  of  early  experiments58  at  the  Naval  Research  Laboratory,  D.C.,  de 
Packh  has  combined  a  betatron  and  solenoidal  lens  fields.  The  excessive  currents 
required  to  raise  the  solenoid  lens  synchronously  so  as  not  to  cross  resonances 
causes  excessive  heating.  He  has  also  tested  a  combination  betatron  and  air-core 
quadrupole  system  that  can  be  programmed  so  that  no  resonance  are  crossed. 
With  a  55  keV  injection  voltage  he  has  injected  and  accelerated  a  1/2  amp  beam 
to  2  MeV. 

J.  D.  Lawson  has  pointed  out  some  early  work  on  particle  orbits  in  a  betatroa 
with  a  toroidal  magnetic  field.59  The  coupling  between  the  radial  and  vertical 
betatron  oscillations  due  to  the  toroidal  field  is  examined.  This  work  has  been 
motivated  by  some  experiments  in  which  the  addition  of  an  8  gauss  toroidal 
magnetic  field  to  a  conventional  betatron  reduces  the  output  by  75%,  largely 
because  the  injected  electrons  follow  the  toroidal  field  lines  and  intercepted  the 
injector  structure  after  one  revolution. 

In  recent  experiments  at  Cornell  University60  an  electron  ring,  confined  in  a 
cold,  partially  ionized  hydrogen  plasma  has  been  transported  into  a  modified 
betatron  field  configuration  and  accelerated.  The  energy  of  the  beam  is  increased 
from  1.1  to  3.3  MeV.  Although  there  is  no  direct  evidence  of  instabilities,  the  ring 
current  decreases  from  3  kA  to  1.5  kA  during  the  acceleration. 

Recently,  a  periodic  magnetic-focusing  system  for  a  high-current  cyclic  system 
has  been  suggested.61  The  toroidal  field  consists  of  a  series  of  magnetic  cusps.  The 
paraxial  equation  for  this  configuration  may  be  obtained  from  Eq.  (20)  taking 
b  — *  «  and  e  — ►  0  with  be  fixed.  Inje"tion  can  take  place  on  the  open  field  lines 
and  the  nonlinear  focusing  of  the  cusp  may  give  a  sufficient  tune  shift  to  detune 
the  single-particle  resonances. 

VII.  SUMMARY 

The  addition  of  continuous  strong  focusing  fields  to  a  betatron  or  modified 
betatron  leads  to  a  new  configuration  that  has  a  large  tolerance  to  mismatch 
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between  the  particle  energy  and  the  vertical  magnetic  field,  thereby  relaxing  both 
mechanical  and  electrical  design  requirements  on  the  accelerator  and  the  injector. 
The  allowed  mismatch  bandwidth  obtained  by  strong  focusing  is  especially 
attractive  fore  high-current  accelerators,  and  is  compatible  with  operation  at  high 
energy  gain  per  pass. 

Strong  focusing  also  offers  the  advantage  of  introducing  a  threshold,  yT,  for  the 
onset  of  the  negative  mass  instability.  For  y  <  yT,  there  is  no  negative-mass 
instability.  The  threshold  is  yT  =  13  for  typical  stellatron  parameters.40  41  The 
addition  of  stellarator  windings  to  a  betatron  has  not  yet  been  analyzed  in  detail 
for  its  effect  on  collective  instabilities.  Other  issues  which  are  as  yet  unresolved 
include  the  evaluation  of  methods  for  avoiding  orbital  resonances  and  the 
demonstration  of  a  beam  injection  technique. 

The  strong  focusing  associated  with  the  stellatron  allows  a  configuration  to  be 
designed  with  n,  »  1/2  at  injection  which  will  operate  without  passing  through  a 
disruptive  diamagnetic-to-paramagnetic  equilibrium  transition26-27  during  acceler¬ 
ation.  The  modified  betatron,  on  the  other  hand,  suffers  an  instantaneous  loss  of 
momentum  compaction  as  n ,  passes  through  the  value,  1/2,  during  acceleration. 

Of  the  various  /-number  stellatrons,  only  the  1=0  and  1  =  2  configurations 
provide  a  finite  transverse  field  gradient  at  the  beam  axis.  Both  configurations 
allow  a  large  mismatch  between  the  average  beam  energy  and  the  equilibrium 
beam  energy  for  the  applied  vertical  field.  Both  are  also  relatively  insensitive  to 
the  vertical  field  index.  The  1  =  0  system  appears  easier  to  construct,  but  stability 
during  acceleration  requires  the  use  of  a  large  number  of  field  periods,  i.e.  m^b. 
The  /  =  2  configuration,  on  the  other  hand,  requires  a  more  complicated  (quad¬ 
ruple)  field,  but  does  not  require  a  large  number  of  field  periods  for  stable 
operation;  the  field  can  be  generated  using  modular  coils.62-63 

An  objective  of  this  study  has  been  to  describe  a  multi-kiloampere  cyclic 
electron  accelerator  concept  which  is  compatible  with  injectors  based  on  Marx 
pulseline  technology.  Such  injectors  are  limited  to  electron  energy  of  a  few  MeV, 
and  the  accelerator  must  therefore  tolerate  significant  beam  self-fields  at  injec¬ 
tion.  These  fields  can  lead  to  substantial  beam  emittance  as  well  as  beam 
mismatch  with  the  vertical  magnetic  field.  In  addition,  the  Marx  pulseline  technol¬ 
ogy  itself  is  limited  in  the  voltage  ripple  and  the  shot-to-shot  voltage  reproducibil¬ 
ity  that  can  be  achieved  at  reasonable  cost.  The  accelerator  must  be  able  to 
tolerate  such  uncertainties  in  the  injected  beam  parameters  on  each  shot.  The 
stellatron  configurations  described  here  have  a  large  tolerance  both  to  beam 
emittance  and  to  beam  mismatch,  which  makes  this  concept  attractive  for  handl¬ 
ing  high  currents  and  for  mating  to  a  Marx  pulseline  injector. 

As  this  paper  was  being  completed,  encouraging  experimental  results  from  a 
group  at  the  University  of  California  at  Irvine  were  reported  to  us.66  A  200-A 
beam  has  been  accelerated  to  2  MeV  in  a  small  /  =  2  stellatron  device;  smaller 
currents  have  been  accelerated  to  4  MeV  in  the  same  experiment. 
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APPENDIX 

This  Appendix  discusses  the  method  of  averaging  Eqs.  (6a-6c)  over  particle  initial 
conditions  so  as  to  find  separate  equations  for  the  particle  position  (rt(r),  zt(r)) — 
"Lagrangian”  variables,  dependent  on  the  individual  particle’s  initial  position  and 
momentum,  and  for  the  location  of  the  beam  center,  (Ar(0,  t),  Az(0, 0) — 
“Eulerian”  or  field  variables  dependent  only  on  time  t  and  the  observation  point 
6.  The  radial  coordinate  is  treated  here;  identical  expressions  hold  for  z. 

The  distribution  function  may  be  defined  in  complete  generality  in  “Klimon- 
tovich”  form65 


fir ,  9,  z,  p„  p*p2;t)=  J  dV°>  dy0)rV0>,  p(0)) 

■  S<3)(r-r(r<0),  p(0),  0) 

•5<3>(p-p(i<0\p<0\  0),  (A-f) 

where  r(0)  and  pl0)  are  the  initial  position  and  momentum  of  a  particle,  r  and  p  are 
the  solutions  to  the  equations  of  motion  for  a  particle  of  given  initial  position  and 
momentum,  and  the  integral  extends  over  all  initial  conditions.  At  a  given 
azimuth  9,  then 

I  r  dr  dz  d3p(r-  r0)/(r,  9,  z,  p„  pe,  pz ;  t) 

A r(0,  t)  = - - - 

J  r  dr  dz  d3p/(r,  9,  z,  p„  pe,  pz ;  t) 

J dV0'  d3pl0)f ‘’'(r'0’,  p<01)r,  5(0  -  0) 

| d3r<0)  d3pmrV°\  P(0>)  8(9  -  9) 

s<r,),  (A-2) 

where  rt  =  r  -  r„  and  the  arguments  of  r,  and  0  have  been  suppressed.  Next 
consider  the  quantity  {/•,)  and  relate  it  to  derivatives  of  Ar(0,  t).  Defining  for 
convenience  the  operator  f  as 
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one  finds 

at 

<'ri'=  f  ■  8(8-9) 

ot  dt 

f -8(8-8) 

f  (/  ■  r<8(d  -  8))  +  (f  ■  rM»  ~  »)) 

ot  0(7 

=  }■  8(8-6) 

_a  |7-ris(fl-fln  d  r/T^Cfl-gn 

dt  L  f  ■  8(8-8)  J  d0  L  /■  8(6-8)  J 


(A-3) 


( (f  ■  r ,8(0  -  *))  •  j-  f  ■  98(8  -«)-(/■  r&W  -9))~f-8(6-  9) 

OU  ou 


To  this  point,  the  expression  for  (rt)  is  formally  exact.  To  make  further  progress 
the  linearization  approximation  is  utilized  in  order  to  make  the  replacement 
To  linear  order,  then,  from  (A-3) 

<fi> = (~ + a0  o  (a-4) 

Similarly  one  may  show,  to  the  same  order  of  approximation, 

<Fl)  =  (^  +  ni0^)  Ar(d,t)'  (A"5) 

Equations  (A-4,  5)  and  the  corresponding  equations  for  <z,}  and  (z,>  are  used  in 
the  text  to  obtain  (7a,  b)  from  (6a,  b). 
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Collective  effects  in  the  /  =  2  steilatron  [Phys.  Rev.  Lett.  50, 507  (1983)],  a  high-current  electron 
accelerator,  are  studied.  A  thin-beam  model  is  employed,  and  only  long-wavelength,  low- 
frequency  modes  are  considered.  The  eigenvalue  problem  is  formulated  in  general  and  solved 
analytically  in  the  Smooth  approximation;  comparison  with  a  numerical  solution  of  the 
eigenvalue  problem  is  excellent.  It  is  found  that  the  dispersion  relation  in  the  smooth 
approximation  is  identical  in  form  to  that  for  the  modified  betatron,  under  a  simple  substitution. 
An  analytical  expression  for  the  transition  energy,  obtained  earlier  by  a  simple  dynamical 
argument,  is  confirmed.  The  stellarator  field  is  found  to  reduce  the  growth  rate  of  the  negative 
mass  instability. 


I.  INTRODUCTION 

The  1  =  1  steilatron,  a  betatron  with  an  added  /  =  2  stel¬ 
larator  winding,  was  originally  proposed1  in  order  to  im¬ 
prove  the  tolerance  of  a  betatron  or  modified  betatron2  to  a 
mismatch  between  the  average  beam  energy  and  the  vertical 
magnetic  field.  Such  a  tolerance  would  be  a  desirable  feature 
of  any  high-current  accelerator,  since  injectors  inevitably 
produce  a  beam  with  measurable  energy  spread  as  well  as  a 
variable  average  energy.  Though  the  original  analysis  in¬ 
cluded  space-charge  effects  in  a  simple,  non-self-consistent 
way,  the  argument  for  improved  momentum  compaction 
can  be  based  on  single  particle  orbit  considerations  alone. 

Treatment  of  space-charge  effects  in  the  steilatron  is 
complicated  by  the  lack  of  symmetry  in  the  fields.  Self-con¬ 
sistent,  monoenergetic,  KV-like3  equilibria  can  be  construct¬ 
ed,4  however,  and  these  give  explicit  relations  among  the 
beam  radii  (the  beam  is  elliptical,  in  general),  emittance, 
current,  energy,  and  externally  applied  fields;  the  shifts  in 
the  betatron  frequencies  due  to  space  charge  are  also  ob¬ 
tained  from  these  equilibrium  studies. 

In  the  present  paper  the  stability  of  a  beam  confined  in  a 
steilatron  is  considered.  This  work  was  motivated,  in  part,  by 
certain  experimental  results  on  the  /  =  2  steilatron  at  the 
University  of  California  at  Irvine.5  Though  beam  currents 
have  been  increased  in  that  experiment  by  the  addition  of  the 
1  =  2  windings,  currents  are  still  not  as  large  as  one  might 
expect  from  the  amount  of  injected  charge.  An  instability  of 
some  kind  may  be  responsible  for  the  current  loss.  An  at¬ 
tempt  to  discover  possible  current  limiting  mechanisms  in 
the  steilatron  initiated  the  work  presented  here. 

The  present  stability  calculation  makes  the  following 
approximations:  All  first-order  fields  are  treated  in  the  long- 
wavelength  (X>as  minor  radius  of  toroidal  chamber), 
low-frequency  (a<c/a)  limit;  toroidal  corrections  to  the 
first-order  fields  are  neglected.  We  note  in  particular  that  the 
long-wavelength  approximation  may  omit  the  potentially  in¬ 
teresting  case  X  ~  iflTj/m,  where  r0  is  the  major  radius  of  the 
torus  and  m  is  the  number  of  stellarator  field  periods  around 
the  device.  The  beam  is  treated  as  a  structureless,  thin, 
threadlike  charge  distribution  so  that  the  effects  of  betatron 
oscillations  and  energy  spread  on  the  equilibrium  are  ig¬ 
nored;  the  magnitude  of  these  effects  can  be  roughly  estimat¬ 
ed,  however.*  The  beam  environment,  which  has  a  critical 


effect  on  beam  stability,  is  treated  simply:  the  beam  is  as¬ 
sumed  to  be  contained  in  a  perfectly  conducting  toroidal 
chamber.  Only  the  longitudinal,  or  negative  mass  mode  is 
unstable  in  these  circumstances.  Several  comments  are 
below,  however,  regarding  the  effects  of  walls  of  arbitrary 
impedance  on  both  longitudinal  and  transverse  modes. 

Under  the  above  approximations  we  have  the  following 
findings  to  present:  ( 1 )  the  eigenvalue  problem  governing 
the  longitudinal  and  transverse  modes  of  a  beam  in  a  steila¬ 
tron;  (2)  a  dispersion  relation,  obtained  from  an  approxi¬ 
mate  solution  to  the  eigenvalue  problem,  which  is  identical 
in  form  to  that  obtained  for  the  modified  betatron,7  with  the 
redefinition  of  a  certain  symbol;  ( 3 )  confirmation  of  an  ana¬ 
lytical  expression,*  obtained  by  a  simpi*  dynamical  argu¬ 
ment,  for  the  so-called  transition  energy,  below  which  the 
longitudinal  mode  is  stable  (in  the  absence  of  dissipative 
effects);  (4)  a  quite  favorable  numerical  comparison 
between  the  solution  to  the  eigenvalue  problem  of  ( 1 )  and 
the  dispersion  relation  of  (2)  for  the  longitudinal  mode. 

This  paper  is  organized  as  follows.  Section  II  describes 
the  derivation  of  the  linearized  equations  of  motion,  assum¬ 
ing  the  applied  fields  are  given  near  the  axis  of  the  torus. 
Section  III  outlines  the  derivation  of  the  first-order  charge 
and  current  densities  in  terms  of  the  single  particle  variables 
using  a  summation  over  initial  conditions.  In  Sec.  IV  the 
approximate  solution  to  Maxwell’s  equations  is  obtained, 
and  the  total  first -order  fields  to  be  used  in  the  equations  of 
motion  are  derived.  In  Sec.  V  the  eigenvalue  problem  is  for¬ 
mulated  and  is  shown  to  reproduce  the  modified  betatron 
dispersion  relation  in  the  absence  of  the  twisted  quadrupole 
field.  An  approximate  solution  based  on  a  two-time  scale 
analysis  is  presented  and  some  comments  are  made  on  the 
form  of  the  resulting  dispersion  relation.  Finally,  in  Sec.  VI, 
a  numerical  solution  to  the  full  eigenvalue  problem  is  pre¬ 
sented  and  compared  to  results  from  the  dispersion  relation 
of  Sec.  V. 

Section  VII  summarizes  these  results  and  makes  a  few 
additional  remarks. 

II.  EQUATIONS  OF  MOTION 

We  shall  use  a  standard  (r,  6,  z)  cylindrical  coordinate 
system.  In  the  unperturbed  state  we  have  a  thin  beam  en¬ 
closed  in  a  perfectly  conducting  torus  of  major  radius  r0  and 
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minor  radius  a.  The  position  of  an  electron  which,  in  the 
absence  of  ac  excitation  would  be  at  (rQ,  9Q(t),  0)  is,  in  the 
presence  of  the  excitation,  at 

r=  (ro  +  i)?o  +  y#o  +  (D 

where  §,  y,  and*,  known  as  polarization  variables,’10,11  are 
functions  of  90  and  t  and  where  the  unit  vectors  in  ( 1 )  refer 
to  the  unperturbed  position  of  the  particle  at  time  t.  The 
applied  fields  seen  by  a  particle  at  (r0  +  r„  9,  z, )  are  taken  to 
be 

—  nB^zx/r0  +  k,3s(zx  cos md  —  rx  sin/nd),  (2a) 

o,  (2b) 

1  —  nrx/r0)  +  k,B,(rx  cos m9  +  zx  sin m9), 

(2c) 


<5J(r,  t ) 

=  £p</  +  (  -  &>o  ~  rPo  +  ^—Po  -  xrao  ^-2.)  9 

+  XP<2-  (6; 

The  perturbed  line-charge  density,  which  will  be  needed  be¬ 
low,  is  obtained  by  integrating  (5)  over  the  beam  cross  sec¬ 
tion: 


&i=a0(— i-l.ii), 

W0  r0  39  J 


(1) 


where  A0  is  the  zero-order  line-charge  density.  In  (7)  the 
first  term  is  due  to  radial  expansion  of  the  ring  while  the 
second  term  is  due  to  azimuthal  bunching  of  charge. 


where  3#,  B9Q,  and  k,3 ,  are  constants,  n  is  the  betatron 
field  index,  and  nt  is  an  integer,  the  number  of  quadrupole 
field  periods  around  the  device.  If  the  rf  fields  (as  yet  un¬ 
known)  caused  by  the  perturbation  are  denoted  by  a  super¬ 
script,  the  linearized  equations  of  motion  become 

£ +  col(l -n-rl+ peas  m90)g  -  y 

-  &0oX  +P&1  sin  m9(X 

=  -  (e/m,y0)(E<*  +  0^^),  (3a) 

y  +  o)^^  -  (*/m.^)££°,  (3b) 

X  +  <al  (n  -p  cos  md0)x  +  n80§  +pco\  sin  md <£ 

=  -  (e/m.roHE?  ~0oB™),  (3c) 

where  a  dot  denotes  3  /dt  +o)0d  /390>  o>0  =  eB^/m.YcP, 
Cl90  —  eBt0/m4y<p,  —  e  and  m,  are  the  charge  and  mass  of 
the  electron,  0O  and  y0  are  the  usual  relativistic  factors,  and 
p  39  ksoB./Bn  is  a  dimensionless  quadrupole  strength.  In 
the  absence  of  rf  fields  we  integrate  ( 3b )  once  and  substitute 
in  (3a)  to  obtain  the  usual  equations1  for  betatron  oscilla¬ 
tions.  To  complete  the  calculation  it  is  only  necessary  to 
obtain  expressions  for  the  rf  fields  in  terms  of  £,  y,  and  X- 


III.  SOURCES  FOR  MAXWELL'S  EQUATIONS 

To  express  the  first-order  charge  and  current  densities 
in  terms  of  the  polarization  variables  we  begin  with  the  for¬ 
mal  Klimontovich  expression  for  the  total  charge  density  at 
the  point  r  at  time  t 

p(r,t)=*  -ej  J3r0d3p0F(r0,  p0)<5°’(r  -  r),  (4) 

where  F  is  a  distribution  of  zero-order  positions  and  mo¬ 
menta  and  r  is  the  actual  (zero  +•  first-order)  trajectory 
f  =  r0  +  +  yd0  +  x*o-  Expanding  the  delta  function  in 

(4),  one  finds  for  the  first-order  charge  density 

8p(r,t)  =  -Xi-rpo-  —  ilpQ-xl£a.'  (5) 

r  or  r  99  dz 

where  y,  and  x  •«  evaluated  at  9  and  f,  and  p0  is  the  zero- 
order  density.  A  similar  treatment  for  the  current  density 
yields 


IV.  SOLUTION  TO  MAXWELL’S  EQUATIONS 


As  described  in  the  Introduction  we  will  confine  our 
attention  to  the  long-wavelength,  low-frequency  modes.  In 
particular  we  shall  assume 

(8) 

where  y^  is  the  relativistic  factor  corresponding  to  the 
wave-phase  velocity  and  An  is  the  wavenumber  in  the  direc¬ 
tion  parallel  to  the  beam  velocity.  This  inequality  allows  us 
to  replace  the  wave  operator  V2  +  oF/c2  by  just  Vf  and  so  to 
calculate  transverse  fields  from  statics.  The  longitudinal 
field,  required  in  (3b),  is  then  calculated  from  Far¬ 
aday's  law,  in  this  approximation. 

It  is  convenient  to  think  of  a  general  density  perturba¬ 
tion  as  a  superposition  of  two  types:  ( 1 )  the  beam  does  not 
move  transversely  but  suffers  a  change  in  density;  and  (2) 
the  beam  retains  its  zero-order  density  but  moves  transverse¬ 
ly.  Even  though  we  will  be  mainly  interested  in  the  longitudi¬ 
nal  mode,  we  include  perturbations  of  type  (2),  since  for 
certain  transverse  modes  the  quasistatic  approximation  still 
holds  and  these  are  not  difficult  to  retain  in  the  analysis.  The 
transverse  modes,  however,  will  be  stable  in  the  absence  of 
wall  resistivity  or  other  real  impedance. 

In  solving  Maxwell's  equations  for  either  type  of  pertur¬ 
bations  one  encounters  a  difficulty  in  treating  an  infinitesi¬ 
mally  thin  beam,  since  such  a  beam  produces  infinite  fields 
on  its  surface.  We  will,  consequently,  in  this  section  think  of 
the  beam  as  having  some  small,  but  finite  radius  r».  Our 
results  will  depend  only  very  weakly  on  the  value  of  r„ . 

A  perturbation  of  type  ( 1 )  is  represented  by  the  middle 
term  of  (S)  and  the  middle  term  of  the  9  component  of  (6). 
Calculation  of  the  transverse  fields  for  such  a  perturbation, 
for  perfectly  conducting  wall  boundary  conditions  is  ele¬ 
mentary.  Use  of  the  transverse  fields  in  Faradays's  law  gives 
E  a?  at  the  center  of  the  beam  as 


(9) 


where  we  have  assumed  an  e  ~  time  dependence  and  have 
retained  only  the  second  term  of  (7),  in  the  spirit  of  the 
cylindrical  approximation.  All  transverse  fields  vanish  at  the 
center  of  the  beam  for  perturbation  type  ( 1 ). 

The  dipole  perturbation,  type  (2),  produces  nonzero 
transverse  rf  electric  and  magnetic  fields  at  the  beam  center 
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but  no  contribution  to  (?)  since  transverse  mo¬ 

tion  of  the  beam  simply  compresses,  but  does  not  change,  the 
flux  between  the  beam  center  and  the  wall.  The  transverse 
fields  due  to  this  perturbation,  represented  by  the  first  and 
third  terms  of  (5)  and  the  first  and  third  terms  of  the  9 
component  of  (6),  are  straightforward  to  find,  again  using 
statics 


E,'  =  -2irpog(l  -rf/a2). 

(10a) 

=  -  torpid  -ri/a2). 

(10b) 

(10c) 

K  - 

(10d) 

The  fields  in  (10a)-(10d),  however,  are  not  the  only 
first-order  fields.  There  is  another  contribution  as  may  be 
seen  from  the  following  argument:  Any  field  component 
F(x )  may  be  written  F0(x)  +  Ft(x),  where  F0  is  the  zero- 
order  field  and  is  the  first-order  field  due  to  the  first-order 
sources  defined  at  the  observation  point  x.  In  the  equations 
of  motion  all  fields  are  evaluated  at  the  location  of  the  beam 
x0  +  jc,.  Consequently,  the  correct  first-order  fields  to  be 
used  in  (3a)-(3c)  are  sums  of  the  form 


£ 

dx 


x,  +Fl(x0), 


vhere  Ft  is  given  by  ( 10).  We  find  then  that 


(11) 


E)?  =  2(A0/a2)g,  (12a) 

E^=2(A0/a2)x.  (12b) 

while  ( 10c)  and  ( lOd)  still  hold. 


V.  THE  EIGENVALUE  PROBLEM 

Substitution  of  the  fields  in  (9),  ( 12a),  ( 12b),  ( 10c), 
and  ( lOd)  in  the  linearized  equations  (3a)-(3c)  gives 

D  2g  4-  ( 1  -  n  -  n*  -  +  p  cos  m9)g  —  Yo  Drj  —  b  Dx 

4-  fi  sin  mdx  =  0,  (13a) 

(1  ~a)4z7Tl~  -u20-  +  a0l)rj  +  Dg  =  0, 

dtf  dtf 

(13b) 

D  *x  4-  (n  —  ntb  —  n  cos  m9)x  +  b  Dg  +ps\n  m9g  —  0, 

(13c) 

where  D  —  3  /39  —  ic5,Zi  =  co/co0,  n,b  =  2vri/(y%0la2),  v 
is  Budker’s  parameter,  v  =  -  eA0 /m,c2,  a  =  2v(J  In  a / 

b  and  we  have  dropped  the  sub¬ 

script  0  on  9.  If  we  define  0  to  be  a  column  vector  with 
components  g,  g ',  rj,  if',  y,  X‘,  where  a  prime  denotes  d /dd, 
the  eigenvalue  condition  determining  57  is 

0(0  +2ir)  =  0(0).  (14) 

In  the  case^  =  0  the  condition  for  ( 14)  to  be  satisfied  is 
most  easily  obtained  by  assuming  for  0  the  dependence  e‘10, 
where  /  is  an  integer,  the  azimuthal  mode  number.  For  this 
case  (13)  reduces  to  the  evaluation  of  a  3  x  3  determinant  for 
the  eigenvalue  condition,  with  the  result 

y\  (  1  n2  -  vj  \ 

Jiff  VTT  +  W-  vj7(H2 -  v?)  - b 2n2)  ’  (I5) 


where  v2  -  21 2( v/y0)  (J  +  In  a/rk  )/tf»,  2  =  I  -  aV/ 

(Pc2),  0=a>-/,  ^-1-/1-^,  and  P,=n-n». 
Equation  (15)  agrees  with  the  result  of  Ref.  7,  where  the 
dispersion  relation  for  the  modified  betatron  is  obtained  un¬ 
der  assumptions  similar  to  those  employed  here.  The  solu¬ 
tions  to  ( 1 5 )  are  thoroughly  analyzed,  in  the  approximation 
YPh  =  Yo>  in  Ref.  7,  where  resistive  wall  effects  are  also  con¬ 
sidered;  the  effects  of  a  more  general  wall  impedance  are 
examined  in  Ref.  6. 

Let  us  now  return  to  the  eigenvalue  problem  for  the  case 
ft  0.  For  long-wavelength  modes,  one  expects  each  o(g,  77, 
and  x  to  be  composed  of  a  main,  slowly  varying  part  and  a 
much  smaller,  rapidly  oscillating  piece,  this  latter  pan  due  to 
the  rapidly  oscillating  quadrupole  focusing.  This  suggests 
that  we  write  g  —  g,  4-  g  f ,  and  similarly  for  rj  and  x<  where 
the  subscripts  denote  slow  and  fast  variation.  It  should  be 
said  here  that  we  anticipate  that  |m|>l  in  typical  applica¬ 
tions;  the  UC  Irvine  experiment,  for  example,  uses  m  =  12. 
Substituting  for  g,  r),  and  x  in  ( 13a)— ( 13c)  and  performing 
an  average  over  the  fast  scale,  2ir/m,  gives  the  equations 

D2g,  +  (1  -/i-n* -li)g,  +n(ca*m9gf) 

-Yo  Di),  -bDx,  =0,  (16a) 


—  2i  0-^—  Vi  —  <»2(1  +a0o)rh  A- Dg,  =0,  (16b) 

39 

D2X,  +  (n-na )x.  - /i<coe mdxf  >  4-h Dg, 

+ p(sin  m9g *0,  (16c) 

where  the  angular  brackets,  (— >,  denote  an  average  over  the 
fast  scale,  and  we  have  assumed  (g,)  =  g,,  (gf) 
=  0  =  (g  / )  =  (g  7  >,  etc.  If  we  now  subtract  ( 16a)-(  16c) 
from  the  unaveraged  equations  ( 13a)— ( 13c),  use  \gf 1  «|£,  |, 
keep  all  derivatives  acting  on  fast  quantities,  and  ignore 
quantities  of  the  form  cos  mOg f  —  (cos  mOgf  >,  we  arrive 
at  equations  for  the  fast  variables: 


d*  d 

—  g/~b  —  X/=  -/ucos  mdg,  —  M  sin  m&x„ 


39 


1L 

39 2  Xr 


—  p  sin  m9g,  +  p  cos  m9x„ 


which  are  easily  solved  to  give 


(17b) 


gf  —  [  p/m(m  +  b)\(coim9gt  4- sin m9x,),  (18a) 
Xf  —  [  fi/m(m  4-  b) )(sin  m9g,  —  cos  m0y, ).  (18b) 

Finally,  substitution  of  ( 18a)  and  ( 18b)  into  ( 16a)-(  16c), 
using  (cos2  m9 )  =  (sin2  m9)  gives  the  equations  for 
the  slow  components: 

D2g,  4-  0  —  n  -n,t  +  [p2/m(m  +  b)]  -  Ya}gt 

-  Yo  Dy,  -  b  Dx,  =  0,  (19a) 

(1  -a)  w*71,  ~2iW~hv>  -s?2(1 +a^°)77' 


4-  Dg,  =  0, 

D  2x,  +  0*  -  n*  4-  [  MJ/m(m  4-  b) )}*,  4-  b  Dg, 


(19b) 

=  0. 
(19c) 
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Equations  ( 19a)— { 19c)  are  differential  equations  with  con¬ 
stant  coefficients  and  are  in  fact  precisely  of  the  form  of 
( 1 3a)— ( 13c)  if  is  set  to  zero  in  ( 13a)— ( 13c)  andifn,*  is 
replaced  by 

n,„  — ■ **  ~  +  6)].  (20) 

It  is  seen  that  the  focusing  effect  of  the  twisted  quadru¬ 
ple  has  been  replaced  by  a  smoothed,  effective  focusing 
term.  It  follows,  consequently,  from  ( 19a)— ( 19c)  that  the 
dispersion  relation  for  the  stellatron  for  long-wavelength 
modes  is  just  that  for  the  modified  betatron,  Eq.  ( 15),  under 
the  replacement  (20).  This  is  the  basic  result  of  this  paper. 

A  further  result  may  be  obtained  by  examination  of 
( 1 5 ) .  If  we  approximate  and  assume  |  Cl2 1  <  v? ,  vj  for 

the  longitudinal  mode,  then  solving  the  resulting  quadratic 
equation  gives  the  result  that  the  longitudinal  mode  is  stable 
below  a  transition  energy  given  by 

yu  —  v,  =  {1  -  n  -n*  +  [fi2/m(m  +  b)]}l/2,  (21) 

which  agrees  with  an  earlier  estimate.8 

For  the  special  case4  m  ■+■  b  =  0,  the  analysis  leading  to 

(20)  and  (21 )  no  longer  applies.  In  this  case  one  may  show, 
by  reconsideration  of  ( 17a)  and  ( 17b)  ff,  that  the  stellatron 
dispersion  relation  is  again  (IS)  but  now  under  the  replace¬ 
ment  bft-+bCl  —  y}/bl.  A  transition  energy,  as  defined  in 

(21 ) ,  no  longer  exists;  one  is  always  above  transition. 

VI.  NUMERICAL  ANALYSIS 

In  order  to  test  the  analysis  leading  to  the  dispersion 
relation  described  in  the  previous  section,  a  numerical  solu¬ 
tion  to  the  eigenvalue  problem  for  ( 13a)— ( 13c)  has  been 
carried  cut.  The  solution  is  based  on  the  following  construc¬ 
tion.  We  define  a  6  x  6  matrix  'P  ( 0 )  which  has,  as  columns, 
solution  vectors  satisfying  ( 13a)-(  13c).  Furthermore,  we 
define  ♦  to  satisfy  the  initial  condition  't'(O)  =  I,  the  identi¬ 
ty  matrix.  Since  the  coefficients  in  ( l3a)-(  1 3c)  are  periodic 
with  period  2v/m  ==  0m  it  follows  that  V(0  +  0m  )  must 
also  be  a  solution  to  ( 1 3 )  and  so  the  columns  of  ¥  ( ff  +  6m ) 
must  each  be  expressible  as  a  linear  combination  of  the  co¬ 
lumns  of  'V(d).  Considering  the  initial  condition  on  we 


have  that 

+dm)  =  ¥(ff)*(0,  (22) 

from  which  it  follows,  by  induction 

V(0  +Jc0m)  =  *(0)'t>k(0m)  (23) 

for  any  non-negative  integer  k. 

Suppose  we  find  the  matrix  U  which  diagonalizes 

n0m)U=U\,  (24) 

where  A  is  a  diagonal  matrix,  containing  the  eigenvalues  of 
'V(0m).  The  columns  of  ¥(#)(/ are  then  the  Floquet  solu¬ 
tions,  satisfying 

V{0  +  0m)U=V(0)U\.  (25) 

It  follows  from  (25),  then,  that 

*(0  +  2n)U=V{0)U\m.  (26) 

The  eigenvalue  problem  for  ( l3a)-(  1 3c)  may  now  be  stated 


succinctly  as  follows:  One  must  find  S  such  that  (at  least) 
one  of  the  eigenvalues  of  V{9m )  is  one  of  the  m  roots  of 
unity. 

A  fortran  program  has  been  written  to  solve  this 
problem.  Starting  with  an  initial  guess  for  37,  the  matrix 
^  ( 8m )  is  evaluated  and  its  eigenvalues  computed;  an  adjust¬ 
ment  to  2J  is  then  made  using  a  simple  secant  method.  Con¬ 
vergence  typically  takes  three  or  four  iterations. 

For  comparison  purposes  the  smooth  approximation 
dispersion  relation  ( 15),  (20)  was  also  solved  numerically 
using  the  same  first  guess.  This  was  a  matter  of  finding  the 
appropriate  root  of  a  sixth-order  polynomial  ( the  approxi¬ 
mations  employed  to  obtain  (21)  were  not  made  here  ]  and  is 
very  fast  computationally.  In  fact  the  data  for  the  plots 
shown  below,  of  the  longitudinal  mode  growth  rate  versus 
various  parameters,  originally  took  —40  min  each  to  gener¬ 
ate  on  a  VAX  1 1-780  by  the  first  method.  Each  took  less  than 
3  sec  using  the  smooth  approximation  dispersion  relation. 
For  all  cases  tried,  the  two  methods  agree  to  at  least  three 
significant  figures  and  give  plots  indistinguishable  on  the 
scales  used  below. 

Three  plots  are  shown  here.  All  were  made  using  values 
typical  of  the  device  at  UC  Irvine,  i.e.,  a/r0  =  0.1,  rb/a 
=  0.2,  n  =  0.5,  m  =  12,  and  b  =■  40.  In  the  first  plot.  Fig.  1, 
Im  (37)  is  plotted  versus  current  for  a  beam  energy  of  2  MeV 
and  fi  =>  50.  The  dependence  follows  quite  closely  to  the 
classical  square  root  law.  In  Fig.  2,  the  growth  rate  is  plotted 
versus  energy  for  a  current  of  200  A  and  n  =  50.  The  value 
of  the  transition  energy  predicted  by  (21)  is  =0.52  MeV  and 
is  seen  to  agree  quite  well  with  the  numerical  value.  Finally 
in  Fig.  3  the  growth  rate  is  plotted  veisus  p  for  a  beam 
energy  of  2  MeV  and  a  beam  current  of  200  A.  The  stabiliz¬ 
ing  effect  of  the  quadrupole  field  can  be  seen  quite  clearly. 
This  is  to  be  expected  of  any  field  which  inhibits  transverse 
motion.  In  all  cases  shown  Re(aJ)  —  1. 

VII.  DISCUSSION  AND  SUMMARY 

The  features  and  various  limiting  cases  of  the  dispersion 
relation  for  the  modified  betatron  have  been  extensively  dis- 


FTG.  1.  Normalized  growth  rate  for  the  longitudinal  mode  versus  current. 
Beam  energy  is  1  MeV.  Here  n  -  50.  (Other  perimeters  given  m  text.) 
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FIG.  2.  Normalized  growth  rate  for  the  longitudinal  mode  venus  energy. 
Beam  current  is  200  A.  Here  n  =  30.  The  theoretical  value  for  the  transi¬ 
tion  energy  is  0.32  MeV. 


FIG.  3.  Normalized  growth  rate  for  the  longitudinal  mode  venus  quadru- 
pole  strength.  Beam  energy  is  2  MeV;  beam  current  is  200  A. 


cussed  in  Refs.  6  and  7.  Toroidal  corrections  and  a  careful 
comparison  with  the  results  of  particle  simulations  have 
been  carried  out  by  Godfrey  and  Hughes. 12  We  have  shown 
that,  under  the  substitution  ( 20 ) ,  virtually  all  of  the  analyses 
already  carried  out  on  the  long-wavelength  stability  of  the 
modified  betatron  may  be  carried  over  to  the  stellatron.  In 
the  long-wavelength,  quasistatic  limit  the  net  effect  of  the 
addition  of  stellarator  windings  to  a  modified  betatron,  we 
have  found,  is  to  add  a  smooth,  transverse  focusing  term 
which  can  be  large  compared  to  the  wcak-foc  using  or  space- 
charge  defocusing  terms.  Stellarator  fields  enhance  stability 
in  this  sense.  Other  proposed  stellaratorlike  focusing 
schemes  may  be  treated  in  a  similar  manner.  There  are  two 
resulting,  related  effects  on  the  dispersion  relation.  One  is 
the  introduction  of  a  finite  transition  energy,  given  by  (21 ), 
below  which  the  longitudinal  mode  growth  vanishes  alto¬ 
gether  and  the  second  is  a  simple  shift  (increase)  of  the 
transverse  frequencies,  v,  and  vx .  Of  course,  below  the  tran¬ 
sition  energy,  dissipative  effects  (a  resistive  wall,  for  exam¬ 
ple)  will  cause  both  longitudinal  and  transverse  modes  to 
grow,7,13-14  but  at  a  growth  rate  which  is  smaller  the  larger 
the  quadrupole  focusing.  Such  modes  would  then  presuma¬ 
bly  be  more  easily  stabilized  by  effects  neglected  here,  espe 
daily  azimuthal  velocity  spread. 

The  smooth  approximation  of  Sec.  V,  above,  will  break 
down  for  l>m/ 2.  Also,  the  solution  to  Maxwell’s  equations 
obtained  in  Sec.  IV  becomes  unreliable  for  l^r^/a,  which  is 
typically  of  the  order  of  m.  Consequently,  short- wavelength 
modes  must  be  studied  by  entirely  different  methods  from 
those  employed  here.  In  fact,  there  is  some  evidence13  from 


particle  simulations  for  strong  interactions  of  the  beam  with 
certain  electromagnetic  modes  of  the  toroidal  cavity,  under 
certain  conditions.  These  interactions  are  an  important  topic 
for  further  study. 
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It  is  generally  difficult  to  diagnose  the  energy  of  an  electron  beam  propagating  in  a  focusing 
channel.  Here  we  suggest  observation  of  the  Doppler-shifted  radiation  due  to  the  electrons’ 
betatron  motion  in  the  channel  as  a  sensitive,  nonperturbing  diagnostic.  Sample  spectra  have 
been  calculated  for  a  specific  type  of  channel  used  in  a  cyclic  electron  accelerator  and  their 
features  are  discussed. 


I.  INTRODUCTION 

For  some  recently  proposed  high-current  electron  ac¬ 
celerators  it  would  be  desirable  to  have  a  passive  diagnostic 
of  the  various  beam  parameters.  Beam  current  measure¬ 
ments  may  be  made  by  a  variety  of  means,  but  passive  mea¬ 
surements  of  beam  density  and  energy  are  less  easy.  One 
might  consider  observation  of  ordinary  synchrotron  radi-. 
ation  in  cyclic  devices  as  an  energy  diagnostic,  but  this  is 
effective  only  at  moderately  high  energies;  at  the  lower  ener¬ 
gies  attained  so  far  in  these  devices  the  metallic  chamber  acts 
very  effectively  as  a  shield  for  this  long  wavelength  radiation. 
Higher-frequency  radiation,  significantly  above  cutoff  of  the 
beam  tube,  is  emitted  by  the  electron.;  as  they  execute  beta¬ 
tron  oscillations.  This  radiation  is  Doppler  shifted  in  the 
forward  direction  by  a  factor  proportional  to  the  square  of 
the  energy  (for  0—  1 );  hence,  it  is  a  good  candidate  for  an 
energy  diagnostic. 

This  paper  presents  an  analysis  of  the  radiation  from 
betatron  oscillations  for  a  particular  device,  the  stellatron1 
( a  strong-focused  betatron ) ,  one  of  which  has  been  built  and 
operated  at  the  University  of  California  at  Irvine.2  We  will, 
however,  make  a  number  of  general  comments  applicable  to 
many  similar  machines. 

This  paper  is  organized  as  follows:  In  Sec.  II,  we  de¬ 
scribe  the  magnetic  field  configuration  of  the  stellatron, 
write  down  the  paraxial  equations  of  motion,  and  write 
down  their  solution,  giving  the  expressions  for  the  four  beta¬ 
tron  frequencies;  this  section  is  brief  since  the  results  are 
given  elsewhere.' 3  It  is  mainly  included  to  introduce  some 
notation.  Next,  we  examine  the  standard  expression4  for  ra¬ 
diated  intensity  for  prescribed  particle  motion,  and  make 
some  comments  on  what  one  expects  to  observe.  We  then 
give  some  specific  examples  of  spectra  and  discuss  their  fea¬ 
tures,  including  the  effect  of  “geometric  line  broadening,”  an 
inhomogeneous  broadening  mechanism  due  to  orbit  curva¬ 
ture.  The  final  section  draws  conclusions  and  makes  a  few 
additional  remarks. 

II.  PARTICLE  ORBITS  IN  THE  STELLATRON 

The  stellatron  is  a  betatron  with  superposed  strong  fo¬ 
cusing  ir.  the  form  of  an  /  =  2  stellarator  (twisted  quadru¬ 
ple)  field;  we  take  the  accelerator  to  be  in  the  form  of  a 
torus  of  major  radius  rQ  and  a  minor  radius  a.  The  magnetic 


fields  near  the  null  line  (taken  to  be  a  circle  of  radius  r0)  of 
the  quadrupole  are,  in  cylindrical  coordinates  (r,dj), 

Br=x-  nB^y  +  kwB,r0(y  cos  kws-x  sin  k^s),  ( la) 
B9=xBg0,  (lb) 

B,  =*11,0  (1  -nx)  +  kwB,r0(x  cos  kws  +  y  sin  k„s), 

(lc) 


where  BM,  B„  and  Be0  are  constants,  n  is  the  betatron  field 
index,  kw  is  2jt/( period  of  quadrupole  twist),  and  x  and  y 
are  normalized  radial  and  vertical  displacements,  respective¬ 
ly  :  x  =  ( r  —  r0 )  /r^  y  —  z/r q.  If  an  electron  has  energy  mc2y0 
then  its  first-order  equations  of  motion  are 

x"  +  [(1  -  n)/r%  +k\ cosA^r]* 

=  - yk 2  sin  kws  +  k„y',  ( 2a) 

y"  +  [n/rl  -  k]  cosAr„,j].y  =  ~xk\  sin kws  -  kbx',(2b) 

where  k  2q  =  kwB,/r0Bt0,k„  =  B90/r0Bz0,  and  a  prime  de¬ 
notes  d  /ds.  Note  that  in  the  cylindrical  limit  r0  —  oo ,  kq ,  kw , 
and  kb  all  remain  fixed.  Only  the  focusing  terms  due  to  cur¬ 
vature  in  the  vertical  field,  represented  by  the  1  —  n  and  n 
terms  in  (2a)  and  (2b),  vanish. 

The  Floquet  solutions  of  (2a)  and  (2b)  generally  con¬ 
tain  an  infinite  number  of  terms  but  in  the  special  case  n  =  1, 
a  symmetry  exists  in  (2a)  and  (2b)  and  one  may  find  the 
general  solution  as  a  sum  of  four  oscillating  terms: 

x  +  iy  =  e'kjs/\A  +  a +A  *  e  "  '*•  * 


where 


+  A_e,k-’  +  cr_A*_e-lk-s), 


1  ,  k\  (kw  +k„)2 


k2  =  _ _ l  _L  -i- 

*  +  ■  + 


fa+*L) 

L\2rg  4  J 


(kw  +  kb)2  +  k* 


and 


(3) 


(3a) 


at  +  kb)/D2  -  \/2/i  ~  k  ;/4], 

In  (3)  the  complex  constants  A  ±  depend  on  the  four  initial 
conditions.  The  solution  (3)  is  stable  if 
|Jk2  +  k^kg  —  lAo|  >  1 2k  J|,  which  sets  the  limit  on  the 
quadrupole  strength. 

In  the  experimentally  interesting  limit  of  large  kh  ( large 
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toroidal  field)  the  betatron  frequencies  take  the  forms 


+  k+)0ocs;(kb  +  km)0(f,  (4a) 

a>2s(\kw  -k+)0oc~  -kb0<f,  (4b) 

co,  =  ([kw  +  k_)0oc^kw0(fi,  (4c) 

a>4mHkm-k_)fas:fo/l2iikb).  (4d) 


Note  that  kb0cf  is  just  the  electron  cyclotron  frequency  in 
the  toroidal  field,  eB^/my^.  Of  the  four  modes,  only  the 
first  two  typically  emit  radiation  at  frequencies  significantly 
above  the  beam  tube  cutoff,  ~c/a. 

In  the  opposite  limit  kb  —0  (no  longitudinal  magnetic 
field),  we  obtain  the  forms  of  the  betatron  frequencies  in  the 
cylindrical  limit  as3 

[kw/2±(iki±kl)U2]0dc.  (5) 

Of  course  since  the  betatron  frequencies  depend  on 
many  parameters  (magnetic  fields,  field  gradients,  etc.)  one 
may  consider  a  measurement  of  a  betatron  frequency  as  a 
measurement  of  any  of  these  parameters;  most  of  these,  how¬ 
ever,  are  most  easily  and  directly  measured  in  other  ways. 
An  exception  is  a  measurement  of  beam  density.  Although 
we  have  not  included  beam  self-field  effects  in  our  evaluation 
of  the  betatron  frequencies  one  may  do  so1;  if  this  is  done, 
however,  one  finds  that  only  the  lowest  betatron  frequency  is 
at  all  sensitive  to  beam  density  for  typical  values  of  other 
parameters.  This  lowest-frequency  mode  is  typically  far  be¬ 
low  cutoff  of  the  beam  tube  and  so  is  most  efficiently  ob¬ 
served  by  simple  electrostatic  probes. 

III.  RADIATION  SPECTRA 

The  energy  radiated  per  unit  solid  angle  per  unit  fre¬ 
quency  by  a  single  electron  moving  along  a  specified  trajec¬ 
tory  in  free  space  is  given  by4 

dI(ail  =  -L  ( — V  I  fdrfix(nX3)c'"l'"A‘r<')/c,r. 

dO  c  \2ir )  \J  i 

(6) 

where  h  is  a  unit  vector  pointing  from  the  electron  to  a  dis¬ 
tant  observer  and  the  integral  extends  over  times  during 
which  3  is  changing.  Equation  ( 16)  is  applicable  to  the  radi¬ 
ation  from  electron  orbits  in  the  stellatron  for  frequencies  far 
above  the  beam  pipe  cutoff;  it  also  assumes  that  collective 
effects — electron  bunching  on  a  scale  of  a  radiation  wave¬ 
length — may  be  neglected.  In  order  to  calculate  the  spec¬ 
trum  of  radiation  emitted  by  an  electron  circulating  in  the 
stellatron,  it  is  possible,  if  rather  tedious,  to  substitute  (3) 
into  (6),  expand  the  exponentials  in  terms  of  Bessel  func¬ 
tions,  and  perform  the  integral  over  time.  The  result  one  sees 
involves  spectral  peaks  at  the  frequencies 

ntfOQ  4-  n,aj2  +  nz<o2  +  n}co}  +  n4o41  (7) 

where  <u0  is  the  basic  “cyclotron"  frequency  =  0O c/r0  and 
rt0,...,n4  are  integers.  The  amplitude  of  the  radiation  for  any 
givt.i  set  of  n' s  depends  on  the  energy  of  the  electron  and  on 
its  initial  conditions. 

Rather  than  write  out  the  explicit  but  nontransparent 
mathematical  results  for  the  integral  in  (6),  let  us  discuss 
what  one  expects  to  see.  Consider,  first  of  all,  a  particle  mov¬ 
ing  not  in  a  torus  but  along  the  axis  of  a  cylinder  with  veloc¬ 


ity  0ffi,  executing  a  betatron  oscillation  at  a  frequency  aB 
about  the  axis.  An  observer  stationed  on  the  axis,  with  the 
electron  moving  toward  him,  will  see  radiation  at  the 
Doppler-shifted  value  o)e/{  1  -  0O).  If  the  observer  is  locat¬ 
ed  at  some  angle  6  off  the  axis  of  the  cylinder,  he  sees  a 
smaller  Doppler  shift  co0/(  1  —  0O  cos  6).  Therefore,  a  de¬ 
tector  centered  on  the  axis  with  a  small  finite  angular  accep¬ 
tance  A0  will  see  a  spread  in  radiation  frequencies  of  the 
order  of 

0o(M)2a>g/{2{l  -0o)[l-0o  +0o(m1/2  ]}. 

In  addition,  as  the  observer  moves  off  the  axis  he  begins  to 
see  energy  radiated  in  harmonics  of  the  Doppler-shifted  fun¬ 
damental. 

This  “geometric  line  broadening”  is  illustrated  in  Fig.  1 
and  2  for  an  electron  moving  along  the  axis  of  a  cylinder  in 
which  focusing  is  provided  by  both  a  twisted  quadrupole  and 
longitudinal  magnetic  fields.  Figure  1  shows  the  spectrum  as 
seen  by  an  observer  located  exactly  on  the  cylinder  axis  ( i.e., 
6  =  =  0).  Three  peaks  are  seen;  a  fourth  is  actually  pres¬ 
ent  but  at  indistinguishable  amplitude  on  this  scale.  The  four 
peaks  precisely  correspond  to  the  four  Doppler-shifted  beta¬ 
tron  frequencies  co„...^o4  for  this  case.  The  finite  line  widths 
are  solely  due  to  the  fact  that  we  numerically  integrate  (6) 
only  over  a  finite  time;  the  widths  are  seen  to  shrink  as  we 
increase  the  time  of  integration.  For  a  real  cylinder  of  finite 
length,  of  course,  the  radiated  linewidth  is  nonzero.  No  har¬ 
monics  of  the  fundamentals  are  present. 

Figure  2  shows  the  result  of  an  integration  of  ( 6 )  over  an 
angular  range  yO  from  0  to  0.6  (i.e.,  A0  ==  0.6 /y  =  0.06). 
Line  broadening  and  the  appearance  of  energy  radiated  in 
harmonics  is  evident  in  this  case.  The  fundamental  peaks, 
however,  are  still  identifiable. 

The  situation  in  a  torus  is  such  that  it  automatically 
presents  any  observer  with  a  range  of  Doppler-shifted  fre¬ 
quencies,  as  well  as  with  many  harmonics  of  the  betatron  and 
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FIG.  1.  Radiation  emitted  in  forward  direction  by  an  electron  moving  along 
a  cylindrical  focusing  channel,  y  =  10 ,  k„  =  0.  t  cm  ~  B,  =  34  G/  cm 
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FIG.  2.  Integrated  radiation  spectrum  over  an  angular  range  yd  =  0—0.6. 

FIG.  4.  Radiation  emitted  by  electron  moving  in  a  torus  with  r0  =  4 1  x  10* 
cm. 

cyclotron  frequencies.  A  distant  observer  “sees”  an  electron 

moving  toward  and  away  from  him  in  a  periodic  fashion,  radii.  In  all  cases,  the  integration  time  interval  is  held  fixed, 

therefore  one  expects  radiation  in  the  range  cog  /( 1  -  0O)  to  A11  cases  share  the  values  Beo  =  8.6  kG,  ku.  Bs  =  247  G/ 

coe/(  1  +/?<,)•  The  intensity,  though,  at  the  Doppler  up-  cm>  kw  =  0.293  cm"1,  y0  =  5.23,  which  are" typical  of  the 

shifted  value  is  much  larger  than  at  the  down-shifted  value  uc  Irvine  device.  The  values  for  the  four  Doppler  up-shifted 

since  relativistic  effects  tend  to  cause  the  most  radiation  to  be  betatron  frequencies,  for  r0  =  41  cm,  are 

emitted  in  the  forward  direction.  As  a  practical  matter,  when  a . o>«  =  2.03  X  1012,  —  1.57xl0u,  4.62x10",  and 

trying  to  observe  radiation  from  high-frequency  betatron  os-  3.91x10’  rad/s. 

cillations,  one  may  expect  interference  from  very  high  har-  Figure  3  shows  the  result  for  r0  =  41  X  106  cm;  as  ex- 

monics  of  the  cyclotron  frequency  a)0.  This  may  be  mini-  pected,  the  spectrum  is  quite  clean,  consisting  of  narrow 

mized  by  limiting  observation  to  radiation  polarized  normal  ]ines  at  _  ^  and  Bj,  (The  other  two  peaks  are  extremely 

to  the  plane  of  the  torus.  weak  for  the  particle  initial  conditions  we  chose. ) 

In  Figs.  3-5,  we  show  a  series  of  spectra  for  an  electron  When  the  radius  is  reduced  to  41  x  102  cm,  the  spectrum 

moving  in  the  fields  of  a  stellatron  for  different  torus  major  jn  pjg  4  results.  The  expected  line  broadening  is  beginning 


FREQUENCY  (RAO/SEC) 

FIG.  3.  Radiation  emitted  by  electron  moving  in  a  torus  with  r0  =  41 X  10* 

cm.  FIG.  5.  Radiation  emitted  by  electron  moving  in  a  torus  with  r0  =  41  cm. 
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FIG.  6.  As  in  Fig.  5,  but  with  finer  frequency  resolution. 


to  become  evident  at  this  radius. 

Finally,  for  r0  =  41  cm,  the  spectrum  is  dramatically 
more  complicated  as  illustrated  in  Fig.  5.  The  lines  have 
become  quite  broad,  although  the  dominant  one  at  —  co2  is 
still  identifiable,  and  significant  amounts  of  energy  are  now 
being  observed  at  harmonics  and  beats  of  the  betatron  fre¬ 
quencies.  A  calculation  done  at  finer  frequency  resolution 
( Fig.  6 )  shows  the  complicated  structure  of  the  spectrum  for 
this  case. 


IV.  DISCUSSION  AND  CONCLUSIONS 

In  a  straight  section,  we  have  found  that  radiation  emit¬ 
ted  by  the  betatron  motion  of  electrons  may  provide  a  useful 
energy  diagnostic  by  measurement  of  the  forward  Doppler 
shift.  In  a  torus,  however,  geometrical  effects  greatly  compli¬ 
cate  the  form  of  the  spectrum  and  interpretation  of  the  spec¬ 
tral  peaks.  The  widths  of  the  peaks  will  depend  both  on  geo¬ 
metric  effects  and  on  the  actual  energy  spread  in  the  beam. 
Still,  it  may  be  possible  to  identify  the  various  peaks  in  the 
toroidal  case  with  careful  study.  We  have  also  considered  the 
possibility  that  radiation  from  betatron  oscillations  might 
provide  a  useful  diagnostic  for  other  quantities.  Though  this 
may  be  the  case  for  various  field  gradients  and  field 
strengths,  it  does  not  seem  to  be  useful  as  a  beam  density 
measurement,  since  only  the  betatron  frequencies  below  the 
beam  tube  cutoff  are  at  all  sensitive  to  the  value  of  the  den¬ 
sity. 
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Beam  transport  systems  in  which  the  two  transverse  degrees  of  freedom  of  particle  motion  (x  and  y) 
are  coupled  are  often  used  to  transport  high-current,  low-energy  beams.  Consideration  of  the 
matching  problem  for  such  systems  requires  a  formalism  for  the  description  of  the  evolution  of  the 
beam  envelope,  including  the  effect  of  space  charge.  By  defining  a  4  x  4  beam  covariance  matrix,  it  is 
possible  to  give  a  simple  prescription  for  following  the  rms  beam  ellipse  through  an  arbitrary  linear 
transport  system,  allowing  the  length  of  the  axes  of  the  ellipse  and  its  angle  of  tilt  to  vary;  the 
prescription  trivially  reproduces  the  usual  K-V  envelope  equations  in  the  decoupled  case.  The 
matching  condition  for  the  general  coupled  system  is  given,  and  the  stability  problem  for  the  marched 
solution  in  the  presence  of  space  charge  is  formulated  in  general  and  discussed  in  a  particular 
example. 


1.  INTRODUCTION 

The  familiar  K-V  envelope  equations1  are  applicable  to  beams  confined  by 
focusing  systems  that  do  not  couple  the  two  transverse  degrees  of  freedom  of 
particle  motion,  x  and  y.  The  beam  ellipse  remains  upright  in  these  systems,  and 
the  axes  merely  vary  in  length  along  the  direction  of  beam  propagation.  Certain 
focusing  systems  used  for  the  transport  of  high-current,  low-energy  beams, 
however,  do  couple  the  x  and  y  motion.  In  some  simple  cases  with  azimuthal 
symmetry,  e.g.  a  solenoid,  a  transformation  to  a  rotating  frame  decouples  the  two 
transverse  degrees  of  freedom,  but  in  the  general  case  the  problem  is  most 
directly  treated  in  the  coupled,  “laboratory-frame”  variables.  Systems  to  which 
the  following  analysis  may  be  applied  include  the  modified  betatron2  (for  all 
stable  values  of  field  index),  the  bumpy  torus,3  the  reversing  solenoid  lens,4  and 
the  stellarator-focused  betatron,  or  stellatron.5  For  this  last  case,  matched  beam 
solutions  have  been  constructed,6,7  but  the  stability  of  these  solutions  in  the 
presence  of  space  charge  remains  to  be  studied.  The  formalism  outlined  here 
allows  such  an  analysis  to  be  carried  th  ;.tg*i  this  is  done,  as  an  example,  below. 

The  stability  of  K-V  distributions  in  qu  T.  pole  and  solenoid  systems  has  been 
analyzed  by  Hofmann,  Laslett,  Smith,  a  .  i  Haber8  and  by  Struckmeier  and 
Reiser.9  Hofmann,  ft  al.  carried  out  a  thorough  analysis  of  the  electrostatic 
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discussed  here  in  Section  4.]  No  assumption  is  made  here,  of  course,  about  the 
exact  form  of  the  distribution  of  particle  initial  conditions.  Using  Eq.  (1)  it 
follows  that  the  matrix  Z  obeys  the  equation 

Z'(i)  =  MZ  +  (MZ)r,  (3) 

where  T  denotes  transpose. 

The  matrix  M  generally  consists  of  two  parts,  one  due  to  applied,  extenal  fields 
and  another  due  to  self,  space-charge  fields.  For  magnetic  focusing  with 
B  =  (Bt,By,B!), 


(4) 


where  Bab^  -q(dBa/db)/pc,  Bz  =  -qBJpc,  q  is  the  particle  charge,  p  is  its 
momentum,  c  is  the  speed  of  light,  and  cgs  units  are  used. 

The  space-charge  contribution  ot  M(s)  will  depend  on  the  spatial  elements  of 
Z,  viz.,  Z ,ia  Z„  *»  ayy,  and  Z13  =  Z31  m  axy,  that  is,  on  the  size  and 
orientation  of  the  beam  ellipse.  The  space-charge  forces  are  linear,  of  course, 
only  when  the  ellipse  is  uniformly  populated,  as  for  the  K-V  distribution,  which 
here  takes  the  form 

\W\ia 

f{vTWv)  =  '—L -  6[vr(r0)VF(r0)t;(Jo)  -  1],  (5) 

where  W  is  a  positive  definite,  real  symmetric  matrix.  If  we  define  V(s  1  s0)  as  that 
4x4  matrix  solution  to  the  single-particle  equations  of  motion  satisfying 
^0*0 1  Jo)  =  1,  the  identity  matrix,  then 

W(r)  =  V~lT(s  |s„)H'(s0)K-'(s|s„)  (6) 

and10 

Z(i)=4-W-'(i)=  V(s  1 5(,)Z(s„)y T(S  I  5o).  (7) 


Me“(s)  = 


0 

B, 


1 

0 


0 

Bv, 


yx  v  **yy 
0  0  0 
L-B„  Bx  -Bxy 


0  ' 
-B, 
1 
0 


The  space-charge  part  of  M  is  evaluated  in  terms  of  axx,  ayy,  and  axy  for  the 
K-V  distribution  in  the  Appendix.  The  result  may  be  expressed  as 


MK{s) 


(Py)3 


’  0 

0 
•  <?*y 


0  0  0' 
0  qxy  0 
0  0  0 
0  <7yy  0. 


(8) 


where  /  =  beam  current  in  s  direction /(mc3/q),  m  is  the  particle  mass,  /8  and  y 
are  the  usual  relativistic  factors,  q^-Sy/D,  qyy  =  Sx/D,  qxy  =  -oxy/D,  D  = 
Sn(Sx  +  S,),  Sr  =  Sn  +  a„,  Sy  =  So  +  Oyy,  and  =  (oxxoyy  -  o^v),/2. 

With  M{s)  given  by  the  sum  of  Eqs.  (4)  and  (8),  Eq.  (3)  describes  the  evolution 
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of  the  x-y  coupled  K-V  beam  envelope  in  the  presence  of  space  charge.  In  the 
decoupled  system,  with  Af**1  in  2  x  2  block  diagonal  form  (Bxx  =  Byy  =0  =  B.) 
and  aiy  =  0,  some  algebra  shows  that  Eq.  (3)  reproduces  the  usual  K-V  envelope 
equation^,  for  the  radii  rx  =  2aj£  and  ry  =  2 o‘y .  For  axy  #0,  the  K-V  beam  radii 
(rt)  and  the  angle  (a)  that  the  major  beam  axis  makes  with  the  x-axis  are  given 
by 

\r\  =  +  ayy  ±  a0  (9a) 

sin2ar  =  2axy/o0  (9b) 

cos  2 a  =  (a^  -  oyy)/ a0,  (9c) 

where  a0  =  [(oxx  -  ayy)2  +  4o2xy\in. 

In  the  decoupled  case,  the  quantities  oaax  x.  -  afx.  ■  ex  m,,  and  oyyay  y  - 
ofy‘  -  £ynm  are  conserved,  in  the  absence  of  acceleration.  When  coupling  is 
present  it  follows  from  Eq.  (7)  or,  more  generally,  from  Eq.  (3)  that  the 
determinant  of  S  is  conserved.  In  the  decoupled  case,  |S|  =  ex  rm,£y  n,  which 
suggests  that  the  quantity 

1/4  (10) 

may  be  a  useful  definition  of  beam  quality  in  the  coupled  case;  such  a  quantity  is 
easily  calculated  by  particle-tracking  codes. 


3.  MATCHED  SOLUTIONS  AND  THEIR  STABILITY 

We  consider  the  simplest  definition  of  a  matched  solution  to  Eq.  (3)  in  a  periodic 
focusing  system  with  period  L;  for  all  s 

Z(*  +  L)  =  Z(s).  (11) 

A  slightly  more  general  definition,  which  specifies  that  the  beam  ellipse  rotates  by 
some  'e  6  every  period,  is  possible  but  greatly  complicates  the  analysis  in  the 
gen’  e.  This  more  general  definition  is  important  to  consider  for  focusing  by 

a  It.  .nal  field;  in  that  case  one  can  find  matched  solutions  that  rotate  at  the 

Larm ...  (requency.  We  will  assume  here,  however,  that  the  beam  ellipse  returns 
to  its  initial  shape  and  orientation  every  period,  as  specified  by  Eq.  (11).  If  Eq. 
(7)  is  used,  the  matching  condition  for  any  distribution  of  the  form  f{Q),  where 
Q  =  vTWv,  is 

20=V020V0r,  (12) 

where  10  ~  2(j„)  and  V0  =  V(s0  +  L  |  s0 )■  The  solution  of  Eq.  (12)  for  Zo  gives  the 
matched  launching  condition  for  the  beam  matrix. 

The  eigenvalues  of  V0  all  have  complex  modulus  1  (since,  by  assumption,  the 
single-particle  motion  is  stable)  and  occur  in  conjugate  pairs.  When  the 
eigenvalues  are  all  distinct  the  eigenvectors  of  Vo  are  complete,  which  we  assume 
to  be  the  case.  If  we  define  a  matrix  U  that  has,  as  columns,  the  normalized 
eigenvectors  of  V„,  then  any  4x4  matrix,  in  particular  So,  may  be  written 

So  =  UDUT 


(13) 
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for  some  matrix  of  coefficients  D.  Substituting  Eq.  (13)  into  Eq.  (12)  and 
ordering  the  eigenvalues  of  V0  such  that  A,A2  =  A3A4  =  1,  one  obtains  the  result 
that  only  the  elements  Dn,  D2U  D 34,  and  Z)43  of  D  are  non-vanishing.  Requiring 
that  2o  be  symmetric  and  real  means  that  D  is  symmetric  and  real;  that  is,  the 
matched  solution  Z0  depends  on  just  two  arbitrary  real  parameters,  Dn  and  Du. 
In  the  decoupled  case,  of  course,  these  parameters  are  related  to  the  x  and  y 
emittances. 

For  s>  s0  the  matched  solution,  from  Eq.  (7),  is 

I0(s)  =  F(s)DFT(s).  (14) 


where  F(s)  =  V(s  \  sn)U  is  the  matrix  containing  the  four  Floquet  solutions  of  Eq. 
(1);  F  satisfies 

F(s0  +  L)  =  F(s0)  A,  (15) 

where  A  is  diagonal  and  contains  the  eigenvalues  of  V0. 

In  the  presence  of  the  K-V  space-charge  term,  Eq.  (8),  one  cannot  calculate  V0 
from  Eq.  (1)  before  knowing  the  matched  solution  for  2,  and  an  iterative  solution 
is  necessary;  the  most  direct  technique  is  to  neglect  space  charge  at  first,  calculate 
Vo  and  the  matched  solution  for  I,  then  recalculate  Vo.  and  so  on.  This  method  is 
used  successfully  in  the  example  treated  in  Section  4. 

Once  a  matched  solution  is  obtained  one  may  ask  whether  it  is  stable;  that  is,  if 
a  beam  is  launched  with  a  distribution  that  is  close  to  being  matched,  does  it 
remain  close  to  the  matched  solution  for  s  >  s0?  This  is  an  important  question  for 
highly  artificial  distributions  such  as  the  K-V  distribution  that  are  never  exactly 
(but  that  may  be  approximately)  realized  in  practice. 

In  the  absence  of  space  charge  the  matched  solution  is  clearly  stable  with 
respect  to  slight  variations  in  Z0;  replacing  In  by  Z0  +  <520  in  Eq.  (7)  clearly  does 
not  affect  the  stability  of  Z(s)  since  the  single-particle  solutions  V(s  |  s0)  are 
assumed  to  be  stable. 

In  the  presence  of  space  charge  one  examines  stability  by  linearizing  Eq.  (3). 
Writing  Z(s)  =  I„(s)  +  2,(s)  and  M(s)  =  M0(s)  +  Mt(s),  where  a  subscript  1 
denotes  the  perturbation,  one  obtains  the  linearized  equation  for  I,: 

Z',(s)  =  Mol,  +  (MoZ,)r  +  M,I  0  +  (M,Zo)r.  (16) 


Only  the  space-charge  part  of  M,  Eq.  (8),  contributes  to  M,  in  Eq.  (16).  The 
actual  linearization  of  MK  is  completely  straightforward,  if  a  bit  tedious.  The  nine 
quantities  needed  are 


1  3<hx  =  °yy  j 

q**  3axx  2SoSy 
1  dqxx  _  1  +  aJlSp 

q  XX  ^Oyy  Sy 

1  dqxx  -a xy  d 
qxt  daxy  SoSy  Iy 


(17a) 

(17b) 


(17c) 
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1  dqyy  1  +  (Jyy/lSo 

qyy  aoxx  "  s. 

(17d) 

• 

1  _  a**  J 

qyy  dOyy  2SoSt  yy 

(17e) 

1  d<tyy  _  ~a*y  d 

qyy  daXy  So 5, 

(170 

• 

— 

q,y  °oxx 

(17g) 

1  -  d 

_  -3  _  uyy 

<Ixy  V&yy 

1  1  d 
qxy  aOxy  Oxy 

where 

(17h) 

(17i) 

• 

d“  =  25§  +  <?“ 

(18a) 

“yy  *  2Sq  +  ^yy 

(18b) 

• 

d*y  =  ~F?  +  24xy 

0 

(18c) 

Equation  (16)  represents  ten  coupled  linear  equations  with  periodic 
coefficients;  the  ten  quantities  are  the  distinct  elements  of  the  4x4  symmetric 
matrix  By  reassembling  the  ten  quantities  into  a  column  vector  w,  Eq.  (16) 

may  be  formally  written 

• 

»v'(j)  =  N(j)h'(j), 

(19) 

where  /V(s)  is  a  10  x  10  matrix  with  period  L.  In  the  usual  way,  one  now  defines 
another  10  x  10  matrix  T(s  |  j0),  each  column  of  which  satisfies  Eq.  (19);  T  is 
defined  to  satisfy  the  initial  condition  T(s0 1  s0)  =  1,  the  identity  matrix.  The 
eigenvalues  of  T(s0  +  L\s0)  then  determine  the  stability  of  the  system;  an 
envelope  mode  is  (stable,  unstable)  if  the  modulus  of  the  corresponding 
eigenvalue  is  (less  than  or  equal  to,  seems  greater  than)  1 .  There  are  ten  eigenvalues 
of  T(s0  +  L  |  j0),  the  product  of  which  may  be  shown  to  be  1.  If  A  is  an  eigenvalue 
then  so  is  A*;  1/A  is  not  also  necessarily  an  eigenvalue,  as  it  is  in  the  decoupled 
case.9  The  ten  eigenvalues  are  accounted  for  as  follows:  In  the  decoupled  system, 
Eqs.  (3)  aind  (16)  both  reduce  to  two  2x2  matrix  equations,  one  each  for  x  and 
y.  Each  2x2  submatrix  of  £  is  symmetric  so  there  are  2x3  =  6  independent 
variables,  giving  6  eigenvalues  in  the  decoupled  case.  (Struckmeier  and  Reiser 
restrict  perturbations  to  those  that  do  not  change  the  emittances  e,  and  ey,  and  so 
find  6-2  =  4  eigenvalues.)  The  remaining  four  eigenvalues  are  attributable  to  the 
x  -  y  coupling. 
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4.  AN  EXAMPLE:  BEAM  TRANSPORT  IN  AN  1  =  2  STELLARATOR 
FIELD 


It  has  been  suggested  that  an  /  =  2  stellarator  field  may  have  some  advantages  for 
the  transport  of  high-current,  relativistic  electron  beams.5  Such  a  field  consists  of 
a  continuously  twisted  magnetic  quadrupoie  superimposed  on  a  constant  longitudi¬ 
nal  field.  The  analysis  of  the  form  of  the  matched  K-V  beam  in  the  continuous 
quadrupoie  case  has  been  treated  by  Gluckstem;6  inclusion  of  a  longitudinal  field 
requires  only  a  redefinition  of  certain  quantities  in  his  analysis.  Here  we  sketch 
out  the  solution  in  our  own  notation7  and  go  on  to  study  the  envelope  stability 
problem. 

The  quadrupoie  field  takes  the  form 

Bt(x,  y,  s)  ~  kBo(-x  sin  Iks  +  y  cos  2ks)  (20) 

By( x,  y,  s)  **  kBo(x  cos  2 ks  +  y  sin  2ks),  (21) 

where  k  and  &>  are  constants.  In  the  presence  of  the  space-charge  fields  Eqs. 
(A- 1,-2)  the  paraxial  equations  for  a  particle  near  the  axis  may  be  written 
compactly  as  .  . 

F  +  +  +  (22) 

where  §  =  x  +  i>,  b^  =  -eBJpc,  fi,  =  -kBoe/pc,  n,t  =  (cu6//3yc)2/2,  o)b  =  beam- 
plasma  frequency,  —e  and  p  are  the  electron  charge  and  momentum,  and  a,  b, 
and  a  are  the  radii  and  orientation  angle  of  the  ellipse  (see  Appendix). 

The  matched  beam  consists  of  an  ellipse  of  fixed  radii  that  rotates  with  the 
quadrupoie  field;  that  is,  a  and  b  are  constants,  and  a  =  ks.  The  substitution 
|  =  ye,ks  then  allows  the  following  solution  to  Eq.  (22): 

ip  —  A+e'*'*  +  o  y.A\e~iK^s  +  A~e‘*~J  +  o  ~A\e~iK~s , 
where  A  ±  are  arbitrary  complex  numbers, 


(23) 


+  M2 


1/2 


O*  = 


k2±  +  bzK±  +  k\ 


(24) 


(25) 


p2  k\  -  b2K±  +  k\ 

where  k\  =  k2  +  kbx  +  ns],  b2  =  b,  +  2k,  and  p2~  +  nri(a  ~  b)/(a  +  b ).  In  the 

interesting  special  case  where  b2  =  0,  corresponding  to  the  situation  in  which  the 
period  of  die  quadrupoie  field  equals  the  electron  cyclotron  wavelength,7  one 
finds  o ±  =  ±sgn  (p2). 

The  betatron  oscillations  are  stable  (*r±  are  real)  when 


b\{^b\-ktj  +  p\>0 


^b2-k22>0 
l*i|  -  iMa!  >0. 


(26a) 


Zob) 


(26c) 
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FIGURE  1  Stellarator  stability  plane.  The  quantities  u  and  u  are  defined  in  the  text  following  Eq. 
(26c). 


These  conditions  may  be  illustrated  using  the  auxiliary  quantities 

u  -  1  -  4*1/61,  v  *  \p2\lb\, 

as  shown  in  Fig.  1.  In  the  absence  of  space  charge,  only  Eq.  (26c)  represents  a 
nontrivial  constraint. 

We  construct  a  K-V  distribution  that  is  independent  of  the  phases  of  A±  : 

f  =fo6(f+  \A*\2  +/-  \A.\2  -  1],  (27) 

where  f0,  f±  are  constants.  The  amplitudes  may  be  expressed  as  functions  of  the 
real  and  imaginary  parts  of  rp: 


1 


l^4*|2  =  ^3 [*■-(  1  —  o_)y,-(l  +o_)V'IT  +  7^[*_(l  +  o_)V>,  +  (1  -  cr_)vr;i 


Di 


(28a) 


=  7^3 [#c^( 1  -  <7.)^-(l  +  a+)y,’]2  +  73[fc*(l  +  o~)v,  +(l  -  a+)y;)2, 

Is  2 


Di 


(28b) 


where  D,  =  -ir*(l  -  cr+)(l  +  ct_)  +  #r_(l  +  cr+)(l  -  o_),  and  D2  -  jc*( I  +  o„)  x 
(1  -  a.)  -  fr_(l  -  ct»)(  1  +  <7_).  One  may  alternately6  express  the  amplitudes  in 
terms  of  \pr,  <p,.  and  the  canonical  momenta  Pr  =  \p'r  -  2xp,  and  P,**\p',+  \b2xpr. 

Using  Eqs.  (28)  in  Eq.  (27)  and  integrating  over  xp'r  and  i pi  one  finds  that  the 
beam  radii  in  the  real  and  imaginary  directions  are 

(I  j.  »  \2 

(29a) 


2  (1  +  aJ*  ,  u  +  o-r 

a  — - ; - H  - 


A 


A 


b2  = 


(1  -o+f  .  (1-a-)2 


A 


A 


(29b) 
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Similarly,  one  may  calculate  the  emittances  in  the  real  and  imaginary  directions 
by  integrating  Eq.  (27)  over  (\pr,  rp'r)  and  (ip,,  ip!): 

‘ici(l  +  a+)2  ,  k2_(1  +  0_)2] 


A  +  A 
*d(l-o+)2  ,  *■-(!  ~  O-)2 


A 


A 


(30a) 


(30b) 


In  the  quadrupole-only  case,  Gluckstern  has  defined  the  emittances  as  areas  in 
the  Pr  -  ipr  and  P,  -  ip,  planes  instead.  For  the  present  problem  these  are  given  by 


e2cr  =  a 2  ( 

:H) 

(31a) 

eh  =  b2 1 

tbt)- 

(31a) 

where  c*  =  **(1  +  a*)  +  ib2(l  -  a*), 

c_  =  ir_(l  +  ct_)  +  £fi2(l  -  O-), 

<A  = 

ir+(l  -  a *)  +  ^fi2(l  +  a+),  and  d.  =  jc_(1  -  o_)  +  $h2(l  +  o'-). 

Using  the  definitions  of  Z  and  ip  and  Eqs.  (27)  and  (28),  the  value  of  the 
matched  Z  matrix  may  now  be  calculated.  The  nonzero  elements  at  s  =  0  are 

V  -1  2 

I  —  ~  a 


2|4  —  Z41  — 


1 


1  +  q„ 

~T~ 


[*(1  +  <7  +  )  +  K  +  (1  -  CT+)]  +  (+-»-) 


~ a+) +  **0  +  +  (+-*  “) 


z23  =  z12  =  - 


2„  =  r  b- 

4 


t  - 


O'* 


(32) 


A 


[/c(l  -  a*)  +  **(1  +  o+)j  +  (+-*-) 


Z44  =  —  [*(1  +  o„)  +  **(1  -  o+)j2  +  (  +  -*  -  ), 

J  + 

where  the  notation  +  (4 — *  - )  means  add  the  previous  term  with  all  +  subscripts 
replaced  by  -  The  matched  Z  matrix  is  dearly  divided  into  a  sum  of  two 
contributions,  one  from  the  fast  (+)  betatron  oscillation  mode  and  one  from  the 
slow  (-)  mode,  corresponding  to  the  decomposition  in  Eq.  (13).  The  determinant 
of  Z,  a  constant  of  motion,  is  given  by 

.  2 


(33) 


The  specification  of  Z  depends  on  the  two  unspecified  constants  f±,  which  may 
be  determined  in  a  variety  of  ways.  One  reasonable  choice  would  seem  to  be  to 
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specify  a  value  for  |2|  and  to  require  that  e,  =  e,  in  Eq.  (30).  This  choice, 
however,  gives  the  unsatisfactory  result  that  the  matched  beam  radii  are 
decreasing  functions  of  beam  current  unless  the  operating  point  is  in  the  "second 
stable  region”  of  Fig.  1,  where  u  >  1.  Curiously,  instead  of  specifying  e,  =  f,,  the 
choice  £gt  -  fc.»  from  Eq.  (31),  avoids  this  problem,  leading  to  beam  radii  that 
are  increasing  functions  of  current;  we  make  this  choice  in  the  numerical 
example,  below.  The  corresponding  condition  on  f±  is,  using  Eq.  (31), 

=  (34) 


except  in  the  special  case7  b2  -  0,  for  which  a±  =  ±  sgn  (/t2);  in  that  case  Eq.  (34) 
is  replaced  by  the  condition  ic_/+  =  x*/_. 

Figures  2a,  3a,  4a,  and  5a  illustrate  the  dependence  of  beam  radii  on  beam 
current,  beam  energy,  longitudinal  field  strength,  and  quadrupole  field  strength, 
respectively.  All  cases  have  k  =  -2n/l8  cm-1  and  (iy  |2|1/4  =  2.78  x  10'3rad-cm, 
which  correspond  to  typical  values  for  a  beam-transport  experiment  presently 
being  carried  out.11  Each  point  on  the  plots,  except  those  for  zero  current, 
requires  an  iterative  procedure  to  obtain  because  the  betatron  wavenumbers  k± 
depend  on  current  density,  that  is  on  the  beam  radii,  which  makes  Eqs.  (29) 
implicit  relations  for  the  beam  radii,  for  nonzero  current.  Once  the  beam  radii  are 
found,  the  elements  of  the  matched  2  matrix  at  s  —  0  follow  from  Eq.  (32). 

One  may  study  the  stability  of  the  matched  solution,  using  Eq.  (16).  One  might 
expect  that  since  the  matched  solution  is  a  constant  in  the  frame  rotating  with  the 
stellarator  field,  one  could  carry  through  the  stability  analysis  to  calculate 
oscillation  frequencies  and  growth  rates  analytically  for  this  example.  Actually, 
even  in  the  rotating  frame  [the  frame  of  the  ^-variable  of  Eq.  (23)]  the  two 
degrees  of  freedom,  ip,  and  ipit  are  coupled,  requiring  the  calculation  of  the 
eigenvectors  of  a  full  4x4  matrix  to  find  the  matched  solution.  Note  however 
that  in  the  special  case  where  b2  =  0,  described  following  Eq.  (25),  the  coupling  in 
the  rotating  frame  vanishes,  and  the  formalism  described  here  may  be  shown  to 
reproduce  the  K-V  envelope  equations,  with  constant  focusing  terms,  in  the 
rotating  frame;  in  this  case  the  envelope  oscillation  frequencies  are  easily 
calculated'*  analytically,  and  agree  with  those  obtained  numerically  from  Eq.  (16). 
To  study  the  linear  stability  problem  in  the  coupled  case,  Eq.  (16)  must  be 
integrated  simultaneously  with  Eq.  (3)  for  the  matched  solution,  using  Eq.  (32)  as 
me  initial  condition.  One  then  forms  the  matrix  T(LiO),  defined  following  Eq. 
(19)  and  finds  its  eigenvalues  A,-,  j  =  1,  2 . 10.  We  define  the  growth  rate  per 


period  as 


■ln|A,|. 


(35) 


Figures  2b-5b  show  plots  of  the  growth  rates  for  the  unstable  (I",  >  0)  envelope 
modes  as  functions  of  system  parameters.  One  sees  that,  depending  on  parameter 
values,  the  dominant  mode  has  a  fairly  large  peak  growth  rate,  which  suggests  the 
probability  of  large  emittance  growth  for  a  real  system  at  this  point.  The 
growth-rate  curves  of  Figs.  2b-5b  are  typical  of  an  instability  due  to  a  system 
resonance,  which  is  detuned  when  system  parameters  are  varied  slightly;  in  fact, 
this  is  an  example  of  a  "confluent  resonance”  for  the  stellarator  system.9  We  note 
that  calculation  of  linear  growth  rates  for  the  K-V  beam,  although  not 
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FIGURE  2  Equilibrium  beam  radii  and  growth  rate*  for  unstable  envelope  modes  versus  current  for 
*  1  MeV,  B  =  5  kG,  and  kB 0  =  500  G/ctn. 


40 


D.  CHERNIN 


BEAM  RADII 


ENCRCY  ( MeV) 


FIGURE  3  Equilibrium  b«am  radii  and  growth  rates  for  unstable  envelope  modes  versus  beam 
energy  for  =  300  A.  B,  =  5  kG,  and  fcflo  =  500G/cm. 


FIGURE  4  Equilibrium  beam  radii  and  growth  rates  for  unstable  envelope  modes  versus 
longitudinal  field  strength  for  =  300  A,  =  l  MeV,  and  kBo  »  500G/cm. 


QUAD  GRAD  (C/cm) 


FIGURE  5  Equilibrium  beam  radii  and  growth  rates  for  unstable  envelope  modes  versus 
quadrupole  field  strength  for  »  300  A,  -  1  MeV,  and  5,-5  IcG. 
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quantitatively  important  for  real  beams,  may  suggest  values  of  parameters 
(external  field  strengths,  beam  currents  and  energies,  etc.)  required  to  avoid 
significant  emittance  growth  in  real  beams.  Numerical  simulations  need  to  be 
done  to  test  (his  conjecture  for  the  stellarator  system,  though  simulations  in  the 
decoupled  case  show  that  emittance  growth  is  observed  in  simulations  of  beams 
with  realistic  profiles  when  the  growth  rates  for  corresponding  K-V  beams  are 
large.8,9 
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APPENDIX 

Space-Charge  Contribution  to  Matrix  M  for  K-V  Distribution 


We  consider  the  self-fields  of  an  elliptical  beam  of  axes  a  and  b;  axis  a  is  tilted 
at  an  angle  a  to  the  x-axis.  The  electro-  and  magneto-static  potentials  are  simply 
obtained  by  rotation  of  coordinates: 

<p{x,  y)  =  -  ^  ( qxxx 2  +  2 qxyxy  +  qyyy2) 

(A-l) 

A(x,  y)  =  0<ps, 

(A-2) 

where  A  is  the  line-charge  density  and 

4  /cos2  or  sin2  a\ 
q**~a  +  b\  a  +  b  ) 

(A-3a) 

4  /cos2  a  sin2  ar\ 
q,y~a  +  b\  b  +  a  ) 

(A-3b) 

4  .  /I  1\ 

qxy  =  .  sin  ar  cos  ar  1 

y  a  +  b  \a  b) 

(A-3c) 

The  values  of  the  spatial  ct’s  for  a  uniformly  populated  ellipse 

are 

<7„  =  ^  (a2  cos2  a  +  b2  sin2  a) 

(A-4a) 

oyv  =  7  (b2  cos2  a  +  a2  sin2  a) 

4 

(A-4b) 

cr  =  -  (a2  -  b2)  sin  a  cos  or. 

4 

(A-4c) 

It  is  now  simply  a  matter  of  algebra  to  eliminate  a,  bt  and  a 

from  Eqs.  (A-3)  in 
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favor  of  o„,  ayy,  and  axy,  using  Eqs.  (A-4).  The  intermediate  results 

ab  =  4 (a^a.y  -  (A-5) 

a2  +  b2  =  4(o„  +  oyy)  (A-6) 

are  useful  for  this  exercise.  The  results  for  qxx,  qyy,  and  qxy,  cited  in  the  text, 
follow. 
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APPENDIX  C 

Computer  Codes 


BBUSH: 

CBDRP: 


EMNAG: 

ORBITT: 

RPLANE: 

SBSC: 

SPIRAL: 

SSRS: 

STELMAT: 

STELTUNE: 

XYENV: 


Tracks  a  beam  pulse  through  an  arbitrary  series  of  gaps  and  transport  sec¬ 
tions,  calculating  the  transverse  deflection  due  to  the  beam  break-up  insta¬ 
bility. 

Solves  the  long  wavelength  dispersion  relation  for  stellatron  in  the  smooth  ap¬ 
proximation  and  plots  real  and  imaginary  parts  of  the  mode  frequencies  ver¬ 
sus  any  of  11  system  parameters  (beam  energy,  current,  quadrupole  strength, 
etc.). 

A  dispersion  relation  solver  for  the  Hughes-Godfrey  instability. 

Calculates  single  particle  orbits  in  the  field  produced  by  four  helically  wound 
wires,  superposed  on  a  uniform  longitudinal  field;  used  to  study  the  effect  of 
nonlinearities  on  single  particle  motion. 

Plots  curves  in  stellatron  parameter  space  on  which  certain  resonances  are 
excited. 

Solves  the  full  eigenvalue  problem  (no  smooth  approximation)  for  long  wave¬ 
length  beam  stability  in  the  stellatron. 

A  3-D  Biot-Savart  field  solver  and  single  particle  tracking  code;  used  for 
beam  matching  studies. 

Calculates  synchrotron  radiation  spectrum,  including  the  effect  of  betatron 
oscillations,  for  a  stellatron. 

Calculates  the  matched  launching  condition  for  stellarator  systems,  including 
space  charge;  also  calculates  growth  rates  of  unstable  modes  of  envelope 
oscillation. 

Calculates  stellatron  resonance  crossings. 

A  code  for  tracking  beam  envelopes  in  coupled  systems,  including  space 
charge  effects. 


Beam  stability  in  a  stellatron 

0.  Chemin 
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Collective  effects  in  the  /  =  2  stellatron  [  Phys.  Rev.  Lett.  50, 507  ( 1 983 )  ] ,  a  high-current  electron 
accelerator,  are  studied.  A  thin-beam  model  is  employed,  and  only  long-wavelength,  low- 
frequency  modes  are  considered.  The  eigenvalue  problem  is  formulated  in  general  and  solved 
analytically  in  the  Smooth  approximation;  comparison  with  a  numerical  solution  of  the 
eigenvalue  problem  is  excellent.  It  is  found  that  the  dispersion  relation  in  the  smooth 
approximation  is  identical  in  form  to  that  for  the  modified  betatron,  under  a  simple  substitution. 
An  analytical  expression  for  the  transition  energy,  obtained  earlier  by  a  simple  dynamical 
argument,  is  confirmed.  The  stellarator  field  is  found  to  reduce  the  growth  rate  of  the  negative 
mass  instability. 


I.  INTRODUCTION 

The  /  ss  2  stellatron,  a  betatron  with  an  added  l  =  2  stel¬ 
larator  winding,  was  originally  proposed1  in  order  to  im¬ 
prove  the  tolerance  of  a  betatron  or  modified  betatron2  to  a 
mismatch  between  the  average  beam  energy  and  the  vertical 
magnetic  field.  Such  a  tolerance  would  be  a  desirable  feature 
of  any  high-current  accelerator,  since  injectors  inevitably 
produce  a  beam  with  measurable  energy  spread  as  well  as  a 
variable  average  energy.  Though  the  original  analysis  in¬ 
cluded  space-charge  effects  in  a  simple,  non-self-consistent 
way,  the  argument  for  improved  momentum  compaction 
can  be  based  on  single  particle  orbit  considerations  alone. 

Treatment  of  space-cnarge  effects  in  the  stellatron  is 
complicated  by  the  lack  of  symmetry  in  the  fields.  Self-con¬ 
sistent,  monoenergetic,  KV-like3  equilibria  can  be  construct¬ 
ed,4  however,  and  these  give  explicit  relations  among  the 
beam  radii  (the  beam  is  elliptical,  in  general),  emittance, 
current,  energy,  and  externally  applied  fields;  the  shifts  in 
the  betatron  frequencies  due  to  space  charge  are  also  ob¬ 
tained  from  these  equilibrium  studies. 

In  the  present  paper  the  stability  of  a  beam  confined  in  a 
stellatron  is  considered.  This  work  was  motivated,  in  part,  by 
certain  experimental  results  on  the  /  =  2  stellatron  at  the 
University  of  California  at  Irvine.5  Though  beam  currents 
have  been  increased  in  that  experiment  by  the  addition  of  the 
1  =  2  windings,  currents  are  still  not  as  large  as  one  might 
expect  from  the  amount  of  injected  charge.  An  instability  of 
some  kind  may  be  responsible  for  the  current  loss.  An  at¬ 
tempt  to  discover  possible  current  limiting  mechanisms  in 
the  stellatron  initiated  the  work  presented  here. 

The  present  stability  calculation  makes  the  following 
approximations:  All  first-order  fields  are  treated  in  the  long- 
wavelength  U>a3minor  radius  of  toroidal  chamber), 
low-frequency  ( a<c/a )  limit;  toroidal  corrections  to  the 
first-order  fields  are  neglected.  We  note  in  particular  that  the 
long-wavelength  approximation  may  omit  the  potentially  in¬ 
teresting  caseA —lirry/m,  where  r0  is  the  major  radius  of  the 
torus  and  m  is  the  number  of  stellarator  field  periods  around 
the  device.  The  beam  is  treated  as  a  structureless,  thin, 
threadlike  charge  distribution  so  that  the  effects  of  betatron 
oscillations  and  energy  spread  on  the  equilibrium  are  ig¬ 
nored;  the  magnitude  of  these  effects  can  be  roughly  estimat¬ 
ed,  however.6  The  beam  environment,  which  has  a  critical 


effect  on  beam  stability,  is  treated  simply:  the  beam  is  as¬ 
sumed  to  be  contained  in  a  perfectly  conducting  toroidal 
chamber.  Only  the  longitudinal,  or  negative  mass  mode  is 
unstable  in  these  circumstances.  Several  comments  are  made. 
below,  however,  regarding  the  effects  of  walls  of  arbitrary 
impedance  on  both  longitudinal  and  transverse  modes. 

Under  the  above  approximations  we  have  the  following 
findings  to  present:  (1)  the  eigenvalue  problem  governing 
the  longitudinal  and  transverse  modes  of  a  beam  in  a  stella¬ 
tron;  (2)  a  dispersion  relation,  obtained  from  an  approxi¬ 
mate  solution  to  the  eigenvalue  problem,  which  is  identical 
in  form  to  that  obtained  for  the  modified  betatron,7  with  the 
redefinition  of  a  certain  symbol;  ( 3 )  confirmation  of  an  ana¬ 
lytical  expression,*  obtained  by  a  simple  dynamical  argu¬ 
ment,  for  the  so-called  transition  energy,  below  which  the 
longitudinal  mode  is  stable  (in  the  absence  of  dissipative 
effects);  (4)  a  quite  favorable  numerical  comparison 
between  the  solution  to  the  eigenvalue  problem  of  ( 1 )  and 
the  dispersion  relation  of  (2)  for  the  longitudinal  mode. 

This  paper  is  organized  as  follows.  Section  II  describes 
the  derivation  of  the  linearized  equations  of  motion,  assum¬ 
ing  the  applied  fields  are  given  near  the  axis  of  the  torus. 
Section  III  outlines  the  derivation  of  the  first-order  charge 
and  current  densities  in  terms  of  the  single  particle  variables 
using  a  summation  over  initial  conditions.  In  Sec.  IV  the 
approximate  solution  to  Maxwell’s  equations  is  obtained, 
and  the  total  first-order  fields  to  be  used  in  the  equations  of 
motion  are  derived.  In  Sec.  V  the  eigenvalue  problem  is  for¬ 
mulated  and  is  shown  to  reproduce  the  modified  betatron 
dispersion  relation  in  the  absence  of  the  twisted  quadrupole 
field.  An  approximate  solution  based  on  a  two-time  scale 
analysis  is  presented  and  some  comments  are  made  on  the 
form  of  the  resulting  dispersion  relation.  Finally,  in  Sec.  VI, 
a  numerical  solution  to  the  full  eigenvalue  problem  is  pre¬ 
sented  and  compared  to  results  from  the  dispersion  relation 
of  Sec.  V. 

Section  VII  summarizes  these  results  and  makes  a  few 
additional  remarks. 

II.  EQUATIONS  OF  MOTION 

We  shall  use  a  standard  (r,  &,  z)  cylindrical  coordinate 
system.  In  the  unperturbed  state  we  have  a  thin  beam  en¬ 
closed  in  a  perfectly  conducting  torus  of  major  radius  r0  and 


see 
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